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ABSTRACT 

Design  studies  are  performed  and  problem  areas 
analyzed  with  respect  to  the  redesign  of  the  MIT  Reactor 
in  order  to  accommodate  a  compact  core  configuration. 

Mechanical  design  studies  are  performed  for  a  DpO 
moderated,  cooled,  and  reflected  compact  core  design. 
Designs  are  proposed  for  the  core  tank  and  reactor  top 
shielding  configuration,  the  core  support  mechanism,  the 
control  rod  drive  mechanisms,  and  the  refueling  mechanism. 

The  necessary  composition  and  thickness  of  reactor  top 
shielding  are  determined  and  shielding  calculations  per- 
formed for  both  a  D20  cooled  core  and  an  H20  cooled  core. 
Shielding  above  the  D20  core  consists  of  an  8  ft.  column  of 
DoO  plus  a  shielding  enclosure  consisting  of  6  in.  of  mason- 
ite,  1/2  in.  of  boral,  and  1  ft.  of  lead,  through  which  a 
maximum  dose  rate  of  0.2  mrem/hr  is  calculated.   Shielding 
above  the  H2O  cooled  core  consists  of  11  ft.  of  H20  topped 
by  a  10  in.  lead  cover,  through  which  a  maximum  dose  rate 
of  0.9  mrem/hr  is  calculated. 

The  fluid  flow  and  heat  transfer  problems  associated 
with  the  design  of  a  compact  core  configuration  are  analy- 
zed and  it  is  determined  that  the  currently  installed  MITR 
process  system  is  not  capable  of  providing  adequate  heat 
removal  from  the  compact  core,  if  smooth  flat  plate  fuel 
elements  are  utilized.   An  experimental  investigation  is 
undertaken  to  determine  the  effectiveness  of  using  fins  as 
a  means  of  augmenting  fuel  plate  heat  transfer.   Heat 
transfer  from  a  smooth  rectangular  channel  2.5"  x  0.090"  x 
25-1/2"  is  first  investigated  and  experimental  results 
confirm  the  validity  of  using  the  modified  Colburn  corre- 
lation, Nuf  =  0.023  Ref°*°  Prf°'3,  to  predict  smooth 
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channel  heat  transfer  values.   Two  finned  channel  test 
sections  are  investigated,  one  with  fins  cut  into  the 
channel  transverse  to  the  coolant  flow  direction,  and  a 
second  with  fins  milled  parallel  to  the  coolant  flow 
direction.   In  each  case  the  fins  are  milled  to  be  20 
mils  high,  20  mils  wide,  and  spread  40  mils  apart.   The 
transverse  fins  are  found  to  augment  heat  transfer  by  a 
factor  of  2.2,  while  a  factor  increase  of  1.6  is  obtained 
from  the  longitudinal  fins.   However,  due  to  safety, 
pressure  drop  and  cost  considerations,  longitudinally 
finned  fuel  plates  are  recommended  for  use  in  the  High 
Flux  MIT  Reactor. 

Supervisor:   Dr.  Theos  J.  Thompson 

Title:   Professor  of  Nuclear  Engineering 
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INTRODUCTION 
Since  its  completion  in  the  Spring  of  1958,  the  MIT 
Reactor  has  served  as  a  vital  research  tool  for  the  MIT 
and  surrounding  scientific  community.   At  the  time  the 
MIT  Reactor  was  designed,  little  was  known  of  compact, 
heavy  water  moderated  cores.   Consequently,  the  MIT  Reactor 
was  designed  with  fuel  element  spacing  sufficient  to  pro- 
vide for  a  high  percentage  of  neutron  moderation  to  take 
place  within  the  core  area.   Results  published  in  1964 
upon  the  completion  of  the  Brookhaven  High  Flux  Beam 
Research  Reactor  provided  the  experimental  information 
needed  on  compact,  undermoderated  core  configurations  in 
order  to  accomplish  an  exaggerated  thermal  flux  peak  in 
the  reflector.   By  incorporating  a  similar  compact  core 
design  into  the  MIT  Reactor,  it  is  reasonable  to  assume 
that  the  thermal  flux  coming  out  of  the  beam  ports  could 
be  increased  by  at  least  a  factor  of  three,  and  possibly 
as  much  as  a  factor  of  four  or  five.   Such  an  increase 
would  markedly  increase  the  value  of  the  MIT  Reactor  to 
research.   In  addition,  if  the  beam  ports  were  to  be 
located  below  the  base  of  the  core,  the  core  gamma  ray 
and  fast  neutron  background  would  be  markedly  reduced  in 
intensity  due  to  the  fact  that  any  core  gammas  and  fast 
neutrons  which  leak  out  must  have  been  scattered  at 
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least  once  previous  to  escaping.   Such  a  design  would 

markedly  increase  the  quality  of  the  thermal  neutron 

beam  emitted  from  the  beam  ports. 

It  is  primarily  for  the  above  reasons  that  a  MITR 
redesign  group  was  formed  under  the  direction  of  Dr. 
T.J.  Thompson.   The  design  group  is  composed  of  both 
graduate  students  and  MITR  staff  personnel.   Besides  Dr. 
Thompson  and  the  author  of  this  thesis,  the  members  of 
the  design  group  who  have  been  active  throughout  the  re- 
design process  are  graduate  students:   Andrews  Addae 
working  on  core  physics,  Dan  Kennedy  working  on  beam  port 
optimization,  and  Bob  Sanders  working  on  a  cold  neutron 
facility  for  the  thermal  column.   Mr.  Ed  Barnett,  working 
on  all  engineering  phases  of  the  redesign,  represents  the 
MITR  staff. 

Although  never  stated  formally  in  writing,  the  design 
criteria  outlined  orally  by  Dr.  Thompson  at  the  outset  of 
the  design  were  the  following: 

1)  The  thermal  neutron  flux  at  the  beam  ports  must 
be  increased  by  a  factor  of  3-5. 

2)  The  redesigned  reactor  must  be  as  safe,  or 
safer,  than  the  existing  reactor. 

3)  Simplicity  of  design  must  be  a  prerequisite  for 
the  design  of  any  component. 


16 

4)  There  must  be  no  buried  mechanisms. 

5)  There  must  always  be  a  means  for  retreat.  To 
quote  Dr.  Thompson,  "We  have  a  going  reactor,  we  don't 
want  to  replace  it  with  a  monument." 

6)  The  redesigned  reactor  must  be  reliable. 

7)  The  cost  of  the  redesign  project  must  be  within 
$500,000. 

The  redesign  project  has  actually  been  divided  into 
two  distinct  design  phases.   The  first  phase  began  in 
earnest  in  the  Fall  of  1967  and  was  centered  around  the 
design  of  a  D^O  moderated,  cooled  and  reflected  reactor. 
Phase  I  concluded  in  the  Fall  of  1968  and  was  succeeded 
by  an  intermediate  phase  in  which  a  mixture  of  H^O  and 
DpO  was  considered  as  the  moderator  and  coolant  for  the 
core,  with  pure  D^O  as  the  reflector.   In  December  of 
1968  phase  II  began.   Phase  II  is  centered  around  the 
design  of  a  light  water  cooled  and  moderated  core  surrounded 
by  a  heavy  water  reflector. 

At  the  outset  of  phase  I,  it  was  apparent  that  there 
would  be  many  challenges  in  the  area  of  mechanical  design. 
The  incorporation  of  a  compact  core  configuration  dictated 
that  new  methods  be  devised  to  support  not  only  the  fuel 
elements,  but  the  control  elements  as  well.   A  change  in 
core  tank  design  seemed  to  be  required.   This  author's 
first  responsibility  as  a  member  of  the  redesign  team  was, 
then,  to  attempt  to  find  solutions  to  some  of  the  many 
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design  problems  confronting  the  redesign  project.   The 
work  done  toward  this  objective  is  summarized  in  Chapter 
2  of  this  thesis. 

As  the  redesign  project  continued  and  a  necked  cool- 
ant tank  was  adopted,  concern  for  the  requirements  of 
reactor  top  shielding  emerged.   This  author's  second  area 
of  responsibility  became,  therefore,  that  of  designing  the 
shielding  to  be  required  for  the  reactor  top,  and  that  of 
performing  shielding  calculations  to  justify  the  proposed 
shield  composition.   The  results  of  work  done  in  this  area 
are  presented  in  Chapter  3  of  this  thesis. 

As  stated  at  the  beginning  of  this  introduction,  one 
of  the  primary  purposes  for  redesigning  the  MITR  was  to 
accommodate  a  compact  core  configuration  in  order  to  obtain 
a  large  thermal  flux  peak  in  the  reflector.   However, 
almost  by  definition  the  achievement  of  a  pronounced  ther- 
mal flux  peak  in  the  reflector  will  be  accompanied  by  a 
sharply  increasing  thermal  flux  gradient  across  the  outer- 
most fuel  plates.   The  outermost  fuel  plates  will  conse- 
quently generate  a  much  higher  power  density  than  the 
inner  fuel  plates  in  the  core.   From  the  results  presented 
in  reference  (1),  it  is  reasonable  to  suspect  that  due  to 
the  sharply  rising  thermal  flux  gradient  in  the  area  of 
the  outermost  fuel  plates,  the  reactivity  coefficient  due 
to  coolant  voiding  in  the  channels  adjacent  to  these  fuel 
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plates  may  be  positive.   It,  therefore,  becomes  imperative 
that  boiling  not  be  allowed  to  occur  in  the  area  of  the 
outermost  fuel  plates,  which  is  equivalent  to  saying  that 
boiling  must  not  occur  anywhere  in  the  core  since  the  : 
outermost  fuel  plates  generate  the  most  power.   The  maximum 
rate  at  which  heat  can  be  removed  from  the  outermost  fuel 
plates  without  allowing  nucleate  boiling  to  occur  thereby 
sets  an  upper  limit  to  the  power  density  allowable  in  the 
outermost  fuel  plates,  which  in  turn  influence  the  maximum 
thermal  flux  peak  which  can  be  obtained  in  the  reflector. 
It  is  this  critical  area  of  core  heat  removal  which  became 
this  author's  third  and  final  area  of  responsibility.   The 
results  of  work  done  in  this  area  are  presented  in  Chapters 
4  and  5  of  this  thesis. 
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Chapter  1 
DESCRIPTION  OF  EXISTING  MITR  COMPONENTS 

1.1  General 

Before  proceeding  with  a  description  of  the  design 
and  calculational  work  which  was  performed  as  a  part  of 
the  High  Flux  MITR  Project,  it  seems  appropriate  to  first 
present  a  brief  description  of  those  existing  MITR  compo- 
nents and  systems  which  are  the  subject  of  redesign  work 
presented  in  this  thesis.   The  following  information  was 
taken  from  the  operating  manual  for  the  MIT  Reactor  (2). 
The  operating  manual  should  be  consulted  if  more  detailed 
descriptions  are  desired. 

1.2  Reactor  Coolant  Tank 

The  reactor  tank  is  48  inches  in  inside  diameter  and 
has  a  wall  thickness  of  1/4  inch,  except  at  the  bottom 
where  it  is  3/8  inch  thick.   The  overall  height  from  the 
dished  bottom  shell  to  the  top  flange  is  7  feet  5-3/8 
inches.   The  reactor  tank  is  made  of  special  high  purity 
aluminum  and  has  a  design  pressure  of  40  psig.   During 
normal  operation,  however,  the  pressure  in  the  tank  above 
the  coolant  is  only  1  inch  of  H20  above  atmospheric  pres- 
sure.  The  tank  is  shown  in  the  vertical  section  through 
the  reactor  of  Figure  1.2.1,  and  is  also  shown  in  Figure 
1.2.2. 
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21 


BACK-UP  RING- 


THERMAL     SHIELD 


UPPER    ANNULAR    RING 


SUPPORT     BLOCKS 


LOWER    ANNULAR     RING 


REACTOR   CORE  TANK    SUPPORT 
FIGURE    1.2.2 


22 
Heavy  water  flows  into  the  reactor  tank  through  an 

8  inch  diameter  aluminum  pipe,  through  the  holes  in  the 

i 
side  of  the  inlet  pipe  cap  and  out  into  the  entrance 

plenum.  It  then  flows  up  through  the  fuel  elements  and 
out  into  the  tank  itself.  The  D20  from  all  the  fuel  ele- 
ments flows  out  of  the  tank  through  a  transition  section 
which  extends  a  short  distance  above  the  plenum  head,  and 
then  through  the  8  inch  aluminum  outlet  pipe.  The  normal 
total  coolant  flow  rate  is  approximately  2000  gallons  per 
minute. 

The  reactor  tank  is  attached  to  the  lower  shielding 
plug  by  24  one  inch  stainless  steel  socket  head  cap  screws. 
The  cap  screws  compress  a  hollow  stainless  steel  0  ring 
between  the  tank  and  the  shield  plug  to  provide  a  helium 
seal  designed  for  80  psig.   Consequently,  if  a  high  pres- 
sure is  developed  in  the  tank,  the  tank  will  rupture 
downward  before  the  shield  plug  gives  way. 
1. 3  Shielding  Above  the  Core 

1.3.1  Lower  Annular  Ring 

The  lower  annular  ring  (shown  in  Figures  1.2.1  and 
1.2.2)  is  a  concrete-filled  steel  structure  which  is  56 
in.  in  inside  diameter,  approximately  8  ft.  in  outside 
diameter,  and  4  ft.  high.   Its  lower  surfaces  are  covered 
with  1/4  in.  of  boral.   A  4  in.  thick,  water-cooled  lead 
thermal  shield  lines  the  Inner  surface  over  the  lower  18 
in.   There  are  conduits  for  thermocouple  leads  and  helium 
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gas  pipes,  and  6  sleeves  are  provided  for  thimbles  which 
slant  into  the  graphite  reflector. 

1.3.2  Lower  Shield  Plug 

The  lower  shield  plug  is  a  large  cylindrical  struc- 
ture, the  bottom  of  which  is  lined  with  a  1/4  in.  sheet 
of  boral.   Above  the  boral  is  a  1-1/8  in.  steel  plate,  4 
in.  of  lead  provided  with  cooling  coils,  29  in.  of  concrete, 
a  steel  cover  plate,  a  6  in.  space  filled  with  paraffin  and 
steel  shot,  and  a  1-1/8  in.  stainless  steel  top  plate.   The 
plug  is  penetrated  by  numerous  sleeves  for  the  concrete- 
filled  steel  shielding  plugs  which  are  part  of  the  fuel  ele- 
ment and  control  rod  assemblies. 

1.3.3  Upper  Annular  Ring 

The  upper  annular  ring,  together  with  the  rotary  plug, 
complete  the  biological  shielding  above  the  core  tank.   The 
three  foot  thick  annulus  is  made  of  dense  concrete  encased 
in  steel.   Plugged  holes  are  provided  for  access  to  the 
vertical  thimbals  which  protrude  through  the  lower  annular 
ring  and  into  the  graphite  region.   The  upper  annulus  is 
spaced  9-1/4  in.  from  the  lower  annulus  by  several  support- 
ing blocks  bolted  to  the  lower  annular  ring  flanges. 

1.3.4  Rotary  Plug 

The  heavy  concrete  rotary  plug  fits  into  the  center 
of  the  upper  annular  ring  and  rests  on  a  ball-bearing 
thrust  ring.   It  is  rotated  by  a  drive  motor  and  friction 
belt.   When  the  plug  is  latched  in  its  operating  position, 
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the  control  rod  racks  can  move  upward  into  cylindrical 
slots  in  the  bottom  of  the  plug.   To  remove  a  fuel  ele- 
ment, the  control  rods  must  be  fully  inserted  and  the 
plug  unlatched  and  rotated  to  the  appropriate  indexed 
position.   The  shielding  plug  directly  over  the  fuel  ele- 
ment position  can  then  be  removed.   The  entire  rotary 
plug  must  be  removed  for  access  to  the  control  rod  drives. 
1.4  Fuel  Element  Assembly 

A  fuel  element  assembly  is  shown  in  Figure  1.4.1  and 
consists  of  a  fuel  section  supported  by  a  lower  adapter, 
an  upper  adapter  and  a  shield  plug.   The  30  positions  into 
which  fuel  assemblies  can  be  inserted  are  indicated  by 
numbers  1  through  30  in  Figure  1.4.2.   There  is  a  central 
element  position  surrounded  by  an  inner  circle  of  6  posi- 
tions, all  spaced  6.375  in.  center-to-center.   These  in 
turn  are  surrounded  by  an  outer  circle  of  12  positions 
spaced  6.858  in.  center-to-center.   The  11  additional 
positions  are  located  around  the  outside  edge  of  the  reac- 
tor tank. 

The  fuel  section  consists  of  18  curved  plates  assem- 
bled between  two  grooved  side  plates  which  are  O.187  in. 
thick  and  3.168  in.  wide.   Two  nonfuel-bearing  outer 
plates  are  0.060  in.  thick  and  are  curved  to  a  5-1/2  in. 
radius.   They  are  spaced  0.117  in.  apart  to  form  17  cool- 
ing water  passages. 
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The  fuel  assembly  shielding  plug  fits  into  a  stepped 
steel  sleeve  embedded  in  the  lower  shield  plug.   The 
weight  of  the  assembly  rests  on  the  step  in  the  lower 
shield  plug.   This  provides  a  convenient  means  of  removing 
the  fuel  assembly  through  the  tank  top,  and  when  in  place 
the  shielding  properties  are  equivalent  to  that  of  the 
lower  shield  plug.   The  shield  plug  is  securely  latched 
into  place  by  two  small  lugs  which  slide  outward  and 
engage  a  groove  in  the  sleeve.   The  lugs  are  moved  by  turn- 
ing a  square  projection  on  a  cam  which  bears  against  the 
lug. 
1.5  Control  Rod  Mechanisms 

The  MITR  has  one  regulating  rod  and  six  shim  safey 
rods  arranged  in  the  core  as  shown  by  the  solid  black 
circles  in  Figure  1.4.2.   The  shim  safety  rod  assembly  is 
shown  in  Figure  1.5.1  and  is  described  below. 

The  absorber  section  is  a  hollow,  double-walled, 
aluminum  cylinder,  28  in.  long  and  2-1/4  in.  in  outside 
diameter,  in  which  is  sandwiched  a  0.040  in.  thick  cadmium 
tube  2.15  in.  in  outside  diameter  and  26-1/8  in  long.   The 
absorber  cylinder  is  connected  to  the  magnet  armature  (a 
nickel-plated  soft  iron  disk)  and  dashpot  by  means  of  a 
spider-like  section,  thus  permitting  D20  to  flow  freely 
through  the  top  of  the  cylinder  as  it  drops.   During  a 
reactor  scram  as  the  absorber  nears  its  fully  inserted 
position,  it  is  decelerated  by  a  cushion  of  heavy  water 
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which  leaks  off  between  a  washer  and  the  dashpot  case 
with  the  assistance  of  coil  springs  in  the  dashpot.   A 
1/4  in.  thick  disk  of  cadmium  on  the  bottom  of  the  dashpot 
shields  the  iron  armature  from  neutrons.- 

The  electromagnet  consists  of  a  silicon-insulated 
copper  wire,  wound  onto  a  transite  bobbin.   It  is  fitted 
into  an  Armco  magnet-iron  housing  which  is  enclosed  by  two 
stainless  steel  end  plates.   The  magnet  wires  are  connected 
to  a  threaded  connector,  which,  when  coupled  with  the  racks 
forms  a  hei  letically  sealed  unit.   Magnet  wires  are  led 
through  a  tube  in  the  rack  rod,  then  down  the  outside  to 
a  take-up  real  mounted  at  the  top  of  the  casing  assembly. 
The  rack  is  guided  by  the  casing,  and  the  escape  of  helium 
or  DpO  around  the  shaft  is  prevented  by  Styrene  chevron 
packing. 

The  rods  are  driven  by  a  rack  and  pinion  arrangement. 
The  drive  shaft  from  the  pinion  gear  extends  to  the  motor 
and  gear  box  assembly  at  the  utility  shelf  (see  Figure 
1.5.2).   The  six  drive  shafts  are  arranged  radially  but 
differ  in  length  and  arrangement  in  order  to  clear  shutters 
and  experimental  facilities.   All  drive  shafts  are  enclosed 

in  sleeves  from  the  utility  shelf  to  the  between-plug  space, 

Li 

and  are  sealed  to  limit  the  amount  of  A   entering  the  air 

stream. 
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Chapter  2 
PRELIMINARY  MECHANICAL  DESIGN  STUDIES  DONE  FOR  PHASE  I 

2. 1  General 

In  the  early  stages  of  the  MIT  reactor  design  project 
it  was  this  author's  primary  responsibility  to  attempt  to 
find  solutions  to  some  of  the  many  design  problems  criti- 
cal to  the  construction  of  a  High  Flux  MITR.   In  so  doing 
many  design  ideas  were  proposed,  few  of  which  gained 
acceptance  by  the  entire  design  group.   It  is  the  purpose 
of  this  chapter  to  relate  some  of  the  design  proposals 
which  can  be  primarily  attributed  to  this  author,  though 
at  times  when  working  with  a  design  group  it  is  admittedly 
difficult  to  distinguish  which  ideas  are  completely  one's 
own  from  those  which  may  have  had  their  roots  in,  for 
example,  a  comment  made  at  a  design  group  meeting.   It  is 
felt  that  although  a  particular  design  scheme  may  not  have 
appeared  best  for  the  MITR  redesign,  it  may  at  some  future 
date  prove  useful  to  another  reactor  design  project,  and 
therefore  is  worthy  of  description  in  this  chapter.   In 
addition,  the  design  proposals  serve  as  a  kind  of  chronology 
of  mechanical  design  progress. 

2.2  Core  Tank  and  Reactor  Top  Shielding  Configuration 
The  basic  necked  core  tank  design  as  finally  adopted 

evolved  from  the  compact  core  concept  and  short  tank  pro- 
posed in  the  original  redesign  sketches  of  Prof.  T.J.  Thompson 
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and  first  formally  reported  in  August  1966  in  a  22.90 
project  report  by  D.M.  Goebel  (3) .    The  original  purpose 
in  designing  a  necked  core  tank  was  to  provide  for  adequate 
shielding  of  the  core  during  refueling  operations.   (Prior 
to  the  introduction  of  an  enclosed  refueling  structure 
atop  the  reactor,  it  was  proposed  that  there  be  a  single 
refueling  plug  which  would  be  removed  during  refueling 
operations.   Without  approximately  11.5  feet  of  D^O  separa- 
ting the  core  from  the  reactor  top,  the  radiation  level 
would  be  prohibitive  to  workers  in  the  area  of  the  open 
refueling  plug.)   In  order  to  accompolish  the  objective  of 
providing  adequate  shielding  during  refueling  operations, 
but  at  the  same  time  to  minimize  the  required  DpO  inventory, 
the  concept  of  two  rotating  shield  plugs  was  conceived  and 
originally  proposed  by  this  author.   This  design  concept, 
outlined  in  Figure  2.2.1  would  utilize  a  1-1/2  ft.  thick 
lower  shield  plug  fabricated  of  lead  which  would  rotate 
around  a  track  of  aluminum  rollers.   Releasable  latches 
would  prevent  the  shield  plug  from  beint  "shot"  upward  by 
a  sudden  increase  in  reactor  pressure.   Inserted  in  the 
lower  shield  plug  would  be  three  6  in.  diameter  removable 
lead  plugs,  so  situated  that  by  rotating  the  lower  shield 
plug  one  of  the  smaller  plugs  would  be  caused  to  line  up 
atop  any  desired  fuel  element.   A  tool  inserted  from  the 
reactor  top  through  the  upper  shield  plug  could  be  used  to 


FIGURE       2.2. 
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unlatch  a  6  in.  diameter  lead  plug,  lift  it  out  of  its 
slot,  and  set  it  aside.   The  upper  rotating  shield  plug 
would  consist  of  3  ft.  of  dense  concrete,  and  would  rotate 
on  stainless  steel  rollers  and  contain  three  6  in.  diameter 
removable  plugs  similar  to  those  of  the  lower  shield  plug. 
An  inflatable  boot  would  be  located  around  the  upper  shield 
which  could  be  deflated  to  allow  the  shield  plug  to  rotate 
and  inflated  to  act  as  a  helium  seal.   (An  inflatable  boot 
could  be  installed  around  the  lower  shield  plug  if  the 
upper  shield  plug  should  be  entirely  removed  in  order  to 
perform  maintenance  on  the  control  rod  drive  mechanisms, 
for  example.)   Both  the  upper  and  lower  shield  plugs  would 
be  cooled  by  duplicate  D^O  cooling  coils  cast  into  the 
shield  plugs.   When  rotated  to  an  indexed  position  for 
reactor  operation,  the  coolant  lines  in  the  shield  plugs 
would  be  caused  to  line  up  with  their  respective  supply 
lines  in  the  annular  rings. 

The  refueling  operation  would  proceed  somewhat  as 
follows :  One  of  three  6  in.  diameter  inserts  in  the  lower 
shield  plug  would  be  unlatched  and  raised  by  a  tool  inser- 
ted from  above.   The  lower  shield  plug  would  then  be  rota- 
ted to  an  indexed  position  aligning  the  opening  in  the 
lower  shield  plug  over  the  fuel  element  to  be  removed. 
The  6  in.  diameter  plug  could  then  be  rested  on  the  lower 
shield  plug  for  convenience.   Next  the  appropriate  6  in. 
diameter  insert  from  the  upper  shield  plug  would  be  removed. 
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At  this  time,  the  insert  would  not  be  aligned  with  the 
opening  in  the  lower  shield  plug.   The  refueling  flask 
would  be  brought  into  position  on  the  reactor  top,  the 
upper  shield  plug  rotated  to  an  indexed  position  so  as  to 
align  its  opening  above  the  fuel  element  to  be  removed, 
and  the  designated  fuel  element  raised  up  into  the  refuel- 
ing flask.   Once  in  the  refueling  flask,  the  fuel  element 
could  either  be  flipped  end  for  end  and  reinserted  into 
the  reactor,  or  set  aside  while  a  new  fuel  element  is 
placed  into  the  core.   The  above  procedure  would  be 
accomplished  in  reverse  to  complete  the  refueling  operation, 

By  utilizing  the  concept  of  a  double,  rotating  shield 
plug,  the  reactor  DpO  tank  would  remain  shaped  similar  to 
the  present  DpO  tank  (see  Figure  1.4).   Two  important 
differences  would  be  apparent,  however.  Although  4  ft.  in 
diameter,  the  DpO  tank  would  be  shortened  from  7-1/2  ft. 
to  6-1/2  ft.  in  height.   The  volume  of  DpO  in  the  system 
would  thereby  be  decreased  by  14$.  Contained  within  the 
4  ft.  diameter  outer  tank  would  be  an  inner  tank  of  diame- 
ter just  sufficient  to  contain  the  core.   This  inner  tank 
would  extend  the  full  length  of  the  outer  DpO  tank  and 
would  be  constructed  of  sufficient  thickness  of  aluminum  to 
be  able  to  contain  a  sudden  DpO  vaporization  resulting  in 
40  psi  overpressure  in  the  core.   As  a  consequence,  the  4 
ft.  diameter  DpO  tank  need  only  be  of  sufficient  strength 
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to  support  the  coolant  contained  within  it.   There  are  two 
major  reasons  for  utilizing  an  inner  tank  to  separate  the 
core  region  from  the  reflector  region.   First  is  the 
possibility  of  maintaining  the  reflector  DpO  at  a  lower 
temperature  than  the  coolant  DpO,  and  thus  make  the  reflec- 
tor DpO  a  better  moderator.   Second  is  that  by  separating 
the  reflector  from  the  core,  the  entire  reflector  region 
could  be  rapidly  dumped  of  its  DpO,  thus  creating  a  very 
important  mechanism  for  reactor  safety. 

Both  the  upper  and  lower  annular  rings  would  remain 
the  same  as  in  the  present  MITR,  and  the  core  tank  likewise 
would  be  supported  as  it  is  now.   Such  mechanical  similari- 
ties with  the  existing  reactor  would  tend  to  reduce  the 
cost  of  converting  to  the  high  flux  MITR  core  configuration. 

The  control  rod  mechanisms  for  the  above  described 
reactor  tank  configuration  would  extend  through  the  lower 
annular  ring  in  a  manner  similar  to  the  present  MITR 
design.   The  absorber  section  would  protrude  into  the  re- 
flector region  adjacent  to,  but  outside  of  the  inner  tank 
wall.   The  control  rod  drive  motors  and  connecting  shaft 
work  could  remain  and  be  used  in  their  current  positions. 
Maintenance  could  be  performed  on  the  control  rod  mechanisms 
by  removing  the  upper  shield  plug,  thus  allowing  a  man  to 
enter  the  area  between  the  upper  and  lower  annular  rings. 


37 

From  the  above  description  it  would  appear  as  if  the 
above  described  core  tank  design  with  a  double,  rotating 
shield  plug  mechanism  has  no  drawbacks  and  should  have  been 
adopted.   Where  the  above  design  fails  is  that  it  does  not 
sufficiently  conform  to  two  of  the  design  criteria  mentioned 
in  the  thesis  introduction  which  have  played  a  major  role  in 
the  selection  of  final  mechanical  designs  for  the  High  Flux 
MITR.   These  criteria  are  that  the  design  selected  should  be 
the  ultimate  in  simplicity,  and  should  have  no  hidden 
mechanisms.   The  mechanics  involved  in  the  concept  of  a 
double  rotating  shield  plug  with  an  inflatable  helium  seal 
are  not  what  could  be  described  as  simple,  and  the  mechanisms 
associated  with  the  lower  rotating  shield  plug  could  be  con- 
sidered as  hidden. 
2.3  Core  Support  Mechanism 

2.3.1  General 

Of  major  concern  to  the  redesign  of  the  MITR  core  area 
was  the  method  by  which  the  core  is  to  be  supported.   It  is 
possible  that  the  core  receive  its  primary  support  either 
from  below  or  by  being  hung  from  above.   The  core  support 
mechanism  selected  for  use  in  the  phase  I  High  Flux  MITR 
(designed  by  Mr.  E.  Barnett)  supports  the  core  from  below. 
The  overwhelming  reason  for  the  choice  of  this  mechanism  to 
support  the  core  from  below  was  safety.   In  the  design  of  a 
reactor  in  which  the  sudden  removal  of  a  fuel  element  from 
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its  core  position  would  cause  a  reactivity  addition  on  the 
order  of  3%   in  AK/K,  as  is  the  case  for  the  High  Flux  MITR, 
it  is  imperative  that  each  fuel  element  be  supported  in  an 
absolutely  foolproof  manner  (4). 

There  is,  however,  one  very  compelling  reason  for 
wanting  to  support  the  core  from  above.   The  aluminum  struc- 
ture necessary  to  support  the  core  from  below  results  in  a 
loss  on  the  order  of  20^  in  thermal  neutron  flux  in  the 
area  below  the  core  (5> ) .   It  is  this  area  just  below  the 
core  at  which  the  beam  ports  "look".   In  addition,  having 
the  aluminum  structure  below  the  core  results  in  an 

increased  gamma  flux  seen  by  the  beam  ports  due  to  the  ther- 

27 
mal  neutron  activation  of  Al   . 

2.3.2  Grid  Support 

The  first  design  proposed  as  a  means  of  supporting  the 
core  from  above  is  shown  in  Figure  2.3.2.1.   In  this  design 
the  core  is  hung  from  a  grid-like  structure,  the  grid  being 
supported  either  directly  by  the  lower  annular  ring  in  the 
case  of  the  double  shield  plug  design,  or  from  the  base  of 
the  necked  region  of  the  core  tank  in  the  case  of  the  necked 
core  tank  design.   An  adaptor  which  is  capable  of  being 
attached  to  either  end  of  the  fuel  element  would  be  connec- 
ted by  four  aluminum  rods  to  the  corners  of  a  rectangular 
plate.   The  plate  would  rest  on  top  of  the  grid  support  bars 
when  the  fuel  element  is  inserted  into  the  core.   A  hinged 
cover  would  be  lowered  over  the  grid  structure  and  latched 
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in  order  to  hold  the  fuel  element  support  plates  down. 
There  is  both  a  major  design  concept  advantage  and  a  major 
disadvantage  in  this  method  of  core  support.   The  advantage 
is  that  with  such  a  method  both  the  core  and  the  control 
elements  are  supported  from  the  same  structure,  the  lower 
annular  ring.   Therefore,  should  a  major  catastrophe  (such 
as  an  airplane  crashing  into  the  reactor  or  an  earthquake) 
occur  in  which  the  reactor  structure  was  disturbed,  having 
the  core  and  control  elements  supported  from  the  same 
structure  would  markedly  increase  the  probability  of  their 
staying  together;  i.e.,  one  should  not  fall  without  the 
other.   The  disadvantage  to  the  above  method  of  core  sup- 
port is  that  there  must  be  a  space  between  each  fuel  element 
of  width  equal  to  the  width  of  a  grid  support  bar.   The  grid 
support  bars  must  be  approximately  1/4  in.  wide  in  order  to 
provide  a  sufficient  area  on  which  to  rest  the  fuel  support 
plates.   This  spacing  would  increase  the  core  diameter  by 
approximately  1-1/2  in.,  thereby  decreasing  the  core  power 
density,  thereby  decreasing  the  thermal  flux  peaking  in  the 
reflector. 

2.3.3  Sliding  Bar  Support 

In  order  to  support  the  core  from  above,  but  without 
the  spacing  between  fuel  elements  required  by  the  method 
of  the  preceeding  section,  a  design  utilizing  sliding 
support  bars  was  developed  (see  Figure  2.3.3.1).    (A 
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design  of  this  general  type  was  proposed  earlier  by  Dr. 
Thompson  and  Mr.  Barnett.)  Assuming  for  the  sake  of 
illustration,  a  24  element  core  configuration,  a  system  of 
14  east-west  and  14  north-south  tracks  would  be  constructed 
around  the  core.   The  tracks  would  be  attached  at  one  end 
to  the  outer  tank  wall  and  supported  at  the  inner  end  by 
the  inner  tank  wall  plus  a  brace  hung  from  the  necked 
region  of  the  reactor  tank.   On  these  tracks  would  slide 
the  core  support  bars,  seven  bars  for  each  quadrant  of  the 
core.   Special  frame  type  structures  with  keyhole  indenta- 
tions would  lock  either  end  of  the  fuel  elements  and  accept 
the  sliding  core  support  bars  (figure  2.3.3.2).   In  order 
to  refuel  the  support  bars  would  be  withdrawn  in  such  a 
manner  as  to  allow  only  one  fuel  element  to  be  released  at 
a  time.   An  interlock  mechanism  would  prevent  the  with- 
drawal of  support  bars  in  any  manner  which  would  cause 
more  than  one  fuel  element  to  be  released.   The  fuel  ele- 
ment support  structures  being  frame-type  structures  with 
openings  on  all  sides,  the  DpO  coolant  would  flow  up  through 
the  core  and  out  the  sides  of  the  support  structure. 
Although  not  a  necessity,  the  top  of  the  support  structure 
could  be  made  solid  in  order  to  help  prevent  the  DgO  cool- 
ant  (with  its  N   activity)  from  being  circulated  into  the 
upper  portion  of  the  core  tank  after  leaving  the  core.   As 
shown  in  Figure  2.3.3.3,  the  sliding  support  bar  design 
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.,'ould  require  that  the  control  elements  be  angular  shaped 

to  conform  to  the  outline  of  the  core. 

The  major  shortcoming  of  the  sliding  bar  support 
mechanism  is  the  necessity  to  design  an  interlock  which 
would  be  absolutely  foolproof  in  preventing  the  accidental 
release  of  a  fuel  element.   It  can  also  be  argued  that 
this  design  certainly  does  not  meet  the  "simplicity"  design 
requirement. 

2.3.4  Safety  Supports 

From  a  safety  point  of  view,  both  the  designs  proposed 
in  Sections  2.3.2  and  2.3.3  would  be  enhanced  if  some  backup 
means  of  support  were  provided  below  each  fuel  element.   For 
this  purpose,  a  1/4  in.  aluminum  rod  set  into  the  coolant 
entrance  plenum  would  suffice.   On  top  of  the  1/4  in. 
aluminum  rod  would  be  a  structure  resembling  a  crow's  foot 
which  would  act  both  to  support  the  fuel  element  and  to 
turbulate  the  coolant  as  it  approaches  the  entrance  to  the 
fuel  element.   The  support  bar  could  be  removed  and  replaced 
by  a  bar  of  different  length  should  it  be  desired  to  change 
the  height  of  any  element  in  the  core.   The  aluminum  support 
bars  serve  only  one  function,  that  of  acting  as  a  backup  to 
ensure  that  no  fuel  element  is  able  to  fall  through  the 
core. 
2.4  Control  Rod  Drive  Mechanisms 

2.4.1  General 

One  of  the  most  difficult  design  problems  encountered 
in  the  redesign  of  the  MIT  Reactor  was  the  location  of  the 
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control  rod  drive  mechanisms.   Since  it  was  desired  to  re- 
tain the  present  concept  of  a  magnetic  coupling  between 
the  control  absorber  and  withdrawal  shaft,  it  became  a 
design  premis  that  the  control  rod  drive  mechanism  be  de- 
signed in  such  a  manner  as  to  allow  the  magnet  to  be  loca- 
ted as  far  from  the  core  as  possible,  and  if  possible  out- 
side the  DpO  tank.   By  doing  so,  it  was  hoped  that  magnet 
failures  due  to  radiation  damage  and  moisture  corrosion 
could  be  eliminated,  or  at  least  reduced  markedly  in  compari- 
son with  the  rate  of  failures  in  the  present  MITR.   The 
magnetic  coupling  between  the  control  rod  drive  motor  and 
the  control  absorber  was  envisioned  to  be  either  the  general 
type  used  in  the  present  MITR  or  a  magnetic  clutch  arrange- 
ment.  Some  of  the  proposed  designs  and  locations  for  the 
control  rod  drive  mechanisms  are  discussed  in  the  paragraphs 
which  follow: 

2.4.2  Absorber  and  Drive  Mechanisms  Outside  Inner  Core 
Tank 

The  first  control  rod  drive  mechanism  considered  was 
merely  a  modification  of  the  present  mechanisms  (see  Figures 
1.2.1  and  1.5.1).   In  this  design  (shown  in  Figure  2.2.1), 
the  absorber  would  be  located  in  the  reflector  region  just 
outside  the  inner  core  tank.   The  absorber  would  be  connec- 
ted by  a  shaft  to  the  electromagnet  coupling  in  such  a 
manner  that  the  magnet  is  always  located  inside  the  concrete 
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shielding  of  the  annular  ring  and  therefore  also  above  the 
level  of  the  reflector  DpO.   By  locating  the  magnets  away 
from  the  core  in  a  shielded  area  above  the  DpO  level,  the 
only  weakness  in  the  present  control  rod  design  has  been 
virtually  eliminated.   What  remains  is  a  control  rod 
mechanism  which  has  proven  very  reliable  over  the  past 
eight  years  of  reactor  operation.   The  control  rods  could 
be  removed  for  repair,  or  replacement  of  the  absorber  sec- 
tion by  unlocking  and  pulling  the  entire  control  rod  assem- 
bly up  through  the  annular  ring  and  out  through  the  top  of 
the  reactor  (after  first  removing  the  upper  shield  plug). 
Approximately  ten  control  rods  would  be  required  to  surround 
the  core.   For  a  description  of  the  present  control  rod 
assembly,  see  Section  1.5  of  this  thesis. 

2.4.3  Drive  Mechanism  Through  Annular  Ring 
A  modification  of  the  control  rod  mechanism  described 
in  Section  2.4.2  was  designed  in  order  that  the  control 
absorber  could  be  made  to  completely  surround  the  core  in 
the  area  immediately  adjacent  to  the  fuel  elements,  inside 
the  inner  core  tank  (see  Figure  2.4.3.1).   As  with  the 
design  of  the  preceeding  section,  the  control  rod  magnet  is 
located  up  within  the  shielding  area  of  the  annular  ring. 
However,  unlike  the  preceeding  design,  it  is  necessary  to 
provide  a  seal  in  order  to  prevent  DpO  from  leaking  up  into 
the  magnet  area.   This  seal  could  take  the  form  of  a  helium 
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pressure  in  and  above  the  magnet  area  equal  to  the  water 
head  below  the  magnet  area.   It  would  be  necessary  to 
circulate  the  helium  in  the  magnet  area  to  remove  the  D20 
vapor  which  would  otherwise  build  up  in  the  helium  atmos- 
phere  surrounding  the  magnet  area. 

The  control  rod  drive  mechanism  would  be  identical 
to  the  mechanism  used  in  the  present  MITR  as  far  down  as 
the  magnet.   The  assembly  below  the  magnet  would  consist 
of  an  "L"-shaped  arm  connecting  the  control  absorber  to 
the  magnet  armature.   The  absorber  section  would  be  pinned 
to  the  L-shaped  arm  so  that  in  order  to  replace  the  absor- 
ber section  it  would  be  merely  necessary  to  unpin  the  ab- 
sorber from  the  arm  and  lift  it  out  through  the  necked 
region  of  the  reactor  tank.   The  elbow  of  the  L-shaped  arm 
would  be  hinged  and  contain  a  dual  gear  arrangement  capa- 
ble of  raising  and  lowering  the  "forearm".   With  the  fore- 
arm in  the  fully  raised  position,  the  entire  control  rod 
drive  assembly  could  be  removed  out  through  the  reactor 
top.   It  can  be  noted  at  this  point  that  should  either  the 
elbow  of  the  L-shaped  arm  or  the  pin  connecting  the  arm  to 
the  absorber  section  fail,  the  reactor  would  be  shut  down. 
As  shown  in  Figure  2.4.3.1,  a  spring-type  dashpot  would 
slow  the  absorber  support  arm  as  it  neared  the  base  of  its 
travel,  and  aluminum  shields  would  guide  the  control  ab- 
sorbers as  they  descend  into  the  core.   Eight  arced  sections 
of  absorber  would  be  used  to  surround  the  core. 
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2.4.4  Drive  Mechanisms  Above  DpO  Level 


In  order  to  eliminate  the  complexity  involved  with 
having  the  control  rod  drive  mechanisms  penetrate  through 
the  annular  ring  of  shielding,  designs  were  proposed  which 
would  locate  the  magnets  and  drive  mechanisms  on  the  reac- 
tor top,  entirely  above  the  DpO  level  in  the  tank.   Three 
possible  design  arrangements  were  proposed  by  this  author. 

The  first  design  proposed  called  for  horizontally 
placed  withdrawal  mechanisms  installed  flush  with  the 
reactor  top  (Figure  2.4.4.1).   The  control  absorber  would 
be  attached  to  the  magnet  armature  by  a  cable  extending 
vertically  to  the  surface  of  the  DpO  column,  thence  hori- 
zontally to  the  position  of  the  armature.   Eight  curved 
absorber  plates  would  form  a  cylindrical  shell  completely 
surrounding  the  core.   The  entire  withdrawal  mechanism 
assembly  would  be  easily  accessible  through  removable 
covers  on  the  reactor  top.   The  cable  withdrawal  mechanism 
provides  excellent  protection  for  the  magnet,  but  suffers 
one  serious  disadvantage.   The  absorber  will  only  return 
to  the  core  by  the  force  of  gravity.   It  cannot  be  driven 
into  the  core,  thus  ensuring  that  the  absorber  is,  in  fact, 
completely  inserted  to  its  furthest  point. 

In  order  to  provide  for  a  means  whereby  the  control 
absorber  section  could  be  driven  into  the  core,  a  design 
was  proposed  whereby  the  withdrawal  mechanisms  would  be 
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located  directly  above  the  absorber  sections.   The  with- 
drawal mechanisms  would  be  situated  vertically  atop  the 
reactor  top  shielding  and  connected  to  the  absorber  sec- 
tions by  solid  aluminum  rods  (Figure  2.4.4.2).   As  with 
the  previous  design,  eight  arced  absorber  sections  would 
completely  encircle  the  core.  The  shaft  connecting  the 
magnet  armature  with  the  control  absorber  would  be  compo- 
sed of  two  sections  joined  together  by  a  screw-in  fitting. 
The  two  sections  could  be  disconnected  at  a  position 
parallel  to  the  reactor  top,  thus  allowing  a  cover  shield 
to  be  slid  into  place  over  the  D^O  column  during  fuel 
flipping  and  fuel  changing  evolutions.  A  key  at  the  base 
of  the  withdrawal  mechanism  housing  and  a  gear  just  below 
the  armature  would  mate  to  allow  the  upper  rod  section  to 
be  unscrewed.   The  penetrations  through  the  reactor  top 
shielding  for  the  aluminum  rods  would  have  to  provide 
adequate  clearance  so  as  to  guarantee  that  the  aluminum 
rod  would  always  be  free  to  travel  through  the  shielding. 
A  pressure  equalizing  line  would  be  required  in  order  to 
euqalize  the  pressure  of  the  area  below  the  reactor  top 
shielding  with  that  in  the  area  above  the  magnet.   The 
equalizing  line  would  be  required  to  have  a  cross  sectional 
area  greater  than  the  cross  sectional  area  of  the  clearance 
surrounding  the  aluminum  rod  as  the  rod  passes  through  the 
reactor  top  shielding.   The  purpose  of  the  pressure 
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equalizing  line  would  be  to  ensure  that  a  sudden  increase 
in  reactor  pressure  would  not  exert  itself  on  the  bottom 
of  the  magnet  armature,  thus  using  the  armature  as  a  pis- 
ton to  drive  the  control  rods  upward.   With  the  entire 
control  rod  withdrawal  mechanism  located  on  top  of  the 
shielded  refueling  enclosure,  all  components  would  be 
readily  available  for  maintenance. 

A  third  design  alternative  proposed  would  employ 
features  of  both  the  above  designs.  Vertical  aluminum 
shafts  would  lead  from  the  eight  arced  control  absorber 
sections  to  just  above  the  DpO  level  in  the  necked  tank. 
A  horizontal  shaft  would  be  coupled  to  the  vertical  shaft 
by  matching  gears.   The  horizontal  shaft  would  in  turn 
lead  to  a  magnetic  clutch  and  thence  to  the  control  rod 
drive  motor  (Figure  2.4.4.3).   Upon  interruption  of  elec- 
trical current,  the  magnetic  clutch  would  release  the 
horizontal  shaft,  thus  allowing  it  to  turn  freely  and 
allowing  the  absorber  to  drop  down  into  the  core.   With 
this  arrangement,  the  absorber  could  also  be  driven  into 
the  core  should  the  need  arise  for  such  a  feature.   It  is 
a  requirement  that  all  of  the  gears  connecting  the  horizon- 
tal and  vertical  control  rod  shafts  turn  freely  at  all 
times  in  order  to  be  effective  in  a  scram  condition.   Such 
a  requirement  could  be  presented  as  a  design  flaw. 
However,  it  was  felt  that  the  probability  of  two  or  more 
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gear  couplings  jamming  simultaneously  was  so  small  as  to 
be  nearly  nonexistent.   Any  one  of  the  eight  control  absor- 
bers would  be  sufficient  to  shut  down  the  reactor. 

2.4.5  Drive  Mechanisms  for  Inner  Control  Elements 
In  reviewing  the  safety  requirements  for  the  High 
Flux  MITR,  it  was  decided  by  Dr.  Thompson  that  sufficient 
control  must  exist  to  maintain  the  core  in  a  subcritical 
condition,  even  if  the  core  were  to  be  filled  with  light 
water.   In  order  to  satisfy  this  requirement,  Dr.  Thompson 
suggested  that  the  core  should  contain  an  inner  ring  of 
control  elements.   The  necessity  for  an  incore  ring  of 
control  elements  posed  a  considerable  design  problem.   A 
mechanism  which  would  not  interfere  with  the  removal  or 
replacement  of  any  fuel  element,  yet  would  provide  for 
reliable  control  of  the  inner  absorber  elements  was  re- 
quired.  The  design  proposed  by  this  author  is  shown  in 
Figure  2.4.5.1.   The  inner  control  rod  withdrawal  mechanism 
would  consist  of  two  concentric  shafts.   The  inner  shaft 
would  be  attached  at  its  base  to  four  pivoted  support  bars 
capable  of  swinging  up  or  down.   The  support  bars  would  in 
turn  be  connected  to  the  outer  shaft  by  four  aluminum  rods, 
each  pivoted  at  either  end.   The  connection  to  the  outer 
shaft  would  be  accomplished  by  a  clamp-type  collar  which 
would  ride  on  threaded  grooves  cut  onto  a  portion  of  the 
outer  shaft.   By  rotating  the  outer  shaft,  the  collar  could 
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be  caused  to  move  upward  or  downward,  thus  raising  or 
lowering  the  support  bars.   In  their  lowered  (vertical) 
position  the  support  bars  would  occupy  an  area  above  the 
center  of  the  core  no  larger  than  the  area  of  the  cryostat 
planned  for  the  center  of  the  core.   Consequently,  the 
mechanism  would  not  hinder  the  removal  of  any  fuel  element, 
In  their  raised  (horizontal)  position,  the  support  bars 
would  engage  slots  in  the  upper  section  of  the  inner  con- 
trol elements.   The  entire  mechanism  would  be  raised  in 
order  to  withdraw  the  inner  control  elements  from  the  core, 
(It  should  be  noted  that  the  inner  absorber  elements  would 
require  followers  in  order  to  prevent  the  flux  peaking 
which  would  occur  if  DpO  were  allowed  to  fill  the  space 
vacated  by  the  absorber  as  the  control  elements  are  with- 
drawn. )   The  withdrawal  mechanism  would  utilize  a  magnetic 
coupling  and  be  constructed  identical  to  the  withdrawal 
mechanisms  for  each  of  the  absorber  sections  for  the  outer 
control  elements. 
2.5  Refueling  Mechanism 

2.5. 1  General 

The  initial  design  concept  for  refueling  and  fuel 
flipping  was  to  employ  a  refueling  flask  with  a  self- 
contained  flipping  mechanism.   However,  after  the  concept 
of  a  permanent,  enclosed  refueling  and  fuel  flipping 
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facility  to  be  located  directly  above  the  core  tank  on 

the  reactor  top  was  proposed  by  Mr.  Barnett,  a  new  area 

for  design  was  created.   This  author  set  about  the  task 

of  developing  a  refueling  and  fuel  flipping  scheme  which 

could  be  used  in  such  an  enclosure.   In  the  process,  two 

quite  similar  designs  were  developed,  each  of  which  is 

described  below. 

2.5.2  Fuel  Flipper  Above  Core  Tank 

In  order  to  provide  a  means  of  rotating  a  fuel  element 
end  for  end  a  rotating  table  was  designed  to  be  mounted 
vertically  against  a  side  wall  of  the  refueling  enclosure. 
On  the  rotating  table  would  be  a  bracket  designed  to  accept 
a  fuel  element.   Directly  below  the  rotating  table  would  be 
another  bracket  designed  to  accept  the  lower  adapter  unit 
of  the  fuel  element.   To  one  side  of  the  DpO  column  would 
be  a  horizontal  shute  through  which  fresh  fuel  would  enter, 
and  a  downward  sloping  shute  through  which  spent  fuel  would 
be  discharged  into  a  small  D?0  filled  spent  fuel  pool.   A 
leaded  glass  window  would  be  provided  in  the  side  wall 
opposite  the  rotating  table  in  order  to  provide  an  operator 
with  visual  access  to  the  refueling  cavity  during  refueling 
or  flipping  operations.   A  refueling/fuel  flipping  sequence 
would  proceed  somewhat  as  follows  (assuming  control  rod 
drives  to  be  located  atop  the  refueling  enclosure) :   A  fuel 
grabbing  tool  would  be  inserted  into  the  reactor  tank  from 
the  reactor  top,  latch  onto  the  desired  fuel  element,  and 
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FIGURE     2.5.2.1 
FUEL    FLIPPER    MECHANISM 
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pull  it  up  into  the  refueling  enclosure.   Meanwhile  the 
control  rod  drive  shafts  would  be  disconnected  at  the 
surface  of  the  DpO  and  the  upper  portion  of  the  drive 
shaft  withdrawn  into  the  withdrawal  mechanism  housing.   A 
2  in.  thick  lead  shield  would  next  he  slid  into  position 
above  the  DoO  column.   The  lead  shield  would  serve  to 
prevent  the  fuel  element  from  falling  back  into  the  reactor 
core  area  should  it  accidentally  become  unlatched.   (The 
lead  shield  would  also  assist  in  isolating  the  refueling 
enclosure  from  the  core  while  performing  maintenance  on 
any  of  the  refueling  or  control  rod  drive  mechanisms.) 
With  the  lead  cover  in  place,  a  clamping  tool  would  be  in- 
serted through  the  side  of  the  refueling  enclosure  and 
clamp  the  fuel  element  to  the  rotating  table.   Once  secure, 
the  tool  used  to  raise  the  fuel  element  out  of  the  core 
would  be  disconnected  and  removed.   Another  tool  would  be 
inserted  into  the  refueling  enclosure  to  pull  out  the  pins 
connecting  the  fuel  element  to  its  adapter.   Once  the  pins 
are  pulled,  the  adapter  would  be  allowed  to  drop  down  two 
inches,  thus  giving  the  fuel  element  sufficient  clearance 
to  rotate  away.   If  the  fuel  element  were  to  be  replaced, 
the  table  would  rotate  to  align  the  fuel  element  with  the 
chute  to  the  spent  fuel  pool,  the  element  would  be  dis- 
charged, and  table  rotated  to  receive  a  new  fuel  element 
from  the  horizontal  chute.   If  the  fuel  element  were  to  be 
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flipped,  the  table  would  simply  be  rotated  l8o°.   The 
remainder  of  the  refueling/flipping  operation  would  be 
the  reverse  of  the  above.   The  lower  adapter  would  be 
raised  and  attached  to  the  fuel  element,  the  fuel  retrac- 
tion/loading tool  attached  to  the  fuel  element,  the  fuel 
element  and  fuel  element  adapter  clamps  removed,  the  lead 
shield  cover  removed  from  the  tank  top,  the  control  rod 
drive  shafts  reconnected,  and  the  fuel  element  lowered 
into  its  respective  position  in  the  core. 

After  allowing  a  spent  fuel  element  sufficient  time 
to  cool  (radiation  wise),  the  spent  fuel  pool  would  be 
drained  and  the  fuel  element  removed  through  a  plug  in  the 
refueling  enclosure  shielding,  and  thence  taken  to  the 
main  spent  fuel  pool  in  the  reactor  building  basement. 
2.5.3  Fuel  Flipper  Away  from  Core  Tank  Top 
In  order  to  remove  the  fuel  flipping  operation  from 
the  area  directly  above  the  core,  a  modification  to  the 
design  of  the  preceeding  paragraph  was  proposed.   This 
design  is  shown  in  Figure  2.5.3.1.   The  purpose  in  remov- 
ing the  mechanisms  involved  in  fuel  flipping/refueling  from 
the  area  immediately  above  the  core  to  an  area  which  could 
be  isolated  from  the  space  directly  above  the  core  was  to 
produce  a  design  whereby  maintenance  could  be  performed  on 
the  fuel  flipping/refueling  mechanisms  while  the  reactor 
was  operating  at  full  power.   In  this  manner  it  was 


FIGURE       2.5.3.1 

FUEL   FLIPPER     MECHANISM- 
ALTERNATE   APPROACH 


63 


64 

anticipated  that  reactor  down  time  necessitated  by  mainten- 
ance and  repairs  could  be  minimized. 

The  rotating  table  would  be  separated  from  the  area 
above  the  core  tank  by  a  o  in.  lead  shielding  wall.  A 
removable  lead  shutter  in  the  center  of  the  shielding  wall 
would  raise  to  allow  a  fuel  element  and  its  lower  adapter 
to  pass  from  the  core  tank  area  to  the  rotating  table.   The 
fuel  element  clamping  tool  would  extend  through  the  center 
of  the  shaft  of  the  rotating  table.   The  clamping  tool  for 
the  fuel  element  adapter  would  penetrate  the  shielding  wall 
just  below  the  rotating  table.   Instead  of  being  discharged 
into  a  small  spent  fuel  pool  on  the  reactor  top,  spent  fuel 
would  be  removed  through  a  refueling  door  into  a  small  fuel 
element  transport  flask  and  transported  directly  to  the 
large  spent  fuel  pool  in  the  reactor  building  basement. 
Fresh  fuel  elements  would  also  enter  the  refueling  enclosure 
through  the  refueling  door.   Two  leaded  glass  windows  would 
be  provided  to  enable  visual  observation  and  control  of  the 
entire  refueling/flipping  process.   One  window  would  be  in 
the  wall  adjacent  to  the  DpO  column,  facing  the  rotating 
table.   The  other  would  be  in  the  enclosure  containing  the 
rotating  table,  facing  the  refueling  door. 
2.6  Underwater  Periscope 

While  visual  access  to  the  refueling  enclosure  would 
be  provided  by  one  or  more  leaded  glass  windows  placed  in 
the  walls  of  the  refueling  enclosure,  it  is  also  desirable 
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that  a  means  be  provided  for  looking  down  into  the  core 
itself.   Being  able  to  see  the  core  would  be  a  necessity 
during  refueling  operations,  and  in  addition  would  be 
valuable  for  instructional  purposes.   It  would  be  possible 
to  locate  a  leaded  glass  window  in  the  upper  corner  of  the 
refueling  enclosure,  angled  downward  in  such  a  manner  as 
to  allow  one  to  observe  the  core.   However,  with  such  an 
arrangement  one  would  not  get  a  clear  view  of  the  core  as 
a  result  of  attempting  to  look  through  a  gas -liquid  inter- 
face at  the  core  tank  surface.   During  reactor  operation, 
the  liquid  (D20)  surface  would  be  turbulent  due  to  the 
rapid  circulation  of  the  coolant.   The  turbulence  alone 
would  hamper  vision,  but  in  addition  it  is  necessary  to 
shine  a  light  down  into  the  coolant  tank  in  order  to  see 
the  core.   Some  of  this  light  would  be  reflected  back  from 
the  coolant  surface,  thus  creating  what  would  be  the  major 
hindrance  to  clear  observation  of  the  core. 

In  order  to  eliminate  the  difficulties  mentioned  in 
the  previous  paragraph  associated  with  observing  the  core 
by/  way  of  the  refueling  enclosure,  a  retractable  periscope 
was  designed.   As  shown  in  Figure  2.5.2.1,  the  periscope 
is  mounted  at  an  angle  in  such  a  way  that  when  fully  inser- 
ted into  the  core  the  outside  end  of  the  periscope  remains 
above  the  level  of  the  coolant.   In  this  manner  the  connec- 
tion between  the  periscope  and  the  vertical  viewing  tube  is 
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not  required  to  span  a  liquid  interface.   When  retracted, 
the  tip  of  the  periscope  would  be  flush  with  the  edge  of 
the  coolant  tank  wall,  thereby  eliminating  any  interfer- 
ence to  the  insertion  or  removal  of  components  in  the 
area  of  the  periscope.   Although  the  entire  core  would  be 
visible  when  looking  through  the  periscope,  the  mirror 
which  extends  into  the  tank  would  be  able  to  pivot  in 
order  that  one  might  scan  for  some  distance  up  the  sides 
of  the  core  tank  wall.   The  periscope  would  have  its  own 
self-contained  light  source.   The  light  source  would  be 
located  at  the  base  of  the  periscope,  well  removed  from 
the  intense  radiation  area,  and  project  its  light  onto  the 
core  via  a  mirror  near  the  periscope  tip.   One  would  look 
into  the  periscope  through  a  vertical  viewing  tube  which 
would  extend  to  a  height  of  4  ft.  above  the  periscope. 
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Chapter  3 
SHIELDING  CALCULATIONS 
3. 1  General 

The  reactor  top  shielding  design  and  supporting  calcu- 
lations were  originally  done  for  the  phase  one  reactor 
design.   This  chapter  will  first  present,  therefore,  the 
shielding  design  and  calculational  work  done  for  a  D20 
moderated,  cooled,  and  reflected  core  in  some  detail. 
When  the  core  moderator  and  coolant  were  changed  to  light 
water,  a  new  set  of  shielding  calculations  was  required. 
Since  the  method  and  reasoning  employed  in  the  phase  two 
calculations  was  exactly  the  same  as  that  employed  for  the 
phase  one  calculations,  phase  two  calculations  will  not  he 
reported  in  detail,  hut  rather  a  summary  of  required 
shielding  and  calculated  dose  rates  will  be  given  in  Sec- 
tions 3.7-3.9. 

In  phase  one,  based  on  the  assumption  that  the  reactor 
top  shielding  design  would  consist  of  an  enclosure  similar 
to  that  shown  in  Figure  3.1.1,  it  was  necessary  to  deter- 
mine both  the  composition  and  thickness  of  the  shielding 
which  would  be  required.   As  the  calculations  of  this 
chapter  will  show,  no  major  shielding  problems  were  encoun- 
tered.  Such  was  not  known  in  the  early  stages  of  the 
design  project,  however.   It  was  anticipated  that  shielding 
requirements  during  reactor  operation  might  be  dictated, 


FIGURE     3.1.1 
REACTOR  TOP    SHIELDING  ENCLOSURE 
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not  by  the  radiation  being  emitted  from  the  core  but  by 

the  radiation  dose  due  to  gamma  emission  from  N   decay 

3 
or  the  gamma  emission  accompanying  the  formation  of  H  in 

the  necked  region  of  the  core  tank,  to  mention  but  one 
example.   While  the  calculations  which  follow  serve  to 
justify  the  chosen  composition  and  thickness  of  the  reac- 
tor top  shielding  (shown  in  Figure  3.1.1),  it  should  be 
noted  that  the  shield  composition  and  thickness  were 
originally  determined  by  successive  iterations  of  exactly 
similar  calculations.   It  is  for  the  sake  of  brevity  and 
clarity  that  only  the  final  set  of  calculations  is  presented 
here.   No  shielding  calculations  were  performed  in  the 
radial  direction  since  neither  reactor  power  nor  radial 
shielding  will  undergo  any  change  from  the  present  MITR. 
In  fact,  radial  attenuation  will  be  improved  since  the  core 
boundary  will  be  further  from  the  shielding  walls. 

It  was  required  that  the  DpO  column  combined  with  the 
shielding  of  the  refueling  enclosure  be  so  designed  as  to 
be  capable  of  reducing  the  dose  on  the  reactor  top  suffici- 
ently such  that  a  worker  working  full  time  in  this  area 
would  not  receive  a  dose  in  excess  of  the  maximum  whole 
body  dose  limit  established  by  the  U.S.  Atomic  Energy 
Commission.   The  Atomic  Energy  Commission's  maximum  allow- 
able whole  body  dose  levels  are  100  mrem  per  ^-0  hour  week 
or  2.5  mrem  per  hour  (6).   However,  it  is  desirable  that 
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the  shielding  "be  capable  of  reducing  the  dose  anywhere  on 
the  reactor  top  to  well  below  these  maximums.   The  present 
MITR  was  designed  so  that  the  dose  on  the  reactor  top 
would  be  less  than  1  mrem  per  hour  at  5  megawatts  (8). 
Similarly,  the  goal  in  the  design  of  the  shielding  for  the 
High  Flux  MITR  was  to  reduce  the  dose  anywhere  on  the  reac- 
tor top  to  less  than  1.0  mrem  per  hour  when  the  reactor  is 
at  a  power  of  5  Mw. 

The  dose  rates  during  both  shutdown  and  operational 
conditions  were  determined.   For  the  shutdown  condition, 
the  sources  of  radiation  of  most  concern  are  those  resulting 
from  material  activation  and  fission  product  decay.   For  the 
operational  condition  one  must  consider,  in  addition  to  the 
shutdown  sources,  the  dose  due  to  gamma  ray  production  in 
fission,  fission  neutrons,  photoneutrons,  gamma  rays  from 

-i  r 

the  decay  of  N  ,  and  gamma  rays  from  the  production  of 

tritium. 

3.2   Calculational  Methods 

3.2.1  General 

The  calculational  methods  followed  in  performing  the 
shielding  calculations  for  the  High  Flux  MITR  were  the 
same  as  those  used  in  the  design  of  the  shield  for  the 
Omega  West  Reactor  at  Los  Alamos  Scientific  Laboratory 
(see  reference  7).   This  same  procedure  was  also  used  in 
obtaining  the  initial  shielding  estimates  for  the  present 
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MITR.   Radiation  measurements  recorded  subsequent  to  the 
completion  of  the  Omega  West  Reactor  tend  to  bear  out  the 
validity  of  using  the  calculational  methods  employed  in 
the  shield  design.   As  an  example,  the  calculated  dose 
rate  due  to  gamma  rays  through  the  lid  on  the  reactor  top 
was  ll4  mrem/hr.   The  measured  value  was  107  mrem/hr. 

3.2.2  Gamma  Ray  Attenuation 

The  sources  of  gamma  rays  in  the  High  Flux  MITR  will 
be  treated  as  point  sources.   The  assumption  of  a  point 
source  has  been  estimated  by  Gage  to  produce  results  which 
are  at  most  20$  greater  than  those  obtained  from  calcula- 
tions performed  utilizing  a  distributed  source.   The  error 
involved  in  assuming  a  point  source  of  radiation  is  compen- 
sated for  an  adjustment  to  the  radiation  source  strength, 
as  discussed  in  a  later  section  of  this  chapter.   The  equa- 
tion for  the  attenuation  of  gamma  rays  from  a  point  source 
is: 


0 


— ^y  FBe"^    7's/cm2-sec    ,  (3.2.2.1) 

4ttR 


where,   S   =   source  strength   (7's/sec), 

R  =  distance   from  the  source    (cm), 

F  =   fraction  of  gamma  rays    in  energy  range  of 

interest, 
\±  =   linear  absorption  coefficient    (cm"   ), 
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x  =  thickness  of  absorbing  material  (cm),  and 

B  =  buildup  factor 

The  Omega  West  Reactor  shielding  calculations  were 

done  in  terms  of  decade  lengths  rather  than  in  terms  of 

the  linear  absorption  coefficient,  and  so  decade  lengths 

will  be  used  here  also.   A  decade  length  is  defined  as 

that  thickness  of  material,  d,    required  to  reduce  the 

flux  by  a  factor  of  ten.   Decade  length  is  determined  as 

follows : 

1/10  =  e^d 

d  =  2. 3/|a  cm  =  0.075VM.  ft. 

The  shield  thickness  in  decade  lengths,  D,  is  defined  as 
D  =  -t  ,  where  x  is  again  defined  as  the  thickness  of  absorb- 
ing material.   It  is  also  necessary  to  convert  the  buildup 
factor,  B,  from  its  normally  listed  value  in  terms  of  e- 
folding  lengths  to  decade  lengths.   Buildup  factors  in 
terms  of  decade  lengths  for  water,  iron,  and  lead  were 
taken  from  reference  (7)  and  are  given  in  Figures  3.2.2.1, 
3.2.2.2,  and  3.2.2.3  respectively.   It  should  be  noted 
that  Figure  3.2.2.1  is  for  light  water.   However,  since 
gamma  ray  attenuation  is  most  strongly  dependent  on  the 
atomic  number  of  the  absorber,  and  Hp0  and  DpO  both  have 
the  same  total  Z_  value,  Figure  3.2.2.1  is  assumed  to  be 
applicable  to  heavy  water  as  well.   This  is  a  conservative 
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FIGURE   3.2.2 
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FIGURE  3.2.2.2 
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FIGURE  3.2.2.3 
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assumption  since  D^O  is  about  10$  more  dense  than  Hp0. 
It  should  also  be  noted  that  buildup  factors  for  heavy 
concrete  were  not  available,  so  the  buildup  factors  for 
iron  (Figure  3.2.2.2)  were  used  in  the  calculations 
later  in  this  chapter.   This  substitution  is  reasonable 
since  heavy  concrete  is  about  75$  iron.   The  result  is 
conservative  from  a  shielding  point  of  view  because  the 
buildup  factor  for  iron  is  greater  than  that  of  concrete. 
Equation  3.2.2.1  can  now  be  re-written  as 

0  =  — £-*■  FB   10"D  7»s/cm2-sec    .  (3.2.2.2) 

4ttR 

In  applying  this  flux  attenuation  equation  to  a  shield 

consisting  of  many  layers,  it  is  assumed  that  the  buildup 

factor  for  the  shield  is  the  sum,  and  not  the  product, 

of  the  buildup  factors  for  each  layer.   Although  this  may 

not  be  a  conservative  assumption,  it  has  been  shov/n  by 

Gage  to  be  much  more  accurate  than  if  the  product  of  the 

buildup  factors  were  used. 

3.2.3  Neutron  Attenuation 

As  with  gamma  rays,  neutrons  will  be  considered  as 
being  emitted  from  a  point  source.   The  equation  for  the 
attenuation  of  fission  spectrum  neutrons  is : 

0  =  7   e 

47TIT 

where  N  =  neutron  source  strength  (neutrons /sec) , 
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R  =  distance  from  the  source  (cm),  and 
X  =  relaxation  length  for  removal  of  neutrons  (cm). 
A  decade  length  for  the  attenuation  of  fast  neutrons  is 
determined  as  follows : 

J_         "d^r 
10  ~  e 

d  =  2.3  Xr(cm)  =  0.0754  Ar  (ft). 

Relaxation  lengths  and  decade  lengths  for  some  shielding 
materials  of  interest  were  taken  from  reference  (7)  and 
are  given  in  Table  3.2.3.1. 

The  attenuation  of  thermal  neutrons  is  determined  by 
means  of  diffusion  theory.   The  applicable  diffusion  equa- 
tion for  thermal  neutrons  from  an  infinite  plane  source  in 
an  infinite  homogeneous  media  is : 

Le-x/L  2 

0=0     — go —  neutrons/cm  -sec  ,  (3.2.3.2) 

where  L  =  diffusion  length  (cm), 

D  =  diffusion  coefficient  (cm), 

x  =  thickness  of  attenuating  material  (cm),  and 

0^  =   thermal  neutron  flux  on  the  core  side  of  the 
o 

shield 
3.3  Shielding  Calculations,  Normal  Operation  (Phase  I) 

3.3.1  Fission  Gamma  Rays 

The  prim;,  ry  source  of  radiation  to  be  shielded  against 
is  the  fyimma  rays  produced  by  fissions  in  the  core. 
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Assuming  198.22  MEV  are  released  per  fission  (see  Section 
4.2),  the  source  strength  of  the  High  Flux  MITR  operating 
at  5  megawatts  would  be: 

(5  x  106  watts) (3.1  x  1010  fissions /watt-sec)  =  1.55  x  1017 

fissions/sec. 
With  each  fission  event  there  are,  on  the  average,  7.5 
gamma  rays  emitted  from  the  fissioning  nucleus  (9).   Thus, 
the  source  strength  of  the  fission  gamma  rays  from  a  core 
operating  at  5  megawatts  is : 

M  rr  v  in1?  fissions w>7  ^   gair  as  v   .  17   1f)l8  gammas 
[1.55   *  io   — — ;^.t)  fission^  ~  1'1(   x  1U   ""sec    ' 

Figure  3.3.1.1,  taken  from  reference  (7),    shows  the 
relative  number  of  fission  gamma  rays  versus  energy.   The 
straight  line  extrapolation  shown  in  the  figure  was  used 
to  calculate  the  percentage  of  gamma  rays  above  a  given 
energy.   The  results  obtained  are  conservative.   Gage  has 
shown  that  for  energies  above  1.5  MeV  the  extrapolated 
line  predicts  about  30  percent  more  gamma  rays  than  exist. 
The  excess  in  the  fraction  of  gamma  rays  predicted  above 
an  energy  E  in  the  range  0  to  1.5  MeV  decreases  from  30$ 
to  zero  as  E  decreases  from  1.5  to  0  MeV.   It  is  assumed 
that  the  fission  gamma  rays  can  be  lumped  into  three 
energy  groups.   Table  3.3.1.1  gives  the  calculated  frac- 
tions above  each  of  the  three  energy  groups  assumed  using 


80 


10 


FIGURE  3.3.1.1 
RELATIVE  NUMBER  OF  GAMMA  RAYS 
(FISSION)  vs.  ENERGY  ! 
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the  straight  line  analysis,  and  also  the  assumed  frac- 
tion at  each  energy  that  will  be  used  in  the  shielding 
calculations. 

Table  3.3.1.1 


CalcuL 

ated  Percentage 

Assumed  Percentage 

Energy 

Above  E 

at  E 

0  MeV 

100$ 

1  MeV 

33.4$ 

100$ 

3  MeV 

3.7$ 

33$ 

6  MeV 

0.14$ 

3$ 

Note  that  the  assumed  percentage  of  gamma  rays  at  E  „, 

r  °  calc 

have  been  greatly  overestimated.   For  example,  the 
assumed  percentage  of  gamma  rays  at  6  MeV  has  been  over- 
estimated by  about  a  factor  of  20.   This  overestimation 
of  gamma  source  strength  was  the  method  used  by  Gage  to 
compensate  for  the  nonconservative  point  source  assump- 
tion.  Gage  performed  a  twelve  energy  group  calculation 
with  energy  groups  from  0-12  MeV  in  order  to  check  the 
accuracy  of  the  three  energy  group  model.   For  the  water 
and  concrete  radial  shielding  of  the  Omega  West  Reactor, 
he  found  that  the  three  group  model  as  described  above 
was  conservative  by  a  factor  of  5.   Similar  conservatism 
could  be  expected  for  the  D20  and  lead  shielding  of  the 
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High  Flux  MITR. 

The  height  of  the  coolant  above  the  core  top  is  8  ft. 
Since  the  gamma  source  is  being  treated  as  a  point  source 
of  radiation,  it  is  necessary  to  add  to  this  the  additional 
1  ft.  which  exists  between  the  center  of  the  core  and  the 
core  top.   Hence,  the  total  thickness  of  DpO  to  be  consider- 
ed is  9  ft.   The  helium-filled  space  between  the  DpO  surface 
and  the  bottom  of  the  shielding  enclosure  is  4-1/2  ft.  in 
height.   The  refueling  enclosure  top  shielding  block  con- 
sists of  6  in.  of  masonite  sandwiched  between  two  1/4  in. 
sheets  of  boral  (neglected  for  gamma  ray  calculations),  1 
in.  of  steel  used  for  support,  and  1  ft.  of  lead.   (Note: 
DpO  cooling  coils  will  loop  through  the  lead  and  masonite- 
boral  sandwich.) 

For  clarity  Eq.  3.2.2.2  may  be  rewritten  as: 

qF  -(Dn+DT+DT) 

0  =   ^I -p(Bn+BT+BT)10   D  I   L   ,      (3.3.1.1) 

MVV*us>  D  X    L 

where  the  subscripts  D,  H,  I,  and  US  refer  to  DpO,  Helium, 
Iron,  Lead,  and  Upper  Shield  respectively.   For  calcula- 
tional  purposes,  the  helium  gap  will  be  treated  as  void. 
Decades  of  attenuation  and  buildup  factors  for  materials 
of  interest  as  a  function  of  energy  are  given  in  Table 
3.3.1.1. 
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Substituting  into  Eq.  3.3.1.1  for  the  dose  rate  due 
to  6  MeV  core  gamma  rays  on  the  reactor  top  directly 
above  the  core  yields : 

0,   =   (1-I7xl0l8)(0.03) (4.3+l.24l4.6)lo-(3-59+.26+6.7i) 

°       Hit   (9+5+.083+1)x30.48 

p 

0V  =  7.4-1  gamma  rays /cm  -sec. 

For  6  MeV  gammas  157.7  gam^a  rayS  =  1  mrem/hr.,  the  conver- 

„cm  -sec 
sion  from  gamma  rays/cm  -sec  to  mrem/hr  having  been  obtained 

from  page  7-66  of  reference  (9).   The  R.B.E.  for  gamma  rays 

was  taken  to  be  1.   Therefore,  0r   =  0.047  mrem/hr. 

Following  exactly  the  same  procedure  for  the  3  MeV  group, 

0n   =  5.4-2  gammas/cm  -sec  =  0.021  mrem/hr, 

and  for  the  1  MeV  group, 

O  p  1  -l 

0-.   -   4.4  x  10   gammas/cm  -sec  =  7.6  x  10"   mrem/hr. 
Summing  the  above  calculated  dose  rates,  the  total  dose  rate 
on  the  reactor  top  directly  above  the  core  due  to  fission 
gamma  rays  is  calculated  to  be: 

0   =  0.047  x  0.021  +  7.6  x  10*"11  =  0.068  mrem/hr. 

3.3.2  Material  Activation 
3.3.2.1  General 

The  primary  reaction  which  must  be  considered  is  the 
thermal  neutron  activation  of  the  aluminum  contained  in  the 
core  and  surrounding  structure.   The  reaction  of  concern  is: 
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Al27  +  n  >   Al28  +  7  (3.3.2.1.1) 

The  thermal  cross  section  for  this  reaction  is  0.215  barns. 

Al   emits  an  average  of  two  prompt  gamma  rays  for  each 

27 
neutron  captured  by  Al  '  (10) . 

Also  of  concern  from  a  shielding  point  of  view  are 

the  gamma  rays  released  in  the  decay  of  Nitrogen  16,  which 

is  formed  in  the  reaction: 

016  +  n  >     N16  +  p  (3.3.2.1.2) 

Although  the  activation  cross  section  averaged  over  the 
entire  fission  spectrum  for  this  section  is  only  20  micro 
barns,  it  is  usually  the  most  important  reaction  associated 
with  the  reactor  coolant.   Nitrogen  16  decays  with  a  half- 
life  of  7.35  sec,  liberating  gamma  rays  of  energy  6.1 
MeV  (75$)  or   7.1  MeV  (7$),   plus  one  prompt  proton  which  has 
no  significance  from  a  shielding  point  of  view.  (11) 

The  final  material  activation  reaction  to  be  consider- 
ed in  this  section  will  the  formation  of  tritium,  the  reac- 
tion of  interest  being: 

H2  +  n  >  H3  +  7  (3.3.2.1.3) 

The  thermal  neutron  cross  section  for  this  reaction  is  0.57 
milli-barns,  and  the  resulting  prompt  gamma  has  an  energy 
of  approximately  6.26  MeV  (9). 

3.3.2.2  Al27  (n,7)  Al28 

In  order  to  determine  the  dose  rate  on  the  reactor 

28 
top  due  to  the  decay  of  Al   ,  the  reactor  core  area  was 


! 


divided  into  three  regions  as  shown  by  the  dotted  lin< 

in  Figure  3.3.2.2.1.   Gamma  rays  born  due  to  the  decay 

28 
of  Al   were  assumed  to  be  point  sources  located  at  the 

center  of  each  of  the  three  regions.   The  volume  of 
aluminum  in  each  of  the  three  regions  was  determined  based 
on  a  best  estimate  of  utlimate  core  and  supporting  struc- 
ture composition: 

Region  I  V  =  5500  cm3 

Region  II  V  =  56,700  cm3 

Region  III  V  =  10,100  cm3 

The  aluminum  volumes  by  region  were  based  on  assumed  alu- 
minum volume  percents  within  each  region  of  2$  for  region 
I  above  the  core,  43$  for  region  II  in  the  core,  and  4.3$ 
for  region  III  below  the  core,  plus  a  1/8  in.  thick  100$ 
aluminum  tank  wall  forming  the  boundary  of  each  region. 
The  average  thermal  flux  in  each  of  the  three  regions  was 
conservatively  estimated  based  on  the  preliminary  core 

physics  calculations  described  in  Chapter  4  of  this  thesis: 

12  2 

Region  I  0   =  10   neutrons/cm  -sec 

13  2 

Region  II  0  =   5  x  10   neutrons/cm "-sec 

l4  2 

Region  III  0  =   10   neutrons/cm  -sec 

Only  the  thermal  neutron  induced  flux  is  considered 
here  due  to  the  rapid  decrease  in  the  reaction  cross  sec- 
tion of  aluminum  at  higher  energies.   The  activation  cross 
section  for  0.2  MeV  neutrons  is  down  to  about  10   barns, 


FIGURE     3.3.2.2.1 
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for  example.   Assuming  an  irradiation  of  infinite  duration, 

27 
the  number  of  neutron  captures  per  second  in  Al  '  is 

described  by  the  equation: 

Captures/sec  =  N  <Tact0VA1  (3.3.2.2.1) 

Proceeding  with  the  calculation  for  region  I, 

N<ract0VA1  =  (0.06023) (0.215) (1012) (5500)  = 

12 
7.12  x  10   captures/sec. 

Similarly  for  region  II, 

N<r  4-0V/U  =  3.66  x  10   captures/sec, 

aC  13   .A.-L 

and  for  region  III, 

N(r  ct^vA1  =  1-31  x  10   captures/sec. 

Information  for  Table  3.3.2.2.1  was  taken  from 
references  (7)  and  (9)>  and  lists  the  number  and  mean 

energy  of  the  gamma  rays  emitted  per  thermal  neutron  cap- 

27 
ture  in  Al   ,  along  with  the  gamma  flux  to  dose  rate  con- 
version and  the  decade  lengths  in  DpO,  iron  and  lead 
respectively. 

Equation  3.3.1.1  was  used  to  obtain  the  dose  rate  on 
the  reactor  top  directly  above  the  core  due  to  the  thermal 
neutron  activation  of  Al   .   Table  3.3.2.2.2  lists  the  in- 
put data  not  previously  given  required  to  solve  the  equa- 
tion, and  the  resulting  dose  rate  on  the  reactor  top  due 
to  gammas  from  each  region.   The  calculational  procedure  was 
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exactly  the  same  as  the  sample  calculation  done  in  Section 
3.3.1,  and  so  is  not  repeated  here.   Summing  the  dose 
rates  listed  in  Table  3.3.2.1.1  gives  a  total  dose  rate  on 
the  reactor  top  due  to  the  thermal  neutron  activation  of 
Al2^  to  be  0.06k   rarem/hr. 

3.3.2.3   016  (n,-p)   Nl6 

Although  the  decay  of  N   often  results  in  high  radi- 
ation levels  around  exit  coolant  piping,  there  is  normally 
little  reason  to  consider  the  0   (n,p)  N   reaction  for 
purposes  of  reactor  vessel  shielding  due  to  the  small 
activation  cross  section  of  N  .   However,  the  fluid  flow 
scheme  for  the  High  Flux  MITR  calls  for  the  coolant  to 
circulate  up  to  very  near  the  top  of  the  necked  portion  of 
the  D20  tank.   It  is  for  this  reason  that  the  dose  rate  due 
to  N   decay  is  calculated. 

The  threshold  energy  for  the  0   (n,p)  N  '  reaction  is 
very  high  (approximately  10  MeV) .   Therefore,  it  is  reason- 
able to  assume  that  N   may  only  be  formed  on  the  first 
collision  of  a  fission  neutron.   The  probability  of  a  fis- 
sion  neutron  interacting  to  form  N   is  simply  the  ratio  of 
the  cross  section  for  the  formation  of  N  'to  the  total 
removal  cross  section  seen  by  a  fission  neutron.   Table  1.1 
in  Rockv/ell  gives  the  effective  removal  cross  section  for 
fission  neutrons  in  DpO  as  2.8  barns  per  molecule  ( 12) . 
The  probability  that  an  0    (n,p)  N  '  reaction  occurs  in 
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the  first  collision  of  a  neutron  is  therefore, 

-6 
P(n,p)  =  278 =  7*15  x   10 

where  it  has  been  assumed  that  all  neutrons  have  that 

2     16 
first  collision  with  either  H  or  0   .   This  assumption 

is  conservative  and  leads  to  a  slight  overestimation  of 

P(n,p).   The  rate  of  production  of  N   in  the  entire 

reactor  becomes : 

r  =  $F   fission  neutrons  x  P(n,p) 

#  fission  neutrons  =  (5  x  10  watts) 

n  ih.   v  in10  fissionsUo  2lo  neutrons . 
(J. 14  x  10   watt-secH^^o  THUon^ 

^  fission  neutrons  =  3,9  x  10  neutrons /sec. 

R  =  (3.9  x  1017)(7.15  x  10"6)  =  2.78  x  1012  activations/ 
sec.   It  is  now  assumed  that  all  the  N  '  atoms  formed  are 
formed  in  the  DpO  which  circulates  up  through  the  necked 
core  tank.   This  assumption  is  conservative  because  some 
N   will  be  formed  in  the  reflector  DpO,  but  the  N   formed 
in  the  reflector  will  make  a  negligible  contribution  to  the 
dose  rate  on  the  reactor  top.   The  volume  of  the  DpO  in  the 
inner  coolant  tank  is  about  225  gallons.   Therefore,  at 
1850  gallons  per  minute  the  coolant  takes  7.5  seconds  to 
pass  from  the  inner  tank  inlet  to  exit.   The  external  sys- 
tem DpO  volume  is  approximately  475  gallons.   Therefore,  at 
1850  gallons  per  minute  the  total  cycle  type  for  the  D?0 
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is  approximately  23.3  seconds.   Nov/  assuming  N  '  to  be 
formed  uniformly  throughout  the  height  if  the  inner  core 
tank  (an  obviously  conservative  assumption  since  few 
fission  neutrons  will  undergo  their  first  collisions  near 
the  top  of  the  necked  D20  column),  and  taking  the  activity 
of  the  entire  volume  of  DpO  in  the  inner  tank  to  be  that  of 
the  coolant  as  it  exits  the  upper  portion  of  the  inner  tank, 
(again  a  very  conservative  assumption)  the  saturation 
activity  due  to  N  '  decay  becomes : 

"Xtc 
Asat  =  R    -Att  (3.3.2.2.1) 

1-e    * 

where   X  =  decay  constant  for  N   =  —IS  00„  =  .09-1)3  sec 

f  .  o _)   sec 

t  =  irradiation  time  in  the  inner  core  tank 
c 

t,  =  total  cycle  time 

12   i-e-(0-09l|3H7-5) 
Asat  =  (2*68  x  10   )   1]e-(0.0943)(23.3) 

12 
A  sat  =  1.59  x  10  "  disintegrations /sec.   From  Section 

3.3.2.1,  it  is  recalled  that  N   decays  by  emission  of 

either  a  6.1  MeV  gamma  (75$  of  the  time)  or  a  7.1  MeV 

gamma  (7$  of  the  time).   It  is  assumed  in  the  following 

calculations  that  N   always  decays  by  emission  of  a  6 

MeV  gamma.   (This  assumption  is  conservative  since  in 

lead  a  6  MeV  gamma  is  more  penetrating  than  a  7  MeV  gamma. ) 

The  done  rate  on  the  reactor  top  due  to  N"  decay  is 
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calculated  assuming  the  total  saturated  activity  to  be  a 
point  source  located  at  the  surface  of  the  DJO   column. 
Again  substituting  into  Eq.  3.3.1.1: 

0  =   (i.59*lQ12)(l.O) (i.2+i4.6)  io-(-26+6-7i) 

4ir   (5+.  083+1)  30. 48  d 

0  =  6.32  y's/cm  -sec  . 

Therefore, 

Dose  Rate  =  0.04  mrem/hr. 

16 
The  dose  rate  due  to  N   decay  calculated  above  is 

ultra-conservative  from  a  shielding  point  of  view.   (The 

dose  is  overestimated  by  about  a  factor  of  10.)   The 

primary  purpose  of  the  calculations  was  to  show  that  even 

using  overly  conservative  assumptions  the  decay  of  N 

would  not  make  a  dominant  contribution  to  the  total  dose 

rate  on  the  reactor  top. 

3.3.2.4  H2  (n,y)  H3 
Since  only  thermal  neutrons  are  important  in  the 
formation  of  tritium,  the  gamma  ray  flux  resulting  from 
tritium  formation  will  vary,  as  a  function  of  position  in 
the  core,  directly  with  the  thermal  neutron  flux.   Tritium 
decays  by  emission  of  a  beta  particle,  which,  although 
very  significant  from  a  health  physics  standpoint  should 
it  be  ingested  into  the  body,  is  insignificant  for  purposes 
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or  reactor  top  shielding  calculations.   For  calculational 
purposes,  the  DpO  volume  was  divided  into  four  regions  as 
shown  in  Figure  3.3.2.4.1,  and  the  gamma  rays  from  each 
region  were  assumed  to  originate  from  point  sources  loca- 
ted as  shown  in  the  figure.   An  average  thermal  flux  was 
assumed  for  regions  I  and  II  based  on  a  conservative 
interpretation  of  the  results  shown  in  Chapter  4  of  this 
thesis.   The  average  thermal  flux  in  regions  III  and  IV 
was  calculated  by  considering  the  necked  core  tank  to  be 
a  DpO  thermal  column  in  which  the  calculational  method  to 
be  shown  in  Section  3.3.4.2  applies. 

The  number,  n,  of  gamma  rays  produced  per  second  in 
each  region  due  to  the  formation  of  tritium  is  given  by: 

n-N<j;ct^O'  (3.3.2.4.1) 

24  3 

where  N  =  deuterium  atom  density  =  0.06646  x  10   atoms/cm  , 

2      3 
or  +  =  activation  cross  section  for  H  (n, 7)11  reaction  = 

0.57  x  10"3  barns, 

0   =  average  thermal  flux, 

V^  q  =  coolant  volume  of  interest. 


FIGURE    3.3.2.4.1 
REGIONS  FOR  H2  (n.tf)  H3 


96 


CALCULATIONS 


ET 


* 


in 

sly 
7K 


n 

— ,       ^»^ 

'N 

COR 

E 

1 

/IN 

— — ~ 

97 


Table  3.3.2.4.1 


3  2 

Re  r;ion   Volume  (cm  )    (neutrons /cm  -sec)    S( gammas /see) 


I 

6.0  x 

io5 

14 

1.0  x  10^ 

1.14  x 

lO1* 

II 

5.4  x 

l05 

1? 
1.0  x  1CT^ 

1.02  x 

lb13 

III 

1.95  x 

io5 

1.8  x  1010 

6.65  X 

1010 

IV 

1.95  x 

105 

5.3  x  106 

1.96  x 

107 

Again  the  calculational  procedures  outlined  in  Section 

3.3.1  are  followed  in  order  to  obtain  the  dose  rate  on  the 

2        3 
reactor  top  due  to  the  H   (n,?)  H  reaction.   Since  a  6  MeV 

gamma  is  more  penetrating  in  lead  than  a  7  MeV  gamma,  the 
reaction  gamma  rays  are  taken  to  be  6  MeV  (a  conservative 
assumption).   Table  3.3.2.4.2  lists  the  results  of  the  cal- 
culations for  gamma  flux  and  dose  rate  on  the  reactor  top 

2        3 
due  to  the  H  (n, 7)  H  reaction. 

Table  3.3.2.4.2 

Region (gammas /cm  -sec) Dose  Rate  (mrem/hr) 

I  0.246  1.56  x  10"3 

II  0.026  I.65  x  10""^ 

III  0.0056  3.5^  x  10"5 


IV  0.000036  2.30  x  10 


-7 


Total  O.278  I.76  x  10"3 


As  shown  by  the  above  calculation,  the  dose  rate  on  the  reac 
tor  top  due  to  tritium  formation  is  negligible. 
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3.3.3  Fission  Product  Decay 

As  one  might  expect,  the  dose  rate  due  to  fission 
product  decay  is  small  compared  with  the  dose  rate  from 
prompt  fission  gammas.   An  infinite  irradiation  time  is 
assumed,  thereby  yielding  a  dose  rate  which  contains  the 
contribution  of  both  short  and  long-lived  fission  products. 
The  results  so  obtained  are  conservative  since  in  practice 
the  reactor  will  operate  during  the  week  and  be  shut  down 
on  weekends. 

According  to  Rockwell,  the  rate  of  energy  emission  by 
gamma  rays  after  an  infinite  operating  time  at  5  IV±W  is  1.14 
x  10   MeV  per  second  (12.)*   If  "the  average  decay  gamma  ray 
energy  is  taken  to  be  2  MeV,  the  decay  gamma  ray  source 
strength  becomes  5.7  x  10  '  gammas/sec  ( 13J .   A  point 
source  is  assumed  at  the  core  center,  and  the  calculational 
procedure  outlined  in  Section  3.3.1  is  once  again  used.   The 


dose  rate  on  the  reactor  top  due  to  fission  product  decay 

i       -4 
.lculated  to  be  4.1  x  10   n 

3.3.4  Neutron  Attenuation 


i       -4 
is  calculated  to  be  4.1  x  10   mrem/hr. 


3.3.4.1  Fast  Fission  Neutrons 
The  relation  for  the  attenuation  of  fission  spectrum 
neutrons  is  given  by  Eq.  3.2.3.1,  here  repeated  in  slightly 
different  form: 

Q  =  -^e        l      ,  (3.3,4.1.1) 

47TIT 
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where  S  =  neutron  source  strength  (neutrons /sec) , 
R  =  distance  from  the  source  (cm), 
x  =  thickness  of  attenuating  materials  (cm),  and 
A  =  relaxation  length  for  removal  of  neutrons  (cm). 

The  relaxation  length,  X  ,  is  calculated  by  the  relation: 

Xr  -  Nib-  (cm)  ,  (3.3.4.1.2) 

a   r 

where   Q"*  =  microscopic  removal  cross  section  (barns/atom), 

23 
N  =  Avagardo's  number  (6.04  x  10   molecules/gm- 

mole) , 

3 
p  =  density  (gms/cm  ),  and 

A  =  molecular  weight  (gms/gm-moie) . 

For  Do0, 

Xr  =  P,  20-°3 pip-  =  10.7  cm. 

r   (0.6o4xl0^4) (1.105) (2.8  x  10"^4) 

The  relaxation  lengths  for  lead  and  iron  are  given  in  Table 
4.2.3.1  to  be  8.94  cm  and  5.90  cm  respectively. 

As  calculated  in  Section  3.3.1  there  will  be  1.55  x 

17 
10   fissions/second  in  the  High  Flux  MITR  operating  at  5 

MW.   Taking  the  neutron  yield  per  fission  to  be  2t47,  and 

considering  that  one  of  these  neutrons  is  needed  to  sustain 

the  chain  reaction,  there  will  be  a  total  of  (1.47) (1.57  x 

17  \  17 

10"')  2.31  x  10"  '  fission  neutrons  per  second  available  for 

leakage  or  for  absorption  in  the  various  reactor  materials. 
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The  conservative  approach  is  to  assume  that  all  2.31  x 

17 
10  '  neutrons  leak  from  the  core. 

Substituting  into  Eq.  3.3.4.1.1,  the  fast  fission 

neutron  flux  on  the  reactor  top  is  calculated: 

iv  ,9x30.48  ,  .083x30.48   1x30. 48N 

0   = (2.31xl01Y) e'^  IU7T'  57*90   +  b.9+  } 

47r{]; 9+5+. 083+1)  x  30.48}  2 

-2  2 

0  =  I.69  x  10"  neutrons /cm  -sec. 

From  Figure  2.3  of  reference  (12),  the  conversion  from  fis- 
sion neutron  flux  to  dose  rate  is  taken  to  be  8.6  neutrons/ 
cm  -sec  per  mrem/hr.   Therefore,  the  dose  rate  on  the  reac- 
tor top  due  to  fast  fission  neutrons  is  calculated  to  be, 
-2 
fcjX61Q —  =  1'9^  x  10"3  mrem/hr.   As  one  might  expect,  the 

low  fission  neutron  flux  on  the  reactor  top  is  primarily  due 
to  the  excellent  moderating  characteristics  of  the  DpO  in 
the  necked  core  tank  region.   Using  the  above  calculational 
method,  it  was  found  that  the  fast  neutron  flux  at  the  cool- 
ant  surface  is  only  I.96  neutrons/cm  -sec. 

Not  included  in  the  above  calculation  is  the  moderating 
effect  of  the  masonite.   In  practice,  most  of  the  fast  neu- 
trons reaching  the  masonite  are  thermalized  in  the  masonite, 
and  subsequently  absorbed  in  the  surrounding  boral.   To 
neglect  the  masonite,  therefore,  yields  a  conservative 
result. 
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3.3.4.2  Thermal  Neutrons 
In  order  to  determine  the  dose  rate  on  the  reactor  top 
due  to  thermal  neutrons,  it  is  first  necessary  to  determine 
the  thermal  neutron  flux  at  the  surface  of  the  necked  coolant 
tank.   This  can  be  done  by  assuming  the  necked  coolant  tank 
region  to  be  a  cylindrical  thermal  column.   The  thermal  flux 
at  the  base  of  the  necked  coolant  tank  is  conservatively 
estimated,  based  on  the  results  of  the  flux  calculations  in 

IP  2 

Chapter  4  of  this  thesis,  to  be  10  neutrons /cm  -sec.  For 
the  situation  of  interest  here,  the  applicable  relation  for 
the  flux  as  a  function  of  axial  position  in  the  necked  core 
tank  is : 

0(7,)    =  0oe"7Z  (3.3.4.2.1) 

where  Z  =  axial  distance  above  the  necked  core  tank  base  (cm), 


t 


7   =l/(i)2  +  (%^)2  (cm"1)  , 


L'     v  r 

L  =  diffusion  length  =  100  cm  in  DpO, 

r1  =  extrapolated  tank  radius  =  r  +  0.71  )\  +  . 
XfY   =  transport  mean  free  path  =  3D 

D  a  diffusion  length  =  0.80  cm  in  DpO,  and 
r  =  core  tank  radius  =  25.4  cm. 
Substituting  in  numerical  values, 


i? 


■>  =   1/(150^  +  (f$£)  =  0.0842  cm' 
0(Z)  =  (io12)e-(0-08i;2MZ>  neutrons/cm2-sec. 
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Calculated  values  of  0(Z)  are  given  in  the  following  table: 

Table  3.3.4.2.1 

Z  (cm  above  base  of  P 

necked  core  tank) (Z)  neutrons/cm  -sec 

48  1.8  x  1010 

145  5.3  x  106 

193  (coolant  surface)        9.6  x  10 

Because  lead  is  a  very  poor  attenuator  of  thermal  neutrons, 
and  because  it  is  desirable  to  avoid  the  capture  gamma  rays 
which  are  produced  in  a  Fe^   (n,7)  Fe^ '   reaction,  it  was 
decided  to  line  the  interior  of  the  refueling  enclosure  with 
boral.   A  1/4  in.  thick  sheet  of  boral  will  result  in  a 
maximum  fractional  transmission  of  thermal  neutrons  of  5  x 
10   (l4)»   Consequently,  if  the  refueling  enclosure  is 
lined  with  two  1/4  inch  sheets  of  boral,  the  thermal  neutron 
flux  will  be  effectively  eliminated  before  entering  the 
iron  or  lead  shielding. 

3.3.5  Photodis integration  of  Deuterium 

2        1 
The  threshold  energy  for  the  H  {y,n)   H  reaction  is 

2.2.3  MeV.   The  reaction  cross  section  increases  from  zero 

at  the  threshold  energy  to  a  maximum  value  of  2.24  millibarns 

for  a  gamma  energy  of  4  MeV.   The  cross  section  then  decreases 

to  a  value  of  1.2  millibarns  for  a  13  MeV  gamma  ( j-1 ) .   For 

calculational  purposes,  it  will  be  assumed  that  all  gammas 
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of  energy  greater  than  2  MeV  are  capable  of  initiating  the 
(y,n)  reaction  and  that  the  average  reaction  cross  section 
is  1.5  millibarns. 

In  order  to  determine  the  number  of  neutrons  generated 
as  a  function  of  position  in  the  core,  it  is  necessary  to 
determine  the  flux  of  gamma  rays  above  2  MeV  as  a  function 
of  position  in  the  reactor  tank.   For  this  purpose,  the  tank 
was  again  divided  into  four  regions  as  shown  in  Figure 
3.3.2.3.1.   The  gamma  flux  in  a  given  region  is  that  due  to 
gammas  born  in  the  region  plus  that  due  to  gammas  entering 
the  region  from  each  of  the  surrounding  regions.   Except 
for  fission  gammas,  the  number  of  gamma  rays  per  region  of 
energy  greater  than  2  MeV  was  taken  from  the  results  of 
calculations  performed  in  previous  sections  of  this  chapter. 
Table  3  on  page  50  of  reference  (7)  shows  that  19.3$  of  all 
fission  gammas  have  energies  greater  than  2  MeV,  and  so  this 
percentage  was  used  rather  than  the  assumed  percentage  given 
in  Section  3.3.1.   The  flux  entering  a  region  from  each  of 
its  three  neighbors  is  calculated  using  Eq.  3.3.1.1  where 
the  gamma  rays  born  in  each  region  v/ere  treated  as  point 
sources,  and  attenuation  to  neighboring  regions  is  based  on 
point  source  to  point  source  distance.   Mean  gamma  energy 
is  taken  to  be  3.5  MeV  (a  very  conservative  estimate). 
Results  obtained  are  listed  in  the  following  table: 
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Having  determinedthe  average  gamma  ray  flux  per 
region,  the  number  of  neutrons,  n,  produced  per  second  is 
given  by  the  relation: 

n  =  N(T0  V  neutrons/sec,  (3.3.5.1) 

where 

?  22  3 

N  =  atom  density  of  H  =  6.6  x  10   atoms /cm  , 

-27    2 
(T  =   activation  cross  section  =  1,5  x  10   '  cm  , 

0     =  average  region  gamma  flux 

V  =  volume  of  D^O  in  the  region. 

Table  3.3.5.2 

Region Photoneutrons,  n,  born  per  second 

14 
I  1.09  x  10 

14 
II  1.07  x  10 

III  8.4  x  lo11 

IV  1.25  x  1010 

2         1 
The  neutron  and  proton  released  by  the  H  (y,n)  H  reac- 
tion share  equally  the  difference  between  the  gamma  energy 
and  the  threshold  energy.   Assuming  an  average  reaction 
gamma  energy  of  3.5  MeV,  the  average  energy  of  the  neutron 
released  is  0.64  MeV.   Conservatively,  the  energy  spectrum 
of  the  photoneutrons  born  in  the  coolant  tank  can  be 
approximated  as  a  fission  spectrum.   The  neutron  flux  on 
the  reactor  top  is  therefore  calculated  by  substituting 
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into  Eq.  3.3.4.1.1.   To  the  results  obtained  from  Eq. 


I 

2, 

,8  x  10"13 

II 

6. 

-4 
,5  x  10  ^ 

III 

0.03 

IV 

0.09 

Total 

0.12 

3.3.4.1.1,  which  included  neutron  attenuation  due  to  dis- 

_2 
tance,  D20,  iron  and  lead  was  multiplied  a  factor  of  10 

to  account  for  the  attenuation  due  to  the  6  in.  of  masbn- 

ite  (16).   The  calculated  results  are  summarized  in  the 

following  table: 

Table  3.3.5.3 

Region Neutrons  /cm  ~  -s  ec Dose  Rate  (mrem/hr) 

3.3  x   10"14 

7.6  x  10"5 

0.004 

0.010 

0.014 

As  the  above  tabulated  values  show,  the  addition  of  6 
in.  of  masonite  is  sufficient  to  reduce  the  neutron  dose 
rate  on  the  reactor  top  to  a  very  low  level.   From  a  dose 
rate  point  of  view,  it  would  be  possible  to  reduce  the  thick- 
ness of  the  masonite.   However,  since  masonite  is  relatively 
inexpensive,  the  reduction  of  neutron  induced  activation  of 
reactor  top  structural  materials  resulting  from  the  addi- 
tional masonite  shielding  is  well  worth  the  small  addition- 
al expense.   A  low  activity  level  is  particularly  desirable 
in  planning  for  the  eventuality  that  one  may  want  at  some 
future  time  to  disassemble  the  reactor  top  shielding  enclo- 
sure. 
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3 . 4  Shielding  Calculations,  Reactor  Shutdown  (Phase  I) 

3.4.1  General 

The  predominate  radiation  sources  in  the  reactor  after 
shutdown  are:   (1)  fission  product  decay  gammas,  (2)  gammas 
from  the  induced  activity  of  the  structure  and  the  coolant, 
and  (3)  the  photodis integration  of  deuterium.   Immediately 
after  reactor  shutdown,  radiation  due  to  delayed  neutrons 
should  also  be  considered. 

Quite  obviously,  the  magnitude  of  the  shutdown  radiation 
sources  is  only  a  small  fraction  of  that  for  the  reactor 
operating  at  power.   Consequently,  a  calculation  of  the  dose 
rate  on  the  reactor  top  after  shutdown  would,  under  normal 
circumstances,  yield  a  minute  result  which  would  have  little 
meaning.   There  are,  hov/ever,  two  important  situations  in 
which  the  dose  rate  on  the  reactor  top  due  to  shutdown  radi- 
ation may  be  quite  significant:   (1)  when  reactor  top 
shielding  is  removed  from  the  area  directly  above  the  core 
tank,  and  (2)  immediately  after  shutdown  when  all  coolant  is 
lost  from  the  necked  core  tank. 

3.4.2  Reactor  Ton  Shield  Plug  Removed 

Shielding  directly  above  the  reactor  core  tank  would 
not  normally  be  removed  except  for  purposes  of  reactor  main- 
tenance or  repair.   Assuming,  for  instance,  that  the  reactor 
is  shut  down  at  5  p.m.  on  Friday,  work  in  the  area  of  the 
reactor  top  would  not  commence  until  Saturday  morning,  a 
cooling  time  of  at  least  14  hours.   During  this  time  period, 
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28 
virtually  all  the  activity  due  to  Al   (Ti/p  =  2«3  rain) 

and  N   (Ti/p  ~  7. 35  sec)  decay  will  have  died  away. 

From  Figure  3  of  reference  (13),  the  total  gamma  energy 

decay  rate  after  infinite  operation  followed  by  a  Ik   hour 

decay  period  is  approximately  1.5  x  10  '  MeV/watt-sec. 

16 

Following  operation  at  5  MW  this  figure  becomes  7.5  x  10 

MeV/sec.   Again  taking  the  average  gamma  energy  to  be  2 
MeV,  the  fission  product  decay  source  strength  is  3.75  x 
10   gammas/sec. 

Equation  3.3.1.1  is  again  used  to  determine  the  gamma 
flux  on  the  reactor  top,  where  in  this  case  gamma  attenua- 
tion is  due  only  to  distance  and  the  coolant  in  the  necked 
core  tank.   The  calculated  result  is  0  -   1.03  x  10  gammas/ 
cm  -sec,  which  converts  to  a  dose  rate  of  303  mrem/hr. 

In  order  to  determine  the  number  of  photoneutrons  pro- 
duced by  the  fission  product  decay  gammas,  it  is  necessary 
to  determine  the  flux  of  gammas  in  excess  of  2.23  MeV  as  a 
function  of  position  in  the  coolant  tank.   Figure  3.2b  of 
reference  (12)  shows  that  gammas  in  excess  of  2  MeV  comprise 
roughly  3°/o   of  the  total  decay  gamma  source  strength  l4  hours 

after  shutdown.   The  decay  gamma  source  strength  capable 

2        ] 
of  initiating  the  H  (7>n)  H"  reaction  is  therefore, 

J0^> (1-5  *  1016  MeV/Wj.  =  g-8  x  1Qlk   garamas/sec, 
2.3  MeV/ gamma         -/         & 
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The  core  tank  was  divided  into  three  regions  (shown  in 
Figure  3.4.2.1)  and  the  average  gamma  flux  per  region 
calculated  using  Eq.  3.3.1.1.   Knowing  the  average  flux 
capable  of  inducing  the  (r,n)    reaction,  the  number  of  photo- 
neutrons  produced  in  each  region  was  calculated  using  Eq. 
3.3.5.1*  where  for  fission  product  decay  gammas  the  average 
activation  cross  section  was  taken  to  be  0.001  barns. 
Results  obtained  are  summarized  in  the  following  table: 

Table  3.4.2.1 

Gamma  Flux           Neutrons  Produced 
Region ( gammas /cm^ -sec) (neutrons /sec) 

I  1.87  x  1010  1.41  x  1012 

II  2.23  x  108  2.88  x  109 

III  7.92  x  105  1.02  x  107 

Taking  the  average  activation  gamma  ray  energy  to  be  2.3 
MeV,  the  average  energy  of  the  photoneutron  is  calculated 
to  be: 

2.30  -  2.23  =  35  KeV 

Neutrons  of  such  low  energy  will  be  thermalized  very  near 
their  birthplace  by  the  surrounding  D^O.   The  photoneutrons 
generated  due  to  decay  gammas  are  therefore  considered  to 
be  thermal  neutrons,  and  the  method  of  Section  3.3.4.2  is 
used  to  determine  the  neutron  flux  on  the  surface  of  the 
core  tank.   Having  reached  the  surface,  the  neutrons  are 
assumed  to  radiate  as  a  point  source  out  through  the  opening 
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FIGURE   3.4.2. 
REGIONS  FOR  PHOTO-NEUTRON  CALCULATIONS 


Ill 

in  the  reactor  top  shielding.   The  dose  rate  on  the  reactor 
top  due  to  photoneutrons  is  calculated  to  be  0.03  mrem/hr, 
negligible  in  comparison  to  the  dose  rate  due  to  decay 
gammas . 

3.4.3  Coolant  Lost  from  Necked  Tank 

Immediately  after  reactor  shutdown  coincident  with  the 
loss  of  coolant  in  the  necked  core  tank,  the  primary  sources 

of  radiation  would  be  the  decay  of  fission  products  in  the 

28 
core  and  the  decay  of  Al   in  the  core  and  surrounding 

structure.   From  Section  3.3.2.1,  the  initial  Al  ' '  rate  of 

16 
decay  will  be  4.97  x  10   gammas  per  second,  the  decay 

gammas  having  an  energy  of  1.8  MeV  (taken  as  2  MeV  for  cal- 
culations).  The  total  source  strength  immediately  after 

17 
shutdown  is,  therefore,  6.2  x  10   gammas  per  second.   The 

gamma  flux  on  the  reactor  top  is  calculated  using  Eq. 

3.3.1.1,  where  in  this  case  gamma  attenuation  is  due  only 

to  distance  and  the  iron  and  lead  shielding  on  the  reactor 

top.   The  calculated  result  is  0  =   1.3  x  lOr    gammas /cm  - 

sec,  which  converts  to  a  dose  rate  of  383  mrem/hr.   After 

14  hours,  should  coolant  still  be  absent  from  the  necked 

core  tank,  the  dose  rate  on  the  reactor  top  will  have 

decayed  to  approximately  23.2  mrem/hr. 

17 
In  Section  3.3.4.1  it  was  assumed  that  2.31  x  10 

neutrons  per  second  leak  from  the  core  during  normal 

operation.   Assuming  the  delayed  neutrons  fraction  to  be 
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0.007,  1.63  x  10   of  these  neutrons  are  delayed,  neutrons. 
Therefore,  at  the  instant  of  reactor  shutdown  the  core 
will  continue  to  emit  1.62  x  10  ^  neutrons/sec.   Using  the 
method  of  Sections  3.3.4.1  and  3.3.5,  the  neutron  dose  on 
the  reactor  top  at  the  instant  of  shutdown  is  calculated, 
where  in  this  case  neutron  attenuation  is  provided  only  by 
distance  and  the  masonite,  iron  and  lead  shielding.   The 
result  is  18.5  rem/hr.   Fortunately,  the  half-life  for 
delayed  neutron  emission  is  small  (0.43  -  55.6  seconds)  and 
so  after  about  15  minutes  the  dose  rate  on  the  reactor  top 
would  be  negligible  (15.). 

Needless  to  say,  should  coolant  be  lost  from  the 
necked  core  tank.,  personnel  should  stand  well  clear  of  the 
area  immediately  above  the  necked  core  tank. 
3 . 5  Dose  Rate  Through  Concrete  Shielding  (Phase  I) 

Thus  far  in  this  chapter  all  shielding  calculations 
have  been  performed  for  radiation  passing  up  through  the 
necked  coolant  tank  and  out  through  the  reactor  top  shield- 
ing enclosure.   Although  one  can  reasonably  expect  that  the 
6-1/2  feet  of  dense  concrete  which  shields  the  reactor  top 
work  area  surrounding  the  necked  core  tank  will  allow  only 
minute  radiation  levels  on  the  reactor  top,  one  brief  and 
approximate  calculation  is  in  order.   An  estimation  of  th< 
required  density  for  the  reactor  top  shielding  concrete 
will  also  be  made. 
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All  radiation,  except  the  thermal  neutron  flux,  is 
assumed  to  be  emitted  from  a  point  source  located  at  the 
core  center.   Radiation  source  levels  as  a  function  of 
energy  were  taken  from  the  preceeding  sections  of  this 
chapter. 

Table  3.5.1 

Gamma  Ray  Source  Strength 
Energy  (MeV) (gammas /sec) 

8  1.74  x  1016 


6  4.56  x  10 


7  1.14  x  10" 

,16 


4  3.80  x  1016 

3  3. 07  x  1017 

2  6.50  x  1017 

"I  R 
1  1.17  x  10 

17 
Fast  neutron  source  strength  =  2,31  x  10   neutrons/sec, 

12 
Thermal  neutron  flux  at  the  tank  wall  =  10   neutrons/ 

p 
cm  -sec.   Gamma  ray  attenuation  through  the  D^O  and  concrete 

shielding  is  calculated  utilizing  Eq.  3.3.1.1  where  for  the 

current  case: 

q  -(D+D  ) 

0  = o  (Bn+Bp)10    D  C   ,  (3.5.1) 

4tt(rd+rc)2   D  C 


where  the  subscripts  D  and  C  refer  to  D^O  and  concrete 
respc  tively.   Decade  lengths  as  a  function  of  energy  were 
taken  from  reference  (7).   The  calculated  results  appear  in 
the  following  table: 
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Equation  3.3.4.1.1  is  used  to  determine  the  fast  neu- 
tron flux  on  the  reactor  top.   Converting  the  exponential 
attenuation  coefficients  to  decade  lengths , 

dD  Q  =  0.0754  x  10.7  =  0.806  ft., 

d      ,   =  0.8  ft.  (experimentally  measured  value  for 
c  one  rcDG 

O.W.R.  concrete  shielding,  specific 
gravity  =  4.6). 
Substituting  into  Eq.  3.3.4.1.1, 

0  =        2.31  x  1017       ]0-(3.72  +  8.13) 
kir   [(3+6.5)30.48]  2 

0   =  0.242  neutrons/cm  -sec  . 
therefore,  Dose  Rate  =  0.028  mrem/hr. 

Thermal  neutron  attenuation  through  the  concrete 
shielding  was  calculated  using  Eq.  3.2.3.2,  where  for  con- 
crete with  specific  gravity  =4.5  Gage  has  determined  that 
L  =  1.91  cm.,  and  D  =  0.372  cm.   Substituting  into  Eq. 
3.2.3.2  yields: 

0   =  2.58  0oe"x/1'91  (3.5.2) 

A  thermal  neutron  decade  length  in  concrete  is  equal  to 
0.0754  x  (L)  =  (0.0754) (1.91)  =  0.144  ft.   Therefore, 

0  =  2.58  0o10-x/0'1J|\    or  0  =    (2.58  x   1012)(10-x/0'144) 

0=  (io12-4l2)(io-6-94x) 
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In  order  to  reduce  the  thermal  neutron  flux  to  1  neutron/cm  - 

IP  4 12 
sec,  a  concrete  thickness  of  x  -   — g-'-qir—  =1.8  ft.  is  re- 
quired. Thus,  nearly  all  the  thermal  neutrons  from  the  tank 
will  be  absorbed  in  the  first  2  ft.  of  concrete.   Beyond  the 
first  2  ft.  the  thermal  neutron  flux  results  from  the  slowing 
down  of  the  fast  flux,  and  is  thus  approximately  equal  to 
the  fast  flux  at  each  point. 

Although  the  calculations  of  this  section  indicate  a 
dose  rate  through  the  concrete  shielding  of  0.03  mrem/hr, 
thereby  indicating  that  concrete  of  specific  gravity  =  4.5 
is  actually  more  dense  than  required,  it  is  felt  that  dose 
rate  uncertainties  attendant  with  the  various  penetrations 
through  the  upper  concrete  shielding  are  sufficient  to 
warrant  the  use  of  the  s.g.  =  4.5  concrete. 
3.6  Shielding  Calculations  Summary  (Phase  I) 

3.6.1  Dose  Through  Shielding  Directly  Above  the  Core 
Source  Dose  Rate  (mrem/hr) 

Core  Prompt  Fission  Gammas  0.068 

Al27  (n,y)   Al28  0.064 

N   Decay  0.04 

H2  (n,7)  II3  0.0018 

Fission  Product  Decay  Gammas  0.0004 

Fast  Fission  Neutrons  0.002 

Thermal  Neutrons  0.005 

Photoneutrons  0.014 

Total       0.195  mrem/hr. 
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3.6.2  Dose  Through  Reactor  Top  Concrete 

Source  Dose  Rate  (mrem/hr) 

All  Sources  Combined  0.0018  mrem/hr 

3.6.3  Shielding  Conclusions 

The  maximum  dose  rate  anywhere  on  the  High  Flux  MiTR 
top  is  calculated  to  be  less  than  0.2  mrem/hr.   This  fig- 
ure is  a  factor  of  5  less  than  the  design  goal,  and  a  fac- 
tor of  12.5  less  than  the  limit  set  by  the  AEC.   The 
reactor  top  shielding,  as  described  in  previous  sections 
of  this  chapter,  is  therefore  considered  adequate  and  safe. 
3 • 7  Shielding  Calculations,  Normal  Operation  (Phase  II) 

3.7.1  General 

The  phase  II  design  differs  from  the  phase  I  design, 
for  shielding  purpose.",  in  three  basic  ways:   (1)  the 
water  column  above  the  core  is  light  water  instead  of  heavy 
water,  (2)  the  water  column  above  the  core  is  11  ft.  high 
vice  8  ft.,  and  (3)  the  shielding  lid  on  top  of  the  water 
column  is  simply  a  flat  cover  instead  of  the  "doghouse" 
enclosure  shown  in  Figure  3.1.1.   In  going  from  phase  I  to 
phase  II,  the  magnitude  of  the  shielding  problem  was 
greatly  reduced.   No  longer  need  one  be  bothered  by  tritium, 
either  from  a  health  physics  or  <x   shielding  standpoint. 
Light  water  makes  a  much  better  neutron  shield  than  heavy 
water  due  to  its  much  larger  neutron  absorption  cross  sec- 
tion, and  in  addition  does  not  act  as  a  source  for 
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photoneutrons.   By  iterative  calculations,  it  was  deter- 
mined that  a  lead  shield  10  inches  thick  on  top  of  the 
water  column  would  provide  adequate  shielding.   The 
following  calculations  serve  to  justify  this  choice  and 
to  indicate  the  expected  dose  rate  on  the  reactor  top. 

3.7.2  Fission  lamma  Rays 

Table  3.7.2.1  presents  the  calculations  for  the  leak- 
age of  fission  gamma  rays  through  the  lead  shielding  lid  on 
the  tank  top. 

Table  3.7.2.1 

1  Mev       3  Mev         6  Mev 


Material  Thickness   D 

B 

D     B 

D      B 

H20       12  ft.    11.2 

170 

6.3    12.3 

4.4    5.4 

Lead      10  in.    8.8 

6 

5.2    6.4 

5.6   10.2 

Totals   20.0 

176 

11.5   18.7 

10.0   15.6 

0(r)  (no  attenuation)   1  x 

1012 

3.3  x  1011 

3  x  1010 

7's/cm  -sec 

0   (attenuated)      1.76  x 

10-6 

19.5 

46.8 

2 
7's/cm  -sec 

Conversion  Factor       70C 

) 

270 

158 

Dose  Rate  (mrem/hr)  2.5  x 

io9 

0.072 

O.292 

Total  Dose  Rate  -  0.364  mrem/hr 

3.7.3  Material  Activa'-.i on 

3.7.3.1  Al27  (n,y)   Al28 

Froij  Section  3.3.2.2  it  is  seen  that  the  only  signifi' 
cant  contribution  to  the  reactor  top  dose  rate  is  made  by 
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6  and  8  MeV  Gammas  from  regions  II  and  III.    Collapsing 

27  28 

the  Al  '  (n,?)  Al  source  strength  to  a  point  source  lo- 
cated at  the  core  center,  a  source  of  5  x  10"  gammas/sec 
is  obtained.  Calculated  results  are  given  in  the  follow- 
ing table. 

Table  3.7.3.1.1 


6  Mev  8  Mev 


Material 

Thickness 

D 

B 

D 

B 

H20 

12  ft. 

4.4 

5.4 

3.78 

4.0 

Lead 

10  in. 

5.6 

10.2 

5.75 

11.0 

Totals 

10.0 

15.6 

9.53 

15.0 

0(r)    (no  attenuation)    3.42  x  10   7's/cm  -sec 

Fraction  0.21  0.35 

o 

0   (attenuated)  7's/cm  -sec  11.2  53 

Conversion  Factor  158  128 

Dose  Rate  (mrem/hr)        0.071  0.415 

Total  Dose  Rate  =  0,485  mrem/hr 
3.7.3.2  H1  (n,y)  H2 
The  production  of  gamma  rays  from  thermal  neutron  cap> 
ture  in  hydrogen  was  not  treated  in  the  phase  I  calcula- 
tions, and  so  will  be  described  here.   One  gamma  ray  of 

1  2 

energy  2.23  MeV  is  emitted  in  the  II"  (n,y)   H  reaction, 

the  reaction  cross  section  being  0.33  barns  (7).   In  order 
to  simplify  calculations,  the  water  in  the  core  was  con- 
sidered to  be  a  sphere  1  ft.  in  radius  and  the  average 
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13 
thermal  neutron  flu  :  over  the  sphere  taken  to  be  5  x  10 

2 
neutrons /cm  -sec.   The  number  of  neutron  captures  by  hydro- 
gen was  calculated  by  the  following  equation,  which  is  also 
the  number  of  gamma  rays  emitted.   A  point  source  of 
gamma  rays  was  assumed  at  the  core  center  and  attenuation 
of  gamma  rays  through  the  water  and  lead  shielding  was 
calculated  as  before. 

neutron  captures/cm  -sec  =  NCT  0  =  2  0 

a     a 

neutron  captures/sec     =  VZ  0  =   gamrr.a  rays/sec 
N  =  6.71  x  1022  ^msH^ 


cm 


V  =  119,000  cm3 


2  (H)  =  0.0221  cm  -1 
a 

gamma  rays/sec  =  1.3  x  10  '  (assumed  to  be  2  MeV) 
The  results  of  the  dose  rate  calculation  are  given  in  the 
following  tabic: 

Table  3.7.3.2.1 

Material Thickness D B_ 

Water        12  ft.  7.85  24.0 

Lead         10  in.  5.7  6.2 

Totals    13.55  30.2 


0(E)    (no  attenuation)    =   8.9  x    1010    7fs/cm2-sec 
0(E)    (attenuated)  =  O.076   7's/cm2-sec 

Dose   Rate  =   2.23   x    10        mrern/hr 


I;1 1 


•  3.7.3.3   O16  (n,p)  N16 
In  the  phase  II  design  the  coolant  water  with  its 
N   activity  will  not  circulate  in  the  upper  4  ft.  of 
the  necked  core  tank.   Thus,  to  calculate  the  reactivity 
on  the  reactor  top  for  the  phase  II  design,  the  point 
source  strength  calculated  in  Section  3.3.2.3  will  be 
assumed  to  be  located  4  ft.  below  the  surface  of  the 
water  in  the  necked  core  tank.   Table  3.7.3.3.1  presents 
the  calculated  results : 

Table  3.7.3.3.1 


Material 

Thickness 

D 

B 

Water 

4  ft. 

1.47 

2.5 

Lead 

10  in. 

5.6 

10 

Totals 

7.07 

12.5 

6        2 
0(R)  (no  attenuation)  =  5.45  x  10   7's/cm  -sec 

0(R)  (attenuated)     =  5.8  7's/cm  -sec 

Conversion  factor     =  158 

Dose  Rate  =  0.037  mrem/hr 

3.7.4  Neut rons 

The  attenuation  equation  given  in  Section  3.3.4.1  can 
be  used  to  show  that  the  fast  neutron  flux  from  the  core 
is  nearly  zero  at  trie  reactor  top.   The  equation  is: 
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0  =  ^j  e      neutrons /on  -Gee 

47TIT 

R  =  13  ft. 
(H20)  =  12  ft. 

\.(H2°)  =  9*98  Cm' 
X  (Pb)  =  10  In. 

Xr(Pb)  =  8.94  era. 

S  =  2.31  x  lO1"*7  neutrons /sec 

u  2 

0  =  5.7  x  10   neutrons /cm  -sec 

The  dose  rate  due  to  fast  neutrons  on  the  reactor  top  is 
thus  only  6.6  x  lO-^  rarem/hr. 

The  dosr.  rate  due  to  thermal  neutrons  is  easily  seen 
to  be  negligible.   The  diffusion  equation  for  thermal  neu- 
trons from  an  infinite  plane  source  in  an  infinite  homo- 
geneous medium  is  : 

La"  %/L  2 

0=0     — rjj neutrons  /cm  -sec 

L  =  diffusion  length,  cm 
D  =  diffusion  coefficient,  cm 
K  =  thickness  of  attenuating  material,  cm. 
For  HpO,  L  =  2.88  cm,  D  =  0.183  cm,  and  the  column  height 

is  305  cm.   The  thermal  neutron  flux  is  thus  attenuated  by 

-45 
about  a  factor  of  10    in  the  water  column  above  the  core. 

3 .  8  Shield ing  C alci -ilations,  Reactor  Shutdown  (Phas e  II) 

The  calculations  performed  in  Section  3.4  were  repeated 

for  the  phase  II  reactor  design  in  order  to  complete  the 
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comparison  of  shielding  calculations  for  the  two  design 
phases.   For  the  case  of  the  reactor  top  shield  plug 
removed,  the  calculated  gamma  flux  at  the  coolant  tank 
surface  after  a  l4  hour  decay  period  is  7.5  x  10  gammas/ 
cm  -sec,  which  converts  to  a  dose  rate  of  22  mrem/hr. 
Were  one  fuel  element  to  be  raised  to  a  position  8  ft. 
beneath  the  surface  of  the  water  (such  as  during  a  fuel 
flipping  operation),  however,  the  dose  rate  at  the  surface 
of  the  water  would  be  about  480  mrem/hr.   Thus,  it  is 
apparent  that  personnel  should  be  protected  by  shielding 
when  conducting  refueling  or  fuel  flipping  operations. 

The  dose  rate  shortly  after  shutdown  if  coolant  is 
lost  from  the  necked  core  tank  is  calculated  to  be  12 
rem/hr.   After  l4  hours  the  dose  rate  would  decay  to  about 
700  mrem/hr. 
3 . 9  Shielding  Calculations  Summary  (Phase  II) 

3.9.1   Dose  Through  Reactor  Top  Shield  Plug 
Source  Dose  Rate  (mrem/hr) 

Core  Fission  Gammas  0.364 

Al27  (n,r)  Al28  0.485 

H1  (n,y)    H2  0.0002 

N   decay  0.037 

Fast  neutrons  0.0006 

Thermal  neutrons  negligible 


Total  O.89  mrem/hr 
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3.9.2  .  Dose  Through  Reactor  Top  Concrete 

Same  as   Phase  I  0.0018  mrem/hr 

3.9.3  Shielding  Conclusions 

The  maximum  dose  rate  anywhere  on  the  reactor  top 
is  calculated  to  he  O.89  mrem/hr.   This  figure  is  less 
than  the  design  goal  of  1.0  mrem/hr,  and  thus  the  10  in. 
lead  shield  on  top  of  the  necked  core  tank  is  considered 
to  provide  adequate  shielding  protection. 
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Chapter  4 
HEAT  TRANSFER,  PRELIMINARY  CALCULATIONS 
4. 1  General 

It  was  the  purpose  of  the  preliminary  heat  transfer 
calculations  to  determine  whether  or  not  the  heat  from 
the  proposed  High  Flux  MITR  core  could  he  adequately  re- 
moved utilizing  the  currently  installed  process  system. 
Although  the  calculations  of  this  chapter  were  originally 
done  for  the  phase  I  core  design,  they  are  equally  appli- 
cable to  the  phase  II  design.   Section  4.6  will  indicate 
the  differences  in  the  two  phases  with  respect  to  heat 
transfer. 

In  establish:!    the  limiting  temperature  parameters 
for  the  operation  of  the  High  Flux  MITR,  the  maximum 
design  values  given  in  the  technical  specifications  for 
the  present  MITR  were  adhered  to.   These  criteria  are 
that  the  temperature  at  the  hottest  spot  in  the  hottest 
fuel  plate  not  exceed  450  C  (842  F)  under  any  circum- 
stances, and  that  under  normal  operating  conditions  boil- 
ing must  not  be  allowed  to  occur  anywhere  in  the  core. 
The  450°C  temperature  limitation  was  established  by  con- 
sideration of  the  structural  stability  of  the  aluminum 
fuel  elements.   While  aluminum  does  not  melt  until  it 
readies  approximately  650°C,  it  begins  to  soften  signifi- 
cantly at  about  450°C  ( 18) .   The  requirement  that  no 
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"boiling  be  allowed  to  occur  in  the  core  was  originally 
intended  as  a  safety  factor  which  would  act  as  a  backup 
to  insure  that  the  fuel  plate  hot  spot  would  not  exceed 
450°C,  and  also  to  avoid  the  problem  of  flow  instability 
caused  by  boiling.   Neglecting  the  reactivity  effects  of 
voiding,  if  the  coolant  saturation  temperature  were  ex- 
ceeded sufficiently  to  cause  nucleate  boiling  at  some 
point  on  the  fuel  plate,  the  vapor  bubbles  formed  at  the 
fuel  plate  surface  would  either  collapse  at  the  wall  or 
would  leave  the  wall  and  collapse  after  traveling  a  short 
distance  into  the  highly  subcooled  bulk  fluid.   The  bubble 
formation  would  cause  considerable  agitation  of  the  cool- 
ant and  thereby  lead  to  a  substantial  increase  in  the 
convective  heat  transfer  coefficient.   Heat  would  thereby 
be  conducted  away  from  the  fuel  plate  surface  at  a  greater 
rate,  and  thus  tend  to  reduce  the  wall  temperature. 
However,  as  briefly  mentioned  in  the  introduction  to  this 
thesis,  the  reactivity  effect  of  voiding  in  the  area  of 
the  outermost  fuel  plates  could  be  positive,  although  all 
core  physics  calculations  done  thus  far  indicate  a  nega- 
tive value.   Because  the  possibility  of  a  local  positive 
voiding  coefficient  does  exist,  however,  one  is  provided 
with  added  incentive  to  insure  boiling  does  not  occur. 

Because  the  fuel  plates  are  so  thin  (0.050  in. ) 
and  because  aluminum  is  a  good  conductor  of  heat 
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(K  =  119  BTU/hr-ft-  F  for  pure  aluminum),  the  maximum 
temperature  at  any  point  within  a  fuel  plate  will  be  at 
most  on  the  order  of  3  to  4  degrees  F  above  the  fuel 
plate  wall  temperature  at  the  same  point. 
4. 2   Distribution  of  Reactor  Heat  Load 

The  High  Flux  MITR  will  generate  5  megawatts  of 
thermal  energy  which  must  be  removed  by  the  reactor  cool- 
ant systems.   For  the  purposes  of  fuel  plate  hot  spot 
calculations,  however,  only  that  portion  of  the  generated 
heat  which  is  trans  f erred  to  the  coolant  as  it  flows  up 
through  the  core  is  of  interest.   This  includes  heat 
transferred  to  the  D^O  coolant  both  by  conduction  through 
the  faces  of  the  fuel  plates  and  by  gamma  ray  and  neutron 
interaction  with  the  coolant. 

Wolak  has  determined  the  total  fission  energy  absorbed 
in  the  present  MITR,  operating  at  1  megawatt,  to  be  198.22 
MeV/fission  ( 1_9) .   The  specific  distribution  of  this  absor- 
bed energy  is  given  in  Table  4.2.1. 

Although  the  spectrum  from  the  proposed  high  flux 
MITR  will  be  harder  than  that  for  the  present  MITR,  and 
although  the  amount  of  aluminum  structure  may  be  somewhat 
different,  the  value  of  198.22  MeV/fission  "or  the  total 
fission  energy  absorbed  is  probably  a  good  estimate  f ■  r  the 
High  Flux  MITR.   As  shown  in  Table  4.2.1,  not  all  of  this 
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Table  4.2.1 

Fission  Energy  Absorption  Distribution  (Present 

MITR) 

Fuel  Elements MeV/fission       %   Total 

K.E.  Fission  Fragments  168.00 

K.E.  Beta  Particles  7.00 

Gamma  Energy  4.29 

179.29  90.5 

D20 

K.E.  Fast  Neutrons  5.00 

Gamma  Energy  10.24 

15.24  7.7 

Shield 

Gamma  Energy  3.69  1.8 

Total         198.22  100.0 

energy  is  absorbed  by  the  DpO  coolant  as  it  passes  through 
the  fuel  elements.   1.8$  of  the  total  fission  energy,  in 
the  form  of  gamma  ray  captures  in  the  thermal  shield,  is 
carried,  away  by  the  thermal  shield  coolant.   The  remaining 
98.2^  of  the  fission  energ   is  transferred  to  the  D?0 
coolant  by  conduction  through  the  walls  of  the  fuel  plates 
or  deposited  directly  into  the  D20.   Of  the  15.24  MeV/ 
fission  deposited  directly  into  the  Do0  within  the  reactor, 
DeVoto  has  estimated  for  the  present  MITR  operating  with  a 
19  element  core  that  9.03  MeV/fission  of  the  gamma  ray 
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energy  and  2.5  MeV/fission  of  the  energy  carried  by  the 
fast  neutrons  are  deposited  in  the  D?0  outside  the  fuel 
element  boxes  (20).   Since  the  ratio  of  D?0  outside  the 
coolant  channels  to  DpO  in  the  coolant  channels  will  be 
approximately  the  same  in  the  High  Flux  MITR  as  in  the 
present  MITR,  it  seems  feasible  to  assume  that  the  above 
figures  will  still  be  reasonably  valid.   Following 
DeVoto's  method  of  estimating  the  fraction  of  gamma 
energy  deposited  in  the  DpO  inside  the  fuel  elements  as 
the  ratio  of  D?0  volume  inside  the  fuel  elements  to  total 
DpO  volume,  one  obtains  for  the  High  Flux  MITR 

7T(l0)2(f--|f) 

,jtJ-J- —  =  0.164.   Therefore,  the  percentage  of  total 

Tr(24)d 
fissi(  posited  in  the  D?0  as  it  flows 

through  the  fuel  elements  becomes  (0.164)  (10.24)  +  1QH 
2.1$.    Due  to  the  fact  that  the  close  packed  core  will 
lead  to  greater  particle  energy  dissipation  before  the 
gamma  rays  can  escape  the  core  than  currently  exists ,    the 
above  figures  for  gamma  energy  deposited  in  the  DpO  out- 
side the  coolant  channels  may  be  a  little  high  and  there- 
fore not  conservative.   He     c,    the  difference  is  probably 
quite  small.   Hence,  it  will  be  assumed  that  92.6$  of  the 
total  fission  energy  will  be  deposited  in  the  DpO  coolant 
as  it  passes  up  through  the  fuel  elements,  and  that  90.5$ 
of  the  to1  I  fission  energy  will  be  conducted  through  the 
walls  of  the  fuel  plates  and  into  the  DpO  coolant.   (Note: 
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approximately  11  MeV  per  fission  is  carried  off  by  neu- 
trinos resulting  from  beta  decay.   However,  since  this 
neutrino  energy  is  not  deposited  anywhere  in  the  reactor 
it  is  not  considered  in  thermal  power  calculations.) 

In  the  calculations  which  follow,  the  heat  flux 
affecting  the  coolant  temperature  rise  will  be  reduced  to 
92.6$  of  the  total,  and  the  specific  heat  flux  at  the  hot 
spot  will  be  reduced  to  90.5^  of  its  computed  value  based 
on  5  megawatts . 
4.3  Hot  Channel  Factors 

For  conservative  design,  it  is  necessary  to  anticipate 
uncertainties  in  operating  conditions,  fuel  element  fabri- 
cation, design  relatio:    etc.   These  uncertainties  are 
expressed  in  terras  of  the  hot  channel  factors  which  are 
elaborated  below.   An  overall  hot  channel  factor  has  been 
estimated  for  the  present  MITR  by  DeVoto  to  be  1.25.   This 
estimate,  however,  not  only  includes  inaccuracies  intro- 
duced by  the  manufacturing  tolerances,  but  also  includes  a 
consideration  of  variations  in  flux  distribution  through 
each  element  contained  in  the  core.   DeVoto' s  value  is 
somewhat  lower  than  those  normally  found  in  the  literature 
due  to  the  fact  that  he  has  given  credit  for  several 
aspects  in  the  design  which  tend  to  reduce  the  value  of 
the  heat  flux  at  the  hot  spot.  (20)  -  It  is  because 
DeVoto 's  results  contain  consideration  of  the  neutron  flux 
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distribution,  and  the  proposed  High  Flux  MITR  differs 
radically  from  the  present  MITR  in  this  respect,  that 
direct  application  of  DeVoto's  results  cannot  be  assumed 
to  be  valid  for  the  High  Flux  MITR. 

In  determining  values  of  hot  channel  factors  to  use 
for  the  High  Flux  MITR  in  order  to  account  for  variations 
in  primary  design  variables  caused  by  manufacturing 
tolerances,  it  seems  appropriate  to  compare  the  magnitude 
of  the  tolerances  which  are  specified  for  the  present 
MITR,  and  which  will  also  apply  to  the  High  Flux  MITR, 
with  those  specified  for  the  Brookhaven  High  Flux  Beam 
Reactor  (HF3R)  (21).   In  this  n    r   the  hot  channel 
factors  determin    or  the  Brookhaven  reactor  can  be  adjus- 
ted so  as  to  mal-   I  :m  applicable  to  the  High  Flux  MITR. 
The  contribution  to  the  overall  hot  channel  factor  due  to 
the  variation  in  neutron  flux  along  an  element  will  not 
be  determined  as  a  separate  sub -factor,  but  rather  the 
variation  in  neutron  flux  will  be  used  directly  to  deter- 
mine average  power  produced  per  element  and  hot  spot 
location  on  the  element. 

Since  the  calculations  to  be  performed  in  this  chap- 
ter will  primarily  be  for  the  hottest  spot  in  the  core, 
the  hot  channel  factors  of  interest  will  be  the  heat  flux 
and  heat  transfer  coefficient  hot  channel  factors.   The 
contribution  of  each  individual  item  to  either  of  these 
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overall  hot  channel  factors  are  termed  sub-factors.   The 
sub -factors  are  multiplied  by  each  other  to  determine  the 
overall  hot  channel  factor.   A  listing  and  description  of 
the  hot  char  el  sub-factors  follows : 

1 .  Channel  Dimensional  Tolerance 

The  hot  spot  temperature  rise  due  to  the  channel 
dimensional  tolerance  will  occur  in  a  channel  of  minimum 
average  dimensions  where  the  local  dimension  is  a  maximum. 
For  the  Brookhaven  High  Flux  Beam  Reactor  (KFBR),  the 
outermost  coolant  channels  are  O.189  in.  thick  with  a 
maximum  tolerance  of  +  0.005  in.   The  sub-factor  deter- 
mined for  these  characteristics  was  1.09.   For  the  phase 
one  High  Flu?  MITRj  the  coolant  channels  are  0.117  in. 
thick  with  a  proposed  maximum  tolerance  of  +  0.005  in. 
The  sub-factor  determined  for  these  characteristics  is 

d    d    1/3  1/3 

/p  _  max,  nomN    _  .122, .  117  \      -1  11\ 

^Q  ~  d  .  ^d  .  1  ~    .112^7X12'     ~  Ltll) 

mm  mm 

This  factor  is  applied  to  the  heat  transfer  coefficient 

hot  channel  factor. 

2 .  Flow  Dl s t r ibut i on  Among  the  El cm ents 
Consideration  of  this  item  was  not  reported  in  ref- 
erence (21).   For  the  II'         ITR,  however,  the  maxi- 
mum variation  in  the  flow  distribution  among  the  elements 
is  estir.     to  be  +  10$,  as  it  is  for  the  present  MITR 
(2).   It  is  assumed  that  this  variation  wil3  increase  the 
;   ,t  transfer  coefficient  hot  channel,  facto:  by  1.10. 
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3 .   Fuel  Core  Alloy  Variation 

For  the  HFBR,  local  variations  in  the  fuel  loading 
are  limited  to  +  5$  by  the  specification  of  the  fuel 
alloy  final  thickness  at  0.020  +  0.001  in.   For  the  pre- 
sent MITR,  however,  no  specification  for  the  final  thick- 
ness of  the  fuel  alloy  section  is  made.   Rather,  it  is 
stated  that  the  overall  plate  thickness  shall  not  vary  by 
more  than  +  0.003  in.,  while  the  variation  in  the  clad 
thickness  on  either  side  of  the  core  alloy  must  be  within 
+  0.005  in.   It  is  also  stated  in  the  MIT  specifications 
that  there  is  to  be  no  objectionable  segregation  of  uran- 
ium in  the  core  alloy.   Each  plate  is  subject  to  radio- 
graphic inspection  to  determine  whether  such  segregation 
exists.   Hence,  it  is  doubtful  that  the  fuel  alloy  thick- 
ness of  the  present  MITR  plates  exceeds  the  specification 
as  quoted  for  the  HFBR.   If  the  U-235  segregation  was 
sufficient  to  exceed  this  specification,  the  fuel  plate 
would  probably  be  rejected. 

It  is  stated  in  reference  (21)  that  the  thickness 
specification  made  for  the  HFBR  can  be  met  without  great 
dj  fficulty  by  the  fuel  fabricator,  except  at  the  ends  of 
the  plates.   Here  the  fuel  alloy  develops  a  local  thick- 
ened region  due  to  Lining  process  which  is  called 
the  "dogbone".   Since  the  specj     Ion  for  the  fuel  alloy 
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thickness  (excluding  the- dogbone)  as  made  for  the  HFBR 

appears  sufficiently  easy  to  satisfy,  the  same  variation 
will  be  allowed  for  the  High  Flux  MITR.   Although  the 
dogbone  thickness  in  the  HFBR  is  required  to  be  no  more 
than  10$  greater  than  nominal  thickness,  a  variation  of 
25$  was  allowed  for  in  their  calculations  due  to  (1)  the 
fact  that  the  dogbone  thickness  is  a  sensitive  item  in 
the  hot  spot  analysis,  and  (2)  the  fact  that  this  is  a 
difficult  item  for  the  fabricator  to  control.   A  25$ 
variation  in  the  dogbone  thickness  will  also  be  assumed 
for  the  High  Flux  MITR.   This  will  give  an  added  touch  of 
conservatism  to  the  hot  spot  calculations  performed  To: 
this  design.   (In  th<  phase  II  design,  the  "dogbone"  is  a 
very  sensitive  item  because  the  fuel  plate  hot  spot  occurs 
in  the  region  of  the  dogbone  at  the  bottom  of  the  core. 
Calibrated  densitometer  measurements  will  probably  be 
required  to  determine  the  degree  of  dogboning. ) 

The  fuel  alloy  thickness  is  combined  with  the  dimen- 
sional tolerances  in  width  and  length  of  the  fuel  plates 
(and  hence  heat  transfer  area).   The  allowable  variation 
in  the  core  heat  transfer  area  for  the  HFBR  is  4.6$.   The 
same  variation  is  assumed  to  apply  to  the  High  Flux  MITR. 
Therefore,  in  both  the  HFBR  and  the  High  Flux  MITR,  the 
heat  flux  term  is  increased  by  1.10  (1.05  x  1.046) 
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everywhere,  except  at  the  ends  of  the  fuel  where  it  is 
increased  by  1.31  (1.25  x  1.046). 

4.  Reactor  Power  Measurement 

It  is  estimated  for  the  HFBR  that  the  power  calcula- 
tor will  be  capable  of  measuring  the  reactor  power  to 
within  5$.   The  same  estimate  is  made  for  the  High  Flux 
MITR.   Hence,  the  heat  flux  hot  channel  factor  is  increa- 
sed by  1.05  to  account  for  this  estimated  variation. 

5 .  Heat  Transfer  Coefficient  Deviation 

The  Phillips  Reactor  Safeguard  Committee  recommended 
the  use  of  the  modified  Colburn  correlation  in  the  HFI3R 
with  a  maximum  negative  deviation  of  2C$  (21).   This  same 
correlation  is  used  to  determine  the  heat  transfer  coef- 
ficient for  the  High  Flux  MITR  (Section  4.6.3).   Consequently, 
this  deviation  results  in  increasing  the  heat  transfer 
coefficient  term  by  1.20  for  both  the  HFBR  and  the  High 
Flux  MITR. 

6 .  Summary  of  t he  1 1  ot  Channel  Fa ctors 
Heat  Flux  Factors,  F„ 

q 

Reactor  power  measurement  error    1.05 

Fuel  alloy  variation  1.10   (1,31)* 


Total,  F    1.155  (1.375)* 


*Value  at  the  "dosboi   : 
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Heat  Transfer  Coefficient  Factors,  F, 


Channel  dimensional  tolerance  1.  ID- 

Flow  distribution  among  the  elements       1.10 
Correlation  equation  1.20 

Total,  Fh  1.46 

A  bulk  hot  channel  factor  was  not  calculated. 
Instead  the  heat  flux  hot  channel  factors  will  be  used 
when  numerically  integrating  the  power  along  a  channel 
for  the  purpose  of  calculating  bulk  temperature  rise. 
The  design  of  the  phase  one  coolant  entrance  plenum  is 
such  that  the  flow  of  coolant  up  through  the  outer  cool- 
ant channels  should  be  somewhat  greater  than  the  flow 
through  the  inner  channels.   Consequently,  it  is  argued 
that  this  effect  would  overshadow  any  bulk  hot  channel 
sub -factors  related  to  flow  velocity,  thereby  making  it 
conservative  to  simply  neglect  such  sub-factors.   It  is 
these  coolant  flow  velocity  related  sub-factors  which 
distinguish  the  bulk  hot  channel  factor  from  the  he 
flux  channel  factor.   In  addition,  by  employing  the  heat 
flux  hot  channel  factors  in  determining  the  bulk  tempera- 
ture rise,  the  c    ibution  to  the  bulk  coolant  tempera- 
ture due  to  the     tone"  effect  at  the  fuel  plate  ends 
will  be  taken  into  account.   This  effect  would  be  missed 
were  a  bulk  hoi         factor  calculated  and  used  instead. 
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4.4  Power  Distribution  in  the  Core 

In  order  to  determine  the  power  distribution  in  the 
proposed  high  flux  MITR  core,  the  modified  twenty  grand 
diffusion  code  was  used,  and  a  two  dimensional  (R,Z) 
plot  of  the  core  power  resulted.   The  basic  core  configu- 
ration and  parameters  utilized  are  described  here. 

A  cylindrical  homogenized  core  24  centimeters  in 
radius  and  65  centimeters  high  was  assumed.   The  core 
and  surrounding  reflector  area  were  divided  into  14 
regions  as  shown  in  Figure  4.4.1.   The  material  composi- 
tion of  each  region  is  listed  in  Table  4.4.1.   The  calcu- 
lations were  done  in  four  neutron  energy  groups : 

Group  1 


Group  2 
Group  3 
Group  4 


00   -  821  x  103  eV 


821  x  103   -  5.53  x  103  eV 
5.53  x  103  -  0.625  eV 
0.625      -  0.025  eV 


The  epithermal  group  cross  sections  for  Boron,  Hydrogen, 
Oxygen,  and  Uranium  (235)  were  obtained  from  the  table 
of  "Few  Group  Cross  Sections  for  Metal-H20  Systems  Fueled 
with  u235u  found  on  page  142  of  ANL-5800.   The  thermal 
cross  sections  for  the  above  three  elements  were  taken 
from  the  16  group  table  of  "Multigroup  Constants  for  Fast 
Reactor  Analysis  (LASL)"  found  on  pages  569-574  of  ANL- 
5800.   The  4  group  cross  sections  for  Deuterium  and 
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FIGURE     4.4.1. 

LAYOUT  OF  REGIONS  FOR  TWENTY 
GRAND    CALCULATION. 
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Uranium  (238)  were  obtained  by  collapsing  the  16  group 
cross  sections  to  four  energy  groups.   Fuel  plate  loading 
was  taken  to  be  1.5  times  the  fuel  loading  in  the  present 
MITR  core,  or  15  grams  of  93$  enriched  uranium  per  fuel 
plate. 

It  should  be  noted  at  this  point  that  the  purpose  of 
the  computer  calculation  was  to  obtain  a  rough  estimate 
of  the  power  distribution  in  the  core  in  order  to  have  some 
foundation  upon  which  to  base  preliminary  heat  transfer 
calculations.   The  assumptions  made  concerning  core  size, 
composition,  and  control  rod  position  could  only  be  de- 
scribed as  educated  guesses  based  on  information  available 
at  the  time.   Due  to  the  expense  involved  in  running  the 
twenty  grand  program,  only  one  run  could  be  made.   More 
detailed  and  precise  core  physics  calculations  were  con- 
currently being  performed  by  Andrews  Addae.   The  core  power 
distribution  is  shown  plotted  in  Figure  4.4.2.  As  expected, 
this  plot  shows  that  the  outermost  fuel  plate  generates  the 
greatest  amount  of  power  and  therefore  it  is  on  the  outer- 
most fuel  plate  that  the  core  hot  spot  will  be  found. 
Table  4.4.2  lists  the  power  generated  (in  watts)  at  the 
various  mesh  points  in  the  core.   (Thermal  neutron  flux  in 
both  the  core  and  reflector  has  been  plotted  in  Figure 
4.4.3  in  order  to  illustrate  why  the  core  power  distribu- 
tion takes  the  shape  it  does.) 
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FIGURE   4.4.2 
CORE   POWER    DISTRIBUTION  (watts/cm3) 
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FIGURE     4.4.3 
FLUX   MAP 


4  Sept.  68 
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Numerically  integrating  the  power  produced  by  a  fuel 
element  located  on  the  outermost  edge  of  the  core  over 
the  height  of  the  core  yields  a  total  power  of  watts  pro- 
duced in  one  of  these  "hot"  plates. 

In  order  to  determine  the  hottest  spot  on  the  hottest 
plate  in  the  High  Flux  MITR  core,  the  axial  power  variation 
along  the  hot  channel  was  investigated.   In  most  reactor 
designs  with  upward  flowing  coolant  the  hot  spot,  or  loca- 
tion of  the  maximum  plate  wall  temperature,  usually  occurs 
a  few  centimeters  above  the  point  of  maximum  power  genera- 
tion.  This  is  true,  for  instance,  for  reactors  in  which 
the  axial  flux  distribution  is  a  cosine  function  which 
peaks  at  the  geometric  center  of  the  core.   The  maximum 
wall  temperature  occurs  beyond  the  power  peak  primarily 
due  to  the  fact  that  the  coolant,  coming  in  from  the 
bottom  at  low  temperatures,  helps  keep  the  bottom  portion 
of  the  fuel  element  cooler  than  the  top.   Such  considera- 
tions can  also  be  applied  when  locating  the  hot  spot  in 
the  High  Flux  MITR,  and  based  upon  these  considerations 
we  can  expect  to  find  the  hot  spot  about  25  centimeters 
above  the  bottom  of  the  core  for  phase  one  calculations. 
4.5  Wall  Temperature  Along  Hottest  Channel 

4.5.1  General 

In  arriving  at  values  for  the  wall  temperature 
along  the  hottest  channel,  it  is  assumed  that  the  reactor 
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is  operating  at  a  state  in  which  both  the  scram  condi- 
tions for  high  coolant  outlet  temperature  and  for  low 
coolant  flow  rate  have  "been  reached  simultaneously. 
4.5.2  Formulation  of  Applicable  Equation 
To  determine  the  magnitude  of  the  wall  temperature 
at  various  points  along  the  hottest  channel,  one  may 
proceed  as  follows.   The  temperature  difference  between 
the  maximum  wall  temperature  and  the  reactor  inlet  is 
made  up  of  two  parts.   First  is  the  film  temperature 
drop,  which  is  made  up  of  the  differences  between  the 
wall  temperature  and  the  average  temperature  of  the  cool- 
ant at  that  point  in  the  channel,  Tw  -  Tc.   Second  is  the 
bulk  coolant  temperature  drop,  which  is  made  up  of  the 
difference  between  the  coolant  temperature  and  the  reactor 
inlet  temperature,  Tc  -  Ti.   Denoting  q-,/A  as  the  specific 
heat  flux  at  the  point  of  interest  and  q~  as  the  total 
rate  of  heat  transferred  to  the  coolant  from  the  channel 
entrance  up  to  the  point  of  interest,  these  temperature 

differences  may  be  written  as: 

i 

Tw  -  Tc  =  q^/hA  (4.5.2.1) 

Tc  -  Ti  =  q2/W(ch)Cp(ch)  (4.5.2.2) 

The  symbols  not  already  defined  are: 

h  =  convective  heat  transfer  coefficient 
W(ch)  =  mass  flow  rate  of  the  coolant  in  the  channel 
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Cp(ch).  =  specific  heat  evaluated  at  the  average 
bulk  temperature  of  the  D20  in  the 

channel,  i 

Adding  Equations  4.51  and  4.52,  one  obtains  an  expres- 
sion for  the  maximum  wall  temperature  in  terms  of  reactor 
inlet  temperature,  the  specific  heat  flux  at  the  point  in 
question,  and  the  total  heat  gained  by  the  coolant  from 
the  channel  entrance  up  to  the  point  of  calculation: 

Tw  =  Ti  +  q1/hA  +  q2/W(ch)Cp  +  ATt   ,     (4.5.2.3) 

(AT,   is  an  allowance  for  error  in  temperature  measure- 
ment,  taken  to  be  5°F. )   Following  the  arguments  of 
DeVoto,  it  is  more  desirable  that  wall  temperature  calcu- 
lations be  based  on  reactor  outlet  temperature,  which  is 
recognized  as  an  incompletely  mixed  bulk  temperature,  than 
on  reactor  inlet  temperature  as  above.   This  is  because 
the  outlet  temperature  recorder  is  the  first  monitor  en- 
countered by  the  coolant  after  it  leaves  the  core  and 
hence  is  the  logical  place  in  the  circuit  to  install  the 
various  alarm  and  scram  circuits. 

By  considering  the  total  heat  removed  by  the  coolant 
as  it  passes  through  the  core  of  the  High  Flux  MITR,  an 
expression  relating  the  reactor  inlet  temperature  and  the 
incompletely  mixed  reactor  outlet  temperature,  To,  may  be 
obtained.   Denoting  q  as  the  total  heat  rate  in  the 
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reactor  core,  the  expression  is: 

qr  =  Wt  Cp(r)(To"Ti)  (4.5.2.4) 

The  symbol  W.  refers  to  the  total  coolant  mass  flow  rate 

in  the  reactor  core,  and  C  /  \  is  the  specific  heat  of 

p(r) 

Ti  +  To 
the  DoO  evaluated  at  the  reactor  bulk  temperature  ~ 

Solving  this  expression  for  Ti,  and  substituting  the 

resulting  expression  into  Equation  (4.5.3)  yields  the 

result: 

ql       q2         qr 
Tw  =  To  +  jjj  +  w(ch)Cp(ch)  -  wtCp(r)  +  ATtm 

(4.5.2.5) 

4.5.3  Heat  Transfer  Coefficient 

The  heat  transfer  coefficient,  h,  for  turbulent  flow 
in  confined  rectangular  passages  is  calculated  by  means 
of  the  modified  Colburn  correlation  (see  references  (22) 
and  (23)  for  experimental  verification  of  the  use  of  the 
Colburn  equation) : 

0.8    r    ./l   0.3 


M'i  M 


!jj2£=   0.023      |iSX£|  /^UJ  (4.5.3.1) 


kf  L   **    J   f       I  iV  J    f 

where 

h  =  heat  transfer  coefficient,  BTU/hr-ft°F 
De  =  equivalent  diameter,  0.0188  ft.  (— -) 
K  =  thermal  conductivity  of  the  coolant,  BTU/hr-ft°F, 
Fig.  (4.5.3.1) 
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V  =  coolant  velocity,  11530  ft/hr  for  flow  rate  of 

1850  gpm. 
p  =  coolant  density,  lb/hr-ft,  Fig.  (4.5.3.2) 
u.  =  coolant  viscosity,  lb/ft  ,  Fig.  (4.5.3.1) 
Cp  =  specific  heat,  BTU/lb-°F,  Fig.  (4.5.3.2) 
The  subscript  "f"  means  that  the  properties  are  evaluated 
at  the  film  temperature.   The  film  temperature  is  taken 

to  be  the  average  of  the  heat  transfer  surface  temperature 

Tw  +  Tb 
and  the  bulk  water  temperature,  Tf  =  ^ .   Tw  has  been 

previously  defined,  and  Tb  is  the  bulk  temperature  of  the 

reactor  coolant  at  the  point  of  calculation.   Equation 

(4.5.1.1)  can  be  rewritten  as: 

0.8   r. 

(4.5.3.2) 


0.023K. 
h  =  - 


De 


pvl  °-8  W°- 

_HAtJ    f  {TTJf 


Implicit  in  the  above  expression  is  the  assumption  that 
the  coolant  flow  rate  is  the  same  for  all  channels,  i.e. 
W./A.  =  Wch/Ach.   Hence,  the  Reynolds  number  will  be  the 
same  whether  defined  for  one  interplate  fuel  element  or 
the  entire  core.   Heat  transfer  coefficient,  h,  is  shown 
plotted  versus  film  temperature  in  Figure  (4.5.3.3)  for  a 
constant  velocity  through  the  core  of  11530  ft/hr,  which 
corresponds  to  a  flow  rate  of  1850  gallons  per  minute  at 
a  bulk  water  temperature  of  104  F.   Taking  a  horizontal 
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out  through  the  core,  coolant  flow  area  was  taken  to  be 
the  entire  core  cross  sectional  area  times  the  ratio  of 
coolant  channel  cross  sectional  area  to  fuel  plate  cross 
sectional  area.   Taken  as  such  coolant  flow  area  is  a 
conservative  factor  because  flow  area  would  be  somewhat 
smaller  than  this  in  the  actual  core,  thereby  increasing 
flow  velocity  and  consequently  increasing  heat  transfer 
coefficient  for  the  actual  core. 
4.5.4-  Limiting  Parameters 

4.5.4.1  Reactor  Outlet  Temperature 
Following  the  reasoning  of  Section  III. 3  of  the 

Technical  Specifications  for  the  present  MITR,  the  maxi- 
mum outlet  temperature,  To,  is  taken  to  be  55°C.   (The 
technical  specifications  quote  a  maximum  outlet  tempera- 
ture of  58°C  which  is  the  audit  and  review  level.   The 
scram  point  was  set  at  55°C.) 

4.5.4.2  Maximum  Allowable  Wall  Temperature 
As  stated  in  Section  3.1,  in  order  to  insure  that 

the  hottest  spot  on  the  hottest  fuel  plate  does  not  reach 

450°C  under  any  circumstances,  the  criterion  established 

is  that  cooling  shall  be  sufficient  to  prohibit  boiling 

of  any  type  within  the  coolant  channels.   Therefore,  the 

difference  between  the  plate  wall  temperature  at  the  hot 

spot  and  the  saturation  temperature  of  the  coolant  at  the 

same  location,  Tw  ,  -  T   , ,  must  not  be  permitted  to 

max    sat  r 
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reach  the  point  of  inception  of  nucleate  boiling. 
Several  correlations  are  available  which  are  designed  to 
predict  the  amount  of  superheat  necessary  for  the  incep- 
tion of  nucleate  (24,25j  boiling.   (See  heat  transfer 
texts  by  L.S.  Tong  or  W.M.  Rohsenow.)   Three  of  these 
correlations  are  used  in  Section  4.5.3.4.   However,  such 
correlations  are  by  no  means  exact  and  so  for  the  purposes 
of  determining  the  maximum  allowable  wall  temperature  it 
is  prudent  that  we  attempt  to  design  our  system  such  that 
.  we  simply  rely  on  the  fact  that  boiling  will  not  take 
place  as  long  as  the  fuel  plate  wall  temperature  does  not 
exceed  the  coolant  saturation  temperature  at  any  point  in 
the  core.   The  top  of  the  phase  one  design  core  is  under 
8  ft.  of  D20J  at  this  pressure  the  saturation  temperature 
of  the  D20  is  223. 5°F. 

4.5.4.3  Coolant  Flow  Rate 

Following  arguments  of  Section  III. 2  of  the  techni- 
cal specifications  for  the  present  MITR,  the  minimum 
allowable  flow  rate  through  the  core  of  the  High  Flux 
MITR  will  be  taken  as  1850  gallons  per  minute. 

4.5.5  Calculations 

Utilizing  Equation  4.5.2.5*  here  repeated,  wall 
temperature  is  calculated  as  a  function  of  axial  height 
above  the  core  bottom  for  the  hottest  (outermost)  fuel 
plate. 
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Tw  "  To  +  *l/hA  +  VWp.ch  "  Wp.r  +  ATm 

(4.5.5.13 

The  first  calculation  will  be  done  in  detail  as  a  sample. 
For  a  point  on  the  outermost  fuel  plate  10  centimeters 
above  the  bottom  of  the  core : 

T  =  55°C  =  131°F,  a  limiting  condition  by  definition 


*1  =  FqGqlP(z=10)V 

F  =  heat  flux  hot  channel  factor  =  1.155 
Cr   a  fraction  of  core  power  conducted  through  the 
fuel  plates  into  the  coolant  =  0.905 
P(z=10)  =  power  produced  per  unit  core  volume  =  113.1 
-  watts/cm  (from  Table  4.5.5.1) 
V  =  volume  of  1/2  fuel  plate  plus  1/2  coolant 
channel  per  unit  axial  dimension  =  1.60  cm 
(see  Figure  4.5.5.1) 

A  =  heat  transfer  area  for  one  side  of  one  plate  per 

o 
unit  axial  length  =  7.6l  cm  (see  Figure  4.5.5.1) 

h  m   1265  BTU/hr-ft2-°F  from  Figure  4.5.3.5  for  a 
film  temperature  of  l68°F;  Fh  =  1.46  from  Sec- 
tion 4.3.6 
10 
q2  =  F  G     2  P(Z)V(Z) 
d    q  q2   z=0 


Gn  =  fraction  of  core  power  carried  off  by  the  cool- 
ant as  it  passes  through  the  core  =  O.926 


^2 


FIGURE     4.5.5.1 
CORE    VOLUME    ELEMENT 
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P(Z)  =  power  produced  per  unit  core  volume  at 

axial  position  Z  in  the  hot  channel;  values 
taken  from  Table  4.5.5.1 

V(Z)  =  volume  element  of  core  producing  power 
P(Z)  =  V  x  mesh  spacing 

w  .  =  one-half  the  coolant  mass  flow  rate  through 

ch  ° 

one  coolant  channel  =  l/2wt  ^core)^  = 
965  lb/hr 
C   .  =  average  specific  heat  in  the  hot  channel  up 
to  the  point  of  calculation  =  0.9995  from 
Figure  4.5.3.4 
q  =  total  power  produced  by  the  reactor  which  is 
removed  by  the  coolant  as  it  passes  up 
through  the  core  =  (0.926)  (5x10  )  =  4.62  x 
106  watts  =  14.5  x  106  BTU/hr 
W.  =  total  mass  flow  rate  of  coolant  through  the 
core  =  1,020,000  lb/hr  (1850  gpm  at  104°F) 
C    =  average  specific  heat  in  the  core  sr  1.0005 
for  a  core  bulk  temperature  of  122°F 
Substituting  in  the  numerical  values  and  solving: 

qx/hA  = 


(1.46)(1.155)(0.905)(3.4l2  BTU/hr -watt) (113. 1  watts/cm3) 
(1265  BTU/hr-ft2-c 

(1.60  cm3) (30.48  cm/ft)2  =  91°F 


(1265  BTU/hr-ft2-°F)(7.6l  cm2) 
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V\hCp,ch  = 

( 1. 375) (0. 926) (3.412  BTU/hr-watt) (1015.0  watts/cm2) (1J60  cm2) 

(965  lb/hr) (0.9995  BTU/lb-°F)         | 

=  7.3°F 

yty  C    =  (°'9237)(5xl06  watts) (3.412  BTU/hr-watt) 
*  t   p,t     (1.0005  BTU/lb-°F)(1.02xlOb  lb/hr) 


-  15.5°F  . 


Hence, 


T  =  131  +  91  +  7.3  -  15.5  +  5  =  218. 8°F  . 
w 

Repeating  the  above  procedure  for  the  other  values 
of  axial  distance  above  the  bottom  of  the  core  resulted  in 
the  values  plotted  on  Figure  4.5.5.2.   As  shown  on  this 
figure,  the  maximum  fuel  plate  wall  temperature  is  231°F 
and  occurs  about  23  centimeters  above  the  base  of  the 
core.  At  this  position  in  the  core,  the  saturation  tempera- 
ture of  the  D20  coolant  is  225. 4°F.   Although  the  above 
calculated  hot  spot  value  is  slightly  above  the  coolant 
saturation  temperature  at  that  point,  the  amount  of  super- 
heat is  insufficient  to  cause  boiling  by  calculation  of 
any  of  the  three  correlations  for  the  inception  of  nucleate 
boiling  given  in  Tong,  Boiling  Heat  Transfer  and  Two  Phase 
Flow  (25). 

For  the  calculated  values  of  heat  flux,  pressure,  and 
saturation  temperature  at  the  position  of  the  hot  spot,  the 
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wall  temperature  necessary  to  cause  nucleate  boiling  to 
commence  is  calculated  as  follows: 
Correlation  (1), 


2. 30.0234 
U  "  *""w  **  ^w^sat* 


q-       -  15.60  pia56(T  -T     .)- 


q"  =  124837.8  BTU/hr-ft2 

Tsat  =  225'6°F 
p  =  19.17  psi 

Tw  =  239. 1°P 
Correlation  (2), 


TLB  "  Tsat  +  ATJ+L  "  *"/h 
TLB  =  25^°F 


Correlation  (3), 

Tw  "  Tsat  =  °-015  x  lo6      (qVio6)1/4 

=  12(500  -  0.707  Tsat)10"6  lb/ft 

Tw  =  262. 0°F 

1 
Consequently,  even  if  one  were  to  abide  by  the  re- 
quirement of  Section  4.5.3.2  that  T   „   not  exceed  T   ., 

w,max  sat 

only  minor  modifications  would  be  required  to  the  fluid 
flow  system.   Since  the  main  coolant  pumps  are  capable  of 
supplying  approximately  2100  gpm  of  DpO  through  the  core 
one  could,  for  example,  simply  redefine  the  scram  point 
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flow  rate  to  be  2000  gpm.    In  addition,  one  might  re- 
define  the  outlet  temperature  scram  point  to  be  53°C.  ( 


The  combined  effect  of  redefining  scram  point  flow  rate 

and  outlet  temperature  would  be  to  reduce  T  v  by 

e  w,max  J 

approximately  8°F,  more  than  enough  to  reduce  T 
rjr  "  w,max 

below  T  +.   Preliminary  calculations,  therefore,  have 
indicated  that  fluid  flow  and  heat  transfer  would  not 
pose  a  major  design  problem. 

4.5.6  Revised  Data 

As  Andrews  Addae  began  to  obtain  reliable  results 
from  his  phase  I  core  physics  calculations,  one  fact 
became  very  clear:  the  maximum  core  power  density  would 
be  much  higher  than  the  values  used  in  the  calculations 
of  the  preceeding  section.   Best  estimates  indicated  a 
maximum  homogenized  core  power  density  of  approximately 
320  watts/cm  ,  a  factor  of  2.5  greater  than  the  maximum 
power  density  computed  in  Section  3.4.   It  should  be 
noted  that  the  computer  code  used  by  Mr.  Addae  in  his 
calculations  is  based  on  diffusion  theory,  and  as  such 
values  of  power  density  near  the  edge  of  the  core  where 
the  power  density  gradient  is  steep  are  only  approximate. 
Since  the  film  temperature  rise  term,  q./hA,  is  the 
dominant  term  in  the  calculation  of  fuel  plate  wall 
temperature  (Equation  4.5.2.5),  it  was  clear  that  to  a 
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first  approximation  one  is  faced  with  the  necessity  of 
increasing  the  magnitude  of  the  product  "hA"  by  a  factor 
of  approximately  2.5  for  the  Phase  I  design. 
4.6  Phase  I  vs  Phase  II 

With  respect  to  heat  transfer,  the  major  difference 
between  the  phase  I  and  phase  II  core  configurations  is 
that  the  phase  II  configuration  is  much  smaller.   The 
power  distribution  in  the  two  cores  is  very  similar.   For 
a  constant  mass  flow  rate,  the  flow  velocity  through  the 
phase  II  core  is  a  factor  of  2.1  greater  than  that  through 
the  phase  I  core,  which  increases  the  heat  transfer 
coefficient  by  a  factor  of  1.8l.   Partially  offsetting 
this  advantage  is  the  fact  that  the  maximum  power  density 
in  the  phase  II  core  is  a  factor  of  1.3  greater  than  that 
for  the  phase  I  core.   The  combined  effect  is  that  for  the 
light  water  cooled  core,  instead  of  the  factor  of  2.5 

augmentation  in  heat  transfer  required  by  Section  4.5.6, 

1  3 
one  must  augment  heat  transfer  by  a  factor  of  i  fli(2»5)  = 

1.8.   Maximum  power  density  for  the  phase  II  core  occurs 
along  the  bottom  edge  of  the  core. 

In  addition  to  being  smaller,  the  column  of  water 
above  the  phase  II  core  is  3  ft.  higher  than  that  above 
the  phase  I  core.   The  additional  head  of  water  increases 
the  coolant  saturation  temperature  by  about  5°F,  thus  mak- 
ing it  conservative  to  use  estimates  based  on  phase  I  cal- 
culations to  develop  a  phase  II  heat  removal  analysis. 
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Chapter  5 
AUGMENTATION  OF  FUEL  PLATE  HEAT  REMOVAL 

5.1  General 
Results  of  core  physics  calculations  being  performed  by 

Mr.  A.  Addae  (23J ,  together  with  the  results  of  the  prelimin- 
ary heat  transfer  calculations  performed  in  the  previous 
chapter  of  this  thesis,  have  indicated  that  the  product  "hA" 
must  be  increased  by  a  factor  of  approximately  2.5  (Phase  I) 
or  1.8  (Phase  II)  in  order  to  provide  adequate  heat  removal 
from  the  High  Flux  MITR  core.   It  is  the  purpose  of  the 
paragraphs  which  follow  to  explore  the  various  means  by 
which  the  product  "hA"  may  be  increased,  and  to  experimen- 
tally verify  any  conclusions  reached. 

5.2  Methods  of  Increasing  the  Product  "hA" 
5.2.1  Decreased  Fuel  Plate  Spacing 

The  fuel  plates  in  the  present  MITR  are  spaced  so  as  to 
allow  a  coolant  channel  width  of  0.117  inches.   A  survey  of 
research  reactors  utilizing  flat  plate  MITR-type  fuel  ele- 
ments has  revealed  no  parameter  which  would  prevent  plate 
spacing  from  being  decreased.   The  narrowest  flow  channel 
discovered  was  the  0.050  inch  channel  between  the  fuel 
plates  of  the  Oak  Ridge  High  Flux  Isotope  Reactor  (22). 
For  a  rectangular  channel  with  high  aspect  ratio,  decreas- 
ing the  fuel  plate  spacing  causes  a  proportional  decrease 
in  the  effective  channel  diameter  and  increase  in  flow 
velocity.   Taking  the  McAdams  correlation  for  heat  transfer 
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coefficient  as  typical  of  the  several  available, 

0.8  „n        0.4 
'e        ,Vyue. 

<  -  °-023(^ 


hD  pVD  w'u   |1C 

'  =  0.023(  — -1)    (-g&)      ,  (5.2.1.1) 


it  is  seen  that  for  a  constant  mass  flow  rate  the  decrease 
in  equivalent  channel  diameter  and  corresponding  increase 
in  flow  velocity  cancel  each  other,  thus  leaving  the 
Reynolds  number  unchanged.  As  a  result  the  Nusselt  number 
remains  unchanged,  which  in  turn  requires  that  the  magni- 
tude of  the  heat  transfer  coefficient  h  increase  in  direct 
proportion  to  the  decrease  in  flow  channel  equivalent 
diameter.   (Equivalent  diameter  is  calculated  for  rectangu- 
lar channels  using  the  standard  equation  D  = 

4  X  cross  sectional  area  x   mu^^-r^-^  **  ~~~  w^v.^  +« 

,,^4-j.^a   ~ -!w.^.i- • )  Therefore,  if  one  were  to 

wetted  perimeter      '  ' 

decrease  the  coolant  channel  width  from  the  present  0.117 
inch  to  0.090  inch,  one  could  expect  the  heat  transfer 
coefficient  to  be  increased  by  a  factor  of  1.3.   (Note:  A 
0.090"  coolant  channel  width  was  assumed  for  all  Phase  II 
calculations. ) 

5.2.2  Increased  Mass  Flow  Rate 

Increasing  the  mass  flow  rate  of  coolant  through  the 
core  results  in  a  proportional  increase  in  flow  velocity 
through  the  core.   Referring  to  Eq.  5.2.1.1  it  is  seen 
that  the  heat  transfer  coefficient  is  proportional  to 
velocity  to  the  0.8  power.   If  one  were  to  add  a  third  main 
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coolant*  pump,  thereby  increasing  the  flow  velocity  by 

approximately  50$,  one  could  expect  the  heat  transfer 
coefficient  to  be  increased  by  a  factor  of  1.38. 

When  an  increase  in  mass  flow  rate  is  contemplated, 
one  must  also  investigate  the  associated  increase  in  pres- 
sure loss  around  the  primary  coolant  circuit. 

AP  =  2f  *j  ^L  ft.  of  H20  (5.2.2.1) 

As  shown  by  Eq.  5.2.2.1,  pressure  loss  varies  with  the 
square  of  the  flow  velocity.   Consequently,  for  a  50$ 
increase  in  flow  velocity  one  could  expect  pressure  losses 
to  increase  by  a  factor  of  2.25.   In  the  present  MITR  the 
pressure  loss  through  the  core  is  on  the  order  of  0.3  feet 
of  H20  (2).   From  Eq.  5.2.2.1,  pressure  loss  through  the 
8-inch  piping  is  calculated  to  be  approximately  0.07  feet 
of  HpO  per  foot  of  pipe.   Since  each  main  coolant  pump 
has  a  capacity  for  providing  a  pressure  head  of  46  psi 
while  pumping  1000  gpm,  it  is  seen  that  neither  the  pres- 
sure loss  through  the  core  nor  the  pressure  loss  through 
the  8-inch  piping  would  prevent  a  50$  increase  in  flow 
velocity.   However,  the  pressure  loss  through  the  heat  ex- 
changers is  approximately  15  psi.   Consequently,  if  one 
were  to  add  an  additional  main  coolant  pump  to  increase 
flow  velocity,  an  additional  heat  exchanger  may  be  required 
to  limit  the  pressure  loss  across  the  heat  exchangers. 
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In  summary,  the  extensive  equipment  modifications  and 
additions  necessary  to  obtain  a  50%   increase  in  flow 
velocity  would  be  a  high  price  to  pay  for  a  factor  of 
1.38  increase  in  heat  transfer  coefficient. 

5.2.3  Modifications  to  Fuel  Plates 

The  various  methods  currently  available  to  augment 
convective  heat  transfer  have  been  summarized  in  a  paper 
by  Dr.  A.E.  Bergles  (27).   Dr.  Bergles  has  grouped  the 
augmentive  techniques  available  into  the  categories  of 
surface  promoters,  displaced  promoters,  vortex  flow, 
heated  surface  vibration,  fluid  vibration,  electrostatic 
fields,  and  additives.   Of  these  categories,  all  except 
surface  promoters  and  vortex  flow  can  be  immediately  dis- 
carded as  not  being  feasible  due  to  the  nature  of  the  MITR 
core  heat  transfer  problem.   The  idea  of  using  some  type 
of  twisted  tape  insert  in  order  to  generate  swirl  flow 
(vortex  flow)  in  the  coolant  channel  was  seriously  con- 
sidered.  However,  the  concept  was  finally  discarded  due 
to  serious  doubt  as  to  the  effectiveness  of  a  swirl 
generator  in  a  narrow  rectangular  channel  of  high  aspect 
ratio,  and  due  to  the  possible  danger  of  a  tape  breaking 
loose  and  blocking  the  coolant  flow  through  a  channel.   The 
only  applicable  method  by  which  the  product  "hA"  can  be 
increased  is,  therefore,  by  modification  of  the  fuel  plate 
surface. 
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Unfortunately,  reference  (1)   compares  the  various 
means  of  augmenting  heat  transfer  as  found  in  the  litera- 
ture with  the  corresponding  unaugmented  case  based  on  main- 
taining constant  pumping  power.   For  the  case  of  coolant 
flowing  through  the  MITR  core,  pressure  loss  through  the 
core  appears  to  be  of  minor  concern.   Consequently,  a 
comparison  of  heat  transfer  data  based  on  constant  mass 
flow  rate  is  desired.   The  results  of  a  literature  search 
for  data  reporting  the  effect  of  various  types  of  surface 
promoters  is  presented  in  the  following  section. 
5.3  Literature  Search 

5.3.1  General 

Table  5.3.1.1  was  taken  primarily  from  reference  (1) 
and  lists  the  various  experimental  investigations  of  arti- 
ficial roughness.   One  fact  stands  out  clearly;  there  is 
no  data  for  water  cooled  thin  rectangular  plates.   Results 
of  investigations  which  most  nearly  approximate  the  situ- 
ation in  the  MITR  core  are  presented  below.   No  attempt  at 
completeness  is  made  here.   The  original  papers  should  be 
consulted  if  detailed  information  is  desired. 

5.3.2  Brouillette,  Mifflin,  and  Myers 

The  authors'  purpose  was  to  determine  the  effect  of 
a  type  of  internal  fin  on  the  heat  transmission  and  pres- 
sure drop  characteristics  of  water  flowing  through  copper 
tubing.   The  fins  were  produced  by  cutting  60-degree, 
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Table  5.3.1.1 
EXPERIMENTAL  INVESTIGATIONS  OF  ARTIFICIAL  ROUGHNESS  - 
NONBOILING  HEAT  TRANSFER  AND  FRICTION 


Investigators 
Cope 

Sams 

Kolar 

Lancet 

Dipprey  and  Saber sky 

Nunner 

Koch 

Nagaoka  and  Watanabe 

Saras 

Kreith  and  Margolis 

i 

Boelter  et  al. 

■ 

Edwards  and  Sheriff 
Kemeny  and  Cyphers 


Description 

Cooling  of  water;  internally 
knurled  tubes 

Heating  of  air  at  high  AT;  square 
threads  inside  tubes 

Heating  of  water  and  air;  internally 
threaded  tubes 

Heating  of  air;  rectangular  channel 
with  machined  roughness 

Heating  of  water;  sand-grain  type 
roughness  in  tubes 

Heating  of  air;  rings  of  various 
cross  section  inserted  inside  tubes 

Heating  of  air;  axially  supported 
rings  inserted  inside  tube 

Heating  of  water;  wire  coils  inside 
tubes 

Heating  of  air;  wire  coils  inside 
tubes 

Heating  of  water  and  air;  tubes 
with  wire-coil  inserts 

Heating  of  air;  rectangular  duct, 
one  surface  heated,  with  transverse 
strips 

Heating  of  air;  rectangular  channel, 
one  surface  heated  lined  with  wire 

Heating  of  water;  annulus,  inner 
heated  surface  with  spiral  grooves 
or  protrusions 
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Table  5.3.1.1 
(Continued) 


Investigators 
Brauer 

Bennett  and  Kearsey 

Brouillette,  Mifflin, 
and  Myers 

Durant  et  al. 
Wilkie 

Sheriff  and  Gumley 
Kjellstrom  and  Larsson 
Draycott  and  Lawther 

Kattchee  and  Mackewicz 


Description 


Heating  of  water;  annulus;  inner 
heated  surface  with  grooves 

Heating  of  superheated  steam;  annu- 
lus, inner  tube  heated  with  grooves 

Heating  of  water;  tubes  with  inter- 
nal V-shaped  fins 

Heating  of  water;  annulus,  inner 
heated  tube  with  diamond  knurls 

Heating  of  air;  annulus,  inner  tube 
with  square  or  rectangular  protru- 
sions 

Heating  of  air;  annulus,  inner  heated 
tube  with  coiled  wires 

Heating  of  air;  annulus,  inner  heated 
tube  with  various  protrusions 

Heating  of  air;  annuli  and  tube 
bundles  with  threads,  knurls,  and 
coiled  wires 

Heating  of  nitrogen  (water  for  fric- 
tion tests);  wire  coils  outside  rods 
in  a  bundle 


Sutherland 


Heating  of  air;  wire  coils  outside 
rods  in  a  bundle 
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Nunner  correlated  his  data  for  Reynolds  numbers  ranging 
from  10  to  Kr  with  the  result: 


f/8  RePr     .  ft-   o  o  A 

Nu  =  — ^— ^ s ;         (5.3.3.1) 

1  +  —  £.•?     ..  (pr  J-  -  1) 
Re1/0Pr1/bV   *o    ' 

where  f  is  the  Darcy  friction  factor  (fn  =  2  — -  .  — «■) 
for  the  rough  surface  and  f  is  the  smooth  tube  friction 
factor.   Equivalent  diameter  for  all  rough  tube  calcula- 
tions is  taken  to  be  the  volumetric  mean  diameter,  d  = 

m 

Nunner' s  data  indicated  that  the  larger  the  6/d 
ratio  the  greater  was  the  increase  in  heat  transfer  coeffici- 
ent.  His  data  also  indicated  that  the  optimum  value  of 
the  ratio  € /L,  where  6  is  the  height  of  the  ring  insert  and 
L  is  the  center-to-center  spacing  of  the  insert  in  the  tube, 
is  approximately  10.   However,  since  the  manually  installed 
ring  inserts  in  Nunner' s  experiment  may  not  have  been  per- 
fectly seated  against  the  tube  wall,  the  effect  of  increas- 
ing the  number  of  ring  inserts  per  inch  of  tube  length  may 
have  been  to  decrease  (or  at  least  fail  to  increase)  the 
effective  area  available  for  heat  transfer.   Consequently, 
it  is  not  possible  to  predict  the  optimum  fin  spacing  dis- 
tance based  on  the  results  of  Nunner 's  work.   Nunner 's 
work  did  indicate,  however,  that  ring  inserts  were  capable 
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V-shaped  notches  on  the. inner  surface  of  five  specimens  of 
tubing.   The  number  of  notches  or  fins  was  varied  from  8 
to  32  per  inch  of  length  and  the  ratio  of  fin  height  to 
tube  diameter  ranged  from  0.009  to  0.050  inch.   Data  was 
taken  for  Reynolds  numbers  ranging  from  1.5  x  10  to 
2  x  10  .  Table  5.3.2.1  below  summarizes  the  results  obtained. 
In  addition,  the  authors  stated  that  fin  height  has  a 
greater  effect  in  increasing  the  heat  transfer  coefficient 
than  the  number  of  fins. 

Table  5.3.2.1 


Tube 
No. 

€ 

d 

hAi 
hAi 

ffin) 

[smooth) 

A(fin) 
A (smooth) 

tiffin) 

h( smooth) 

1 

0.009 

1.18 

1.33 

0.89 

2 

0.013 

1.08 

1.18 

0.91 

3 

0.017 

1.55 

1.48 

1.04 

4 

0.027 

1.90 

1.44 

1.32 

5 

0.050 

2.48 

1.48 

1.68 

5.3.3     Nunner 

Nunner  investigated  the  effect  of  various  sizes, 
shapes,  and  densities  of  ring  inserts  on  heat  transfer 
coefficient  and  friction  factor  for  air  flowing  through  a 
50  mm  diameter  brass  tube.   Both  round  and  rectangular 
springy  split  rings  were  manually  inserted  into  the  tube 
at  equal  distances  apart  to  serve  as  roughness  elements. 
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of  increasing  the  Nusselt  number  for  air  flowing  in  a 
cylindrical  duct  by  nearly  a  factor  of  three. 

5.3.4  Kolar 

Kolar  analyzed  the  mechanism  of  heat  transfer  for 
both  air  and  water  flowing  through  smooth  and  rough  tubes. 
Internal  tube  roughness  was  formed  by  machining  a  60° 
triangular  thread  into  the  tube.   The  smooth  tube  was  26 
mm  in  diameter  and  the  rough  tube  thread  depth  varied  from 
0.5  to  1.5  nim.  Kolar 's  smooth  tube  data  fit  the  standard 
Dittus-Boelter  correlation; 

Nu  =  0.023  Re0,8Pr0*^.  (5.3.4.1) 

In  order  to  correlate  his  rough  tube  data  Kolar  proposed 
a  relation  in  the  same  form  as  the  Dittus-Boelter  equation 
for  smooth  surfaces : 

Nu  =  0.0517  RefPr0,5  ,  (5.3.4.2) 

where  Ref  is  the  Reynolds  number  defined  by  the  friction 
velocity, 

u„  =  u  \/T78  , ' 

f    m  v  '      ' 

wherein  u  is  the  mean  velocity  and  u^  is  the  friction 
m  1 

velocity.  Kolar  shows  that  his  data  is  in  general  agree- 
ment with  that  of  Nunner  and  also  with  that  of  Cope  (30). 
(Cope's  data  was  for  water  flowing  in  circular  tubes  with 
the  inside  surface  roughened  by  knurling)   Reynolds  number 
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was  varied  from  4500  to  345,000  and  Prandtl  number  from 
0.71  to  5.52.  Unfortunately,  Kolar  makes  no  mention  of 
the  fin  effect  and  what  influence  this  effect  might  have 
on  the  rate  of  heat  transfer. 

5.3.5  Lancet 

Lancet  investigated  the  effect  of  a  knurled  surface  on 
the  heat  transfer  coefficient  for  air  flowing  through  a 
small  rectangular  duct.   The  rough  tube  was  4.5  inches  long 
and  the  maximum  width  and  height  were  0.145  in.  and  0.045 
in. ,  respectively.   The  diamond-shaped  protrusions  were 
approximately  0.010  in.  x  0.010  in.  x  0.010  in.  and  were 
spaced  approximately  0.010  in.  apart  leaving  an  open  duct 
area  of  0.025  in.  x  0.125  in.   For  a  constant  Reynolds 
number  in  the  range  3110 <  Re <  26900,  Lancet  found  that  the 
Nusselt  number  was  increased  by  about  a  factor  of  two. 
Lancet's  data  is  fit  quite  well  by  Nunner's  relation  for 
Nusselt  number,  Eq.  2.3.1,  but  even  better  by  Nu  =  0.042 
Re  '   Pr  '    .   The  fin  effect  is  not  mentioned  in  Lancet's 
results. 

Lancet's  results  for  a  smooth  rectangular  channel  were 
found  to  fit  the  Colburn  equation  for  smooth  round  ducts : 

Nu  =  0.023  Re°*8Pr1/3  .  (5.3.5.1) 

5.3.6  Conclusions 

From  the  literature  survey,  it  is  apparent  that  the 
increase  in  heat  transfer  from  a  surface  roughened  duct  is 
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proportional  to  the  increase  in  the  duct  roughness  ratio, 

€/d.   When  correlating  cylindrical  tube  data,  the  equiva- 

j 
lent  diameter,  d,  is  taken  to  be  a  volumetric  mean  diameter 

defined  as  d  =  J^-   .   When  attempting  to  apply  the  results 

of  experiments  performed  on  cylindrical  tubes  in  order  to 

predict  the  effect  of  roughness  on  rectangular  ducts  of 

high  aspect  ratio,  it  is  not  clear  how  equivalent  diameter 

should  be  defined.   (Aspect  ratio  is  the  ratio  of  width  to 

height  of  a  rectangular  channel.)  When  applying  smooth 

tube  correlations  such  as  those  of  McAdams  or  Colburn  to 

rectangular  ducts  equivalent  diameter  is  taken  to  be  d  = 

4  x  cross  sectional  area  ,4AN   c      .        .  „  .   _ 
wetter  perimeter   M  •   Synan,  however,  has  shown 

that  the  rough  rectangular  channel  data  of  Lancet  is  best 
correlated  with  the  data  of  other  investigators  by  using 
the  volumetric  mean  diameter,  d  =  J—r»   as  the  equivalent 
diameter  (35).   It  should  be  noted,  though,  that  the 
aspect  ratio  for  Lancet's  rectangular  ducts  was  only  5  :  1. 
For  a  MITR  coolant  channel,  the  aspect  ratio  is  approxi- 
mately 28  :  1,  a  considerable  difference. 

For  the  proposed  MITR  coolant  channel  shown  in  Figure 
5.4.2.1, 

*l  -  Tr  "  °'174  in*   ;  I  =  °-115 

d2  =  It  "  °'496  in'   ^   |  =  0.04 
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Since  i»t  has  been  concluded  that  the  increase  in  heat  trans- 
fer is  proportional  to  the  roughness  ratio,   £/d,  utiliza- 
tion of  the  correct  value  of  d  is  imperative  if  one  is  to 
predict  rectangular  duct  heat  transfer  coefficients  based 
on  cylindrical  tube  experiments.   From  reference  (35)  one 
would  be  tempted  to  make  the  rectangular-cylindrical 
comparison  based  on  d  =  v/W/ttL  .   However,  using  the 
volumetric  equivalent  diameter  implies  that  the  increase  in 
heat  transfer  due  to  roughness  will  vary  with  the  width  of 
the  rectangular  channel.   For  example,  for  a  given  roughness, 

£ ,  a  3-inch  wide  thin  rectangular  channel  would  have  a 
smaller  €/d  ratio  than  a  2  1/2  inch  wide  rectangular  chan- 
nel, and  thus  a  lower  Nusselt  number.   This  is  obviously 
not  correct.   It  would  appear,  therefore,  that  for  an  arti- 
ficially roughened  thin  rectangular  channel  with  high  aspect 
ratio,  the  correct  equivalent  diameter  lies  somewhere 
between  4a/P  and  J-j-   .   Using  the  pessimistic,  J-j-  t 
equivalent  diameter,  the  heat  transfer  coefficient  for  the 
channel  of  Figure  5.4.2.1  can  be  estimated  (using  the  same 
velocity  and  temperature  conditions  as  in  Section  4.5.3): 

Area 
.  .        .      Increase   Ix, 

d     £  nDB  nK       nN       Due 


(in.)   "3    (Dittus-Boelter)   (Kolar)   (Nunner)   to  Fins   DB 
0.496  0.04       1360  2520     2270     I.67    1.83 
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On  the  basis  of  these  estimated  values  it  can  be  concluded 

0 

that  the  finned  channel  of  Figure  5.4.2.1  will  provide  at 
least  a  factor  of  1.5  increase  in  heat  transfer  over  an 
equivalent  smooth  channel.   The  heat  transfer  augmentation 
may,  however,  be  greater  than  a  factor  of  1.5*  possibly  as 
much  as  a  factor  of  3  or  4. 
5.4  Experimental  Program 

5.4.1  General 

Based  on  the  literature  survey  of  the  previous  section, 
it  was  apparent  that  heat  transfer  from  artificially 
roughened  cylindrical  tubes  is  at  least  partially  understood 
by  correlation  of  roughness  ratios,  €/d.   However,  it  was 
not  clear  how  one  should  apply  the  results  of  rough  tube 
experiments  to  predict  rough  rectangular  channel  heat  trans- 
fer.  An  experiment  which  would  duplicate  as  nearly  as  possi- 
ble the  conditions  which  will  exist  in  the  High  Flux  MITR 
core  was  clearly  desirable. 

5.4.2  Coolant  Channels 

After  making  a  brief  survey  of  the  methods  available 
to  heat  the  walls  of  a  test  section  coolant  channel,  it 
was  decided  to  heat  the  coolant  channel  by  passing  an  elec- 
trical current  through  the  channel  walls.   As  described  in 
Section  5.4.3.2  the  heat  transfer  laboratory  power  supply 
is  capable  of  supplying  3000  amperes  at  24  volts. 
Consequently,  if  the  test  section  is  to  take  full  advantage 
of  the  power  supply  available,  it  should  provide  a  resistance 
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of  8  x  -10   ohms.   Since, 

R  =  -^  (5.4.2.1) 

where  R  =  resistance  (ohms) 

L  =  channel  length  (cm),  and 

o 
A  =  cross  sectional  area  for  current  flow  (cm  ), 

p  =  resistivity  (ohm-cm) 

there  are  three  parameters  which  can  be  adjusted  in  order 

to  obtain  a  correct  resistance  value.  A  stainless  steel 

channel  wall  2.5  inches  wide  by  0.020  inches  thick  was 

decided  upon.   For  stainless  steel  of  these  dimensions,  a 

coolant  channel  wall  electrical  resistance  of  6.8  x  10 

ohms  was  calculated.   Maximum  power  input  into  the  coolant 

channel  was  therefore  estimated  to  be  6l,200  watts,  which 

will  provide  a  maximum  heat  flux  into  the  coolant  flowing 

through  the  channel  of  2.5  x  105  — ^^U  # 

hr-ft 
In  order  to  provide  a  basis  for  future  comparison, 

the  first  channel  to  be  constructed  and  tested  was  a 

smooth  surface  rectangular  duct.   A  flat  piece  of  0.020 

inch  thick  stainless  steel  was  folded  three  times  to  form 

a  rectangular  channel  with  inside  dimensions  of  2.5"  x 

0.090".   The  channel  was  closed  by  a  single  weld  which 

ran  the  entire  length  of  the  coolant  channel.   It  was 

expected  that  test  results  obtained  from  this  channel  would 

agree  fairly  well  with  existing  smooth  tube  correlations, 
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such  as*  those  of  Colburn  or  McAdams. 

The  second  channel  constructed  was  designed  to  provide 
a  "proof  of  principle"  for  internally  roughened  rectangular 
ducts.   It  was  decided  to  mill  fins  into  the  stainless 
steel  channel  wall  in  preference  to  other  means  of  wall 
roughening  (such  as  knurling)  for  the  rough  wall  run 
because,  (1)  the  increase  in  heat  transfer  due  to  heat 
transfer  surface  area  versus  that  due  to  increased  turbu- 
lence could  be  accurately  separated,  and  (2)  the  process 
was  one  which  was  within  the  immediate  capability  of  shop 
equipment  and  personnel. 

Fin  height  was  chosen  to  be  0.020  inches  based  on  the 
philosophy  that  in  the  High  Flux  MITR  no  point  of  two  fuel 
plates  should  be  closer  than  the  smallest  existing  separa- 
tion of  reactor  fuel  plates,  the  0.050  inch  separation  of 
the  Oak  Ridge  High  Flux  Isotope  Reactor  fuel  plates.   The 
projection  of  opposing  0.020  inch  fins  into  a  coolant 
channel  whose  unfinned  width  was  0.090  inches  results  in 
a  minimum  spacing  of  0.050  inches.   Fin  width  was  taken 
to  be  0.020  inches  in  order  to  provide  the  fins  with 
structural  stability  and  to  prevent  the  machining  problems 
which  would  accompany  the  milling  of  very  narrow  fins. 
Fin  spacing  was  a  compromise,  and  in  fact,  somewhat  arbi- 
trary.  From  Nunner's  work  (Section  5.3.3)  the  optimum  fin 
spacing  to  fin  height  ratio  appears  to  be  10  :  1,  which 
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would  require  a  fin  spacing  of  0.200  inches  for  a  0.020 
inch  fin.   However,  from  the  work  of  Brouillette,  Mifflin, 
and  Myers  it  appears  as  though  the  optimum  arrangement  is 
to  have  fins  spaced  as  close  together  as  possible.   It  is 
this  author's  opinion  that  the  manually  inserted  rings  of 
Nunner's  experiment  are  not  capable  of  providing  data  on 
which  to  base  conclusions  concerning  optimum  fin  spacing. 
It  is  felt  that  heat  transfer  will  be  increased  by  spacing 
fins  close  together  due  to  the  greater  heat  transfer  area 
increase  over  widely  spaced  fins.   However,  there  must 
exist  some  minimum  spacing  beyond  which  fluid  circulation 
between  fins  will  be  hindered  and  boiling  in  the  inter- 
fin  spaces  will  be  caused  to  occur.   With  the  above  infor- 
mation in  mind,  a  fin  spacing  of  0.040  inches  was  chosen. 
The  optimization  of  fin  spacing  might  well  be  the  object 
of  future  experiments. 

Rectangular  grooves  0.040  inches  wide,  0.020  inches 
deep,  and  spaced  0.020  inches  apart  were  milled  cross-wise 
three  at  a  time  into  a  flat  piece  of  stainless  steel.   Two 
grooves  0.100  inch  wide  were  then  milled  lengthwise  into 
the  stainless  steel  plate  at  positions  which  were  to 
become  the  channel  edges.   The  stainless  steel  plate  was 
next  folded  into  a  rectangular  channel  and  welded  as  with 
the  smooth  channel. 
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5.4.3  Description  of  Apparatus 

5.4.3.1  Hydraulic  System 

The  experimental  facility  used  was  the  low-pressure 

j 
test  loop  located  in  the  MIT  heat  transfer  laboratory. 

The  basic  apparatus  was  designed  and  constructed  in  1961. 
A  schematic  of  the  loop  is  presented  in  Figure  5.4.3.1.1. 
The  pipes  and  fittings  are  made  of  brass  and  stainless 
steel  for  corrosion  resistance.   Rayon  reinforced  rubber 
hose  is  used  where  flexible  connections  are  required. 
Distilled  water  is  circulated  by  a  bronze,  two-stage, 
regenerative  pump  capable  of  providing  a  discharge  pres- 
sure of  260  psig  at  3.6  gpm.   The  pump  is  driven  through  a 
flexible  coupling  by  a  3-hp  Allis -Chalmers  induction  motor. 
A  Fulflo  filter  is  installed  at  the  pump  inlet.   Pressure 
fluctuations  at  the  outlet  of  the  pump  are  damped  out  by 
means  of  a  2.5  gallon  Greer  accumulator  charged  with 
nitrogen.   The  accumulator  contains  a  flexible  bladder- 
type  separator  which  prevents  the  nitrogen  from  being  absor- 
bed by  the  system  water.   After  the  accumulator,  the  flow 
splits  into  a  bypass  line  and  a  test-section  line. 

In  the  test-section  line,  fluid  flows  through  a 
Fischer-Porter  flowrator  followed  by  a  preheater,  thence 
through  a  Hoke  metering  valve  and  the  test  section,  after 
which  it  merges  with  fluid  from  the  bypass  line.   The  flow 
then  goes  through  the  heat  exchanger  and  returns  to  the 
pump.   The  preheater  consists  of  four  Chromalox  heaters  of 


182 


uj  "- 

_J 
_j  or 

v>  < 
5£ 


4 


D 


< 
O 


OL 


^J- 


% 


SU3iV3H| 


V 


XNV1    9NISSV930 
0± 


I 


O 

Is 

ui  a; 

z  a. 
uj  5 
o  55 


syoivaMoii 


^<E) 


o 


5 
O 

u 

< 


o 

< 
u. 


z 

|H    UJ 

^5 

r^  UJ 

•  a. 

0    UI 
W   O 


< 
5 

UJ 

X 

o 
to 


183 

approximately  6  Kw  each.   Three  of  the  heaters  are  con- 
trolled simply  with  "on-off"  switches  while  the  fourth 
can  provide  a  continuous  range  of  0  to  6  Kw  by  means  of 
a  bank  of  two  variacs  mounted  on  the  test  bench.   Quick- 
action  Jamesbury  ball  valves  are  installed  before  the  in- 
let to  the  flowrator  and  after  the  exit  from  the  test 
section.   The  ball  valves  permit  quick  isolation  of  the 
test  section  to  minimize  fluid  loss  when  conducting  burn- 
out tests.   The  exit  valve  can  also  be  used  to  adjust  the 
test  section  pressure. 

Flow  through  the  bypass  line  is  controlled  by  a  ball 
valve  on  each  side  of  which  there  is  a  300  psig  pressure 
gage.   Pump  operating  pressure,  and  hence  the  pressure  up- 
stream of  the  test  section,  is  controlled  by  this  valve. 

The  heat  exchanger  is  a  counterflow  type  with  system 
water  flowing  in  the  inner  tube  and  city  water  in  the 
outer  annulus.   The  heat  exchanger  is  capable  of  maintain- 
ing a  constant  pump  inlet  temperature.   A  Fulflo  filter 
is  installed  on  the  city  water  line  to  reduce  scale  forma- 
tion in  the  exchanger. 

The  distilled  water  was  de-ionized  continuously  by 
passing  a  portion  of  the  flow  through  four  mixed-bed 
resin  demineralizer  cartridges  installed  in  a  Barnstead- 
"Bantam"  demineralixer  unit.   A  4.7  gallon  degassing  tank 
is  provided  with  five  electrical  heaters  (3-200  Vac  and 
2-110  Vac).   The  degassing  tank  also  serves  as  a  surge  tank. 
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A  15-gallon  stainless  steel  storage  tank  for  filling  the 
system  is  mounted  directly  above  the  degassing  tank  and 
can  be  filled  with  distilled  water  from  standard  5-gallon 

bottles  with  a  small  Hypro  pump. 

j 

5.4.3.2  Power  Supply 

Test  section  electrical  power  is  provided  by  two 
36-KW  dc  generators  connected  in  series.   Each  generator 
is  rated  at  12  volts  and  3,000  amperes.   The  power  control 
console  permits  coarse  or  fine  control  from  0  to  24  volts. 
A  water-cooled  shunt  installed  in  parallel  with  the  test 
section  protects  the  generators  against  the  shock  of  the 
sudden  open  circuit  which  would  occur  at  burnout.   Power 
is  transmitted  from  the  main  bus  to  the  test  section  by 
water-cooled  power  leads.   In  the  present  tests,  the  down- 
stream or  exit  power  lead  was  at  ground  potential.   Rubber 
hose  connected  both  inlet  and  exit  chamber  plenums  to  the 
main  loop  in  order  to  electrically  isolate  the  test  sec- 
tion. 

5.4.3.3  Instrumentation 

All  temperatures  were  measured  by  copper -constant an 
thermocouples  made  from  30-gage  Leeds  and  Northrup  duplex 
wire.   The  test  section  inlet  and  outlet  temperatures  were 
measured  by  thermocouples  directly  in  the  fluid  stream. 
The  thermocouples  were  introduced  into  the  piping 
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immediately  before  and  after  the  test  section  through 
Conax  fittings. 

Coolant  channel  inlet  and  exit  pressures  were  measured 
with  Helicoid  8-1/2  inch  gages  with  full  scale  readings  of 
200  psig  and  100  psig,  respectively.   Both  are  specified 
to  an  accuracy  of  +  0.25$  of  full  scale.   Pressure  loss 
across  the  test  section  was  also  measured  by  a  manometer 
located  adjacent  to  the  test  stand. 

A  float  tube  calibrated  for  flows  of  0  to  8  gallons 
per  minute  was  installed  in  the  basic  Fischer-Porter 
flowmeter  housing  and  provided  measurement  of  the  test 
section  flow. 

The  voltage  drop  across  the  test  section  was  read 
directly  on  a  Digatec  multiple  range  dc  voltmeter  with  a 
specified  accuracy  of  +  1/2  percent.   The  current  flow  was 
determined  by  using  a  Wilson-Maeulen  multi-volt  meter 
(having  0-30  and  0-60  millivolt  scales)  to  measure  the 
voltage  drop  across  a  calibrated  shunt  (60.17  amp/mV)  in 
series  with  the  test  section.   Thermocouple  voltages  were 
read  from  a  Honeywell -Brown  continuous  reading  chart 
recorder. 

5.4.3.4  Test  Section 

The  test  section  consisted  of  a  stainless  steel  cool- 
ant channel  (described  in  Section  5.4.2)  which  protruded 
into  two  stainless  steel  end  blocks.  The  end  blocks 
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served, to  act  as  an  entrance  plenum  for  the  coolant  water, 
to  provide  a  point  of  attachment  for  the  electrodes  of  the 
power  supply,  and  to  serve  as  a  point  of  attachment  for 
the  pressure  taps.   Each  end  block  was  fabricated  in  three 
pieces j  a  base  piece  into  which  a  1/2  inch  NPT  was  machined 
in  order  to  provide  a  point  of  attachment  for  the  coolant 
hose,  and  two  entrance  plenum  sections  which  were  machined 
so  as  to  taper  the  water  channel  from  a  1/2  inch  diameter 
cylinder  to  a  0.130  inch  by  2.5  inch  rectangle.   Figure 
5.4-. 3.4.1  shows  both  a  horizontal  and  a  vertical  cross  sec- 
tion through  an  end  block.   Figures  5.4.3.4.2  and  5.4.3.4.3 
show  the  test  section  in  various  stages  of  construction  and 
assembly.   The  rectangular  coolant  channels  protruded  into 
a  tight-fitting  slot  in  the  end  blocks  to  a  distance  of 
3/4  inch.   Once  inserted,  the  coolant  channel  and  end 
blocks  were  bonded  together  with  a  stainless  steel  solder, 
and  pressure  taps  were  drilled  into  the  coolant  channel  at 
the  positions  indicated  in  Figure  5.4.3.4.1.   During  the 
finned  channel  experimental  runs,  it  was  discovered  that 
water  was  seeping  around  behind  the  coolant  channel  inside 
the  end  blocks.   Consequently,  the  pressure  taps,  which 
protruded  into  the  end  blocks,  were  reading  the  pressures 
of  the  inlet  and  exit  plenums  rather  than  the  pressures 
inside  the  inlet  and  exit  of  the  coolant  channel  itself, 
as  intended.   In  order  to  eliminate  the  inlet  and  exit 
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A.  Three-piece  end  block  "before 
assembly.   Coolant  channel  in 
foreground. 


B.  Channel  end  block  showing 
coolant  channel  slot. 


C.  Finned  coolant  channel  inter- 
ior shown  before  bending  and 
welding  plate  into  rectangular 
channel. 


D.  Assembled  coolant  channel 
and  end  blocks. 


FIGURE  5.^.3.4.2 
EXPERIMENTAL  TEST  SECTION 
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A.  Coolant  channel  inserted  into 
support  block. 


B.   Test  section  assembled. 
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C.  Test  section  in  place  in  test    D.   Instrumentation  stand, 
facility. 


FIGURE  5.4.3.4.3 
EXPERIMENTAL  TEST  SECTION 
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effects  and  thus  to  measure  only  the  pressure  drop  across 
the  coolant  channel,  pressure  sensing  tubes  were  inserted 
through  the  endmost  thermocouple  holes  in  the  front 
fiberglass  block  and  through  the  coolant  channel  wall  at 
those  points.   The  existing  pressure  tap  holes  in  the  end 
blocks  were  plugged. 

The  coolant  channel  was  physically  supported  by,  and 
was  electrically  and  thermally  insulated  by  two  2"  x  6" 
Nema  grade  G-5  fiberglass  blocks.  A  rectangular  channel 
2  1/2"  x  0.125"  was  milled  into  one  of  the  fiberglass 
blocks  in  order  to  accept  the  coolant  channel.   The  two 
fiberglass  blocks  were  held  together  by  16  1/8"  bolts 
spaced  along  either  side  of  the  coolant  channel.   Holes 
were  drilled  in  the  fiberglass  blocks  as  shown  in  Figure 
5.4.3.4.1  in  order  to  provide  access  for  the  thermocouples 
to  the  coolant  channel. 

Teflon  tape  was  used  to  electrically  insulate  the 
coolant  channel  from  the  thermocouples.  A  total  of  seven- 
teen thermocouples  were  mounted  on  the  coolant  channel 
(see  Figure  5.4.3.4.4  for  arrangement).   The  thermocouples 
were  pressed  into  place  and  held  onto  the  coolant  channel 
by  the  fiberglass  support  blocks.   Four  thermocouples  were 
mounted  near  the  surface  of  the  fiberglass  blocks  in  order 
to  monitor  the  temperature  gradient  across  these  insulating/ 
support  blocks.   In  order  to  maintain  a  near  zero  temperature 


FIGURE  5.4.3A4 
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THERMOCOUPLE     ARRANGEMENT  ON 
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gradient,  three  separate  flexible  tape  heaters  were  wrapped 
around  the  fiberglass  blocks.   Heat  output  from  these 
three  heaters  was  adjusted  using  three  variacs  to  maintain 
a  nearly  zero  temperature  difference  between  the  outside 
channel  wall  and  the  fiberglass  block  surface,  thus  creating 
an  adiabatic  wall  condition  at  the  channel  outer  wall.   The 
fiberglass  blocks  and  heater  strips  were  blanketed  with 
fiberglass  insulation  in  order  to  reduce  the  required 
heater  output  and  to  cause  the  strip  heaters  to  give  a  more 
uniform  heating  effect. 

Four  steel  support  arms  were  bolted  between  the  fiber- 
glass blocks  and  the  stainless  steel  end  blocks  in  order 
to  join  the  entire  test  section  more  firmly.   By  so  doing, 
the  probability  that  the  solder  joints  between  the  coolant 
channel  and  end  blocks  would  crack  during  handling  was 
reduced. 

The  method  used  to  change  coolant  channels  in  the 
test  section  was  to  heat  the  end  blocks  (not  the  coolant 
channel)  sufficiently  to  melt  the  solder  permitting  removal 
of  the  channel.   Once  removed  and  allowed  to  cool,  the  end 
blocks  were  emissed  in  a  crock  of  nitric  acid.   The  nitric 
acid  successfully  removed  all  the  remaining  solder  from 
the  end  blocks  without  any  damage  to  the  stainless  steel, 
thus  providing  a  clean  smooth  surface  for  installation  of 
the  new  channel. 
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5. "4. 4  Operating  Procedure 

With  all  system  components  installed  and  ready  for 
operation,  and  with  both  the  test  section  inlet  and  exit 
valves  closed,  the  coolant  pump  was  started  and  the  bypass 
valve  adjusted  to  provide  20  psig  backpressure  to  the  test 
section.   Next  the  test  section  exit  valve  followed  by  the 
test  section  inlet  valves  were  opened  and  adjusted  to  pro- 
vide the  desired  flow  through  the  test  section.   All  high 
point  bleeder  valves  were  opened  and  trapped  gasses 
allowed  to  escape.  At  this  point,  the  guard  heaters  were 
turned  on  and  adjusted  to  provide  a  insullating  block  sur- 
face temperature  equivalent  to  the  expected  coolant  channel 
wall  temperature.   The  bus  bar  cooling  water  was  turned  on 
as  well  as  the  city  water  to  the  heat  exchanger.   The 
final  stage  in  the  startup  process  was  accomplished  by 
synchronizing  and  starting  the  generators.   After  allowing 
the  generators  to  warm  up,  they  were  brought  up  to  power 
and  allowed  to  stabilize.   The  flow  of  city  water  to  the 
heat  exchanger  was  adjusted  until  a  constant  inlet  coolant 
water  temperature  to  the  test  section  was  obtained.   With 
all  conditions  stable,  data  was  taken.   This  data  included 
coolant  channel  inlet  and  exit  pressure,  voltage  drop 
across  and  amperage  flowing  through  the  test  section,  and 
all  thermocouple  readings. 
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5%4.5  Assumptions  for  Data  Analysis 

In  analyzing  coolant  channel  heat  transfer  data,  the 
Colburn  correlation   Nuf  =  0.023  Ref •   Prf0,3   was 
used  as  the  basis  for  comparison.   From  the  basic  equation 
for  heat  conduction   Q  =  hA(T  -Tb)    ,  Q,  T  ,  and  Tfe  are 
experimentally  measured,  and  A  is  the  constructed  internal 
surface  area  of  the  coolant  channel.   The  heat  transfer 
coefficient  can  thereby  be  calculated  directly.   For  pur- 
poses of  comparing  smooth  channel  experimental  data  with 

the  Colburn  correlation,  equivalent  diameter  was  calculated 

4A 
using  the  relation,  D  =  -p-  . 

In  order  to  compare  heat  transfer  from  a  smooth  chan- 
nel with  that  from  a  finned  channel,  it  is  necessary  to 
establish  certain  constants  for  comparison.   For  the  pur- 
poses of  reactor  design,  it  is  both  desirable  and  reasonable 
that  both  the  center-to-center  spacing  of  the  fuel  meat  in 
the  core,  and  the  mass  flow  rate  through  an  individual 
coolant  channel  remain  constant.   Consequently,  in  compar- 
ing the  smooth  channel  with  the  finned  channel,  it  was 
assumed  that  the  two  channels  had  both  the  same  equivalent 
diameter  and  the  same  cross  sectional  area  for  coolant  flow. 
Thus,  for  a  given  mass  flow  rate  Reynolds  number  (smooth) 
equals  Reynolds  number  (finned).   As  a  result,  a  relative 
comparison  of  heat  transfer  ability  was  obtained. 

The  only  temperature  correction  incorporated  into  the 
experimental  data  analysis  was  the  correction  for  temperature 
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drop  apross  the  coolant  channel  wall.   Coolant  channel 
outside  wall  temperature  is  measured  whereas  inner  wall 
temperature  is  desired.   The  correction  was  made  using 
the  relation, 

t    —  T      fo 
^wi  "  ±wo  "  2KA  ' 

where  Q/A  =  heat  flux 

K   =  thermal  conductivity  of  the  channel  wall,  and 

t   =  channel  wall  thickness. 
Although  to  a  first  approximation  the  heat  flux  along  the 
coolant  channel  was  uniform,  a  small  correction  was  required 
to  the  heat  flux  as  a  function  of  position  along  the  coolant 
channel  due  to  the  temperature  dependence  of  resistivity. 
Since  resistivity  increases  with  increasing  temperature 
according  to  the  relation  p  =  p  (1+aT),  more  heat  would  he 
generated  near  the  exit  of  the  electrically  heated  coolant 
channel  than  near  the  entrance.   The  temperature  coeffici- 
ent of  resistivity  for  a  piece  of  the  stainless  steel 
used  in  the  fabrication  of  the  coolant  channels  was  measured 
over  a  range  of  temperature  between  that  of  liquid  nitrogen 
and  100°F  water  and  found  to  be  0.002/°F.   Wall  temperature 
rise  between  channel  inlet  and  exit  was  typically  on  the 
order  of  50°F.   Therefore,  the  actual  input  heat  flux  to 
the  channel  varied  typically  from  -5%   at  the  entrance  to 
+5%   at  the  exit  from  the  channel  average.   Since  for  a  given 
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measured  temperature  drop  between  the  channel  wall  and 
bulk  fluid  heat  transfer  coefficient  is  directly  propor- 
tional to  input  heat  flux,  a  corresponding  correction  was 
made  to  the  experimentally  measured  heat  transfer  coeffici- 
ents. 

5.4.6  Experimental  Results 

Smooth  channel  data  is  presented  in  Figure  5.4.6.1 
and  is  shown  to  be  about  12$  above  the  Colburn  correlation, 
Such  increased  heat  transfer  was  expected  because  the 
smooth  channel  tested  was  not  perfectly  smooth.   When  the 
thermocouples  were  pressed  onto  the  channel  outside  wall 
small  indentations  were  formed  which  protruded  through  the 
wall  and  into  the  coolant  flow  path.   These  thermocouple 
"nipples"  acted  as  roughness  elements  in  promoting  heat 
transfer.   The  smooth  channel  data  taken  is  best  fit  by 
the  correlation  Nuf  =  0.027  Ref0,8  Prf0,3  .   (It  should  be 
noted  that  a  leak  developed  at  the  entrance  solder  joint 
between  the  test  section  end  block  and  coolant  channel. 
It  is  felt  that  this  leak  did  not  affect  experimental 
results,  however.) 

Finned  channel  data  is  presented  in  Figures  5.4.6.2, 
5.4.6.3,  5.4.6.4  and  5.4.6.5.   The  upper  curve  of  Figure 
5.4.6.2  was  drawn  through  points  representing  the  average 
value  of  calculations  from  10  centerline  thermocouple 
readings  for  each  experimental  run.   The  curve  is  fit  by 
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the  relation  Nuf  =  0.052  Ref0*8  Prf0,3  ,   The  points 
along  the  lower  curve  of  Figure  5.^.6.2  represent  average 
values  for  each  smooth  channel  experimental  run.   In 
calculating  the  heat  transfer  coefficients  for  the  finned 
channel  experimental  runs,  it  was  observed  that  the  heat 
transfer  coefficient  was  highest  near  the  entrance  of  the 
coolant  channel  and  decreased  toward  the  exit  of  the 
channel.   This  effect  did  not  show  up  in  the  smooth  channel 
runs.   Heat  transfer  coefficient  averaged  over  the  seven 
experimental  runs  is  shown  plotted  versus  thermocouple 
position  along  the  coolant  channel  in  Figure  5.^.6.3. 
Since  one  would  normally  expect  the  heat  transfer  coeffici- 
ent to  increase  slightly  toward  the  exit  of  the  coolant 
channel,  the  validity  of  the  finned  channel  heat  transfer 
data  was  subject  to  investigation.   It  was  suspected  that 
some  asymmetry  existed  either  in  the  thermocouple  readings 
or  in  the  channel  itself.   In  order  to  investigate  this 
possibility  the  test  section  was  inverted.    Everything 
remained  as  it  was  except,  now,  upside-down,  and  an  entire 

new  set  of  experimental  runs  was  made.   Average  values  of 

0  3 
Nu/Pr  *   are  shown  plotted  in  Figure  5.4.6.4  and  are  fit 

by  the  relation  Nuf  =  0.050  Ref0,8  Prf°*3  ,  almost  exactly 

the  same  as  the  relation  which  fit  the  original  data. 

Heat  transfer  coefficient  versus  position  along  the  coolant 

channel  is  plotted  in  Figure  5.4.6.5.   Although  the  slope 

of  this  curve  is  not  as  great  as  that  of  Figure  5.4.6.3, 
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the  same  trend  of  decreasing  heat  transfer  coefficient  with 
distance  along  the  coolant  channel  is  present,  and  the 
shape  of  the  curves  is  identical.   At  present  there 
appears  to  be  no  good  explanation  for  this  effect. 
However,  it  is  felt  that  the  average  value  curves  of  Fig- 
ures 5. 4.6.2  and  5.4.6.4  can  "be  used  with  confidence  in 
heat  transfer  calculations  for  the  proposed  High  Flux  MIT 
Reactor.   The  average  values  are  conservative  when  one 
considers  that  due  to  control  rod  insertion  into  the  reac- 
tor core,  only  the  lower  12-16  inches  of  the  fuel  plates 
will  be  generating  a  significant  amount  of  power,  and  it 
is  in  this  region  near  the  coolant  channel  entrance  that 
the  heat  transfer  coefficient  was  measured  to  be  above 
average. 

In  computing  average  values  for  the  finned  channel 
data,  only  the  readings  from  thermocouples  along  the  chan- 
nel centerline  were  used.   Heat  transfer  near  the  edge  of 
the  finned  channel  was  observed  to  be  much  greater  (by 
about  a  factor  of  2)  than  heat  transfer  along  the  channel 
centerline.   The  most  probable  reason  for  this  effect  is 
that  there  existed  a  small  gap  between  the  channel  side 
walls  and  the  ends  of  the  transverse  fins  in  which 
secondary  (additional)  turbulence  was  generated. 

It  is  recommended  that  the  modified  Colburn  correla- 
tion, Nuf  -  0.023  Ref0,8  Prf0,3  ,  be  used  for  smooth 
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channel  calculations  and  that  the  correlation  Nu~  =  0.050 
Ref  "   Prf    be  used  for  calculations  for  finned  channels 
similar  to  that  of  this  experiment.   Heat  transfer  from 
the  finned  channel  is  thus  a  factor  of  2.2  greater  than 
heat  transfer  from  a  smooth  channel.  The  heat  transfer 
surface  area  of  the  finned  channel  was  a  factor  of  I.67 
greater  than  that  of  the  smooth  channel.   Since  fin 
efficiency  in  this  case  was  nearly  100$,  the  factor  of  2.2 
increase  in  heat  transfer  can  be  taken  as  the  product  of 
a  factor  of  I.67  increase  in  heat  transfer  area  multiplied 
by  a  factor  of  1.3  due  to  increased  turbulence  caused  by 
the  fins.   To  be  more  precise,  the  product  "hA"  was  in- 
creased by  a  factor  of  2.2,  while  "h"  alone  was  increased 
by  a  factor  of  1.3. 

Pressure  drop  data  for  both  the  smooth  and  finned 
channels  is  presented  in  Figure  5.^.6.6.   The  smooth 
channel  pressure  drop  was  much  larger  than  that  normally 
associated  with  a  smooth  channel,  but,  like  smooth  channel 

heat  transfer,  this  can  be  attributed  to  channel  roughness 

1 

caused  by  thermocouple  indentations.   Channel  entrance  and 
exit  pressure  losses  are  included  in  the  plotted  value  of 
smooth  channel  pressure  drop. 

Finned  channel  pressure  drop  was  very  large,  with  AP 

-1  r 

varying  roughly  proportional  to  V  *  .   This  large  pressure 
loss  is  a  limiting  factor  on  the  coolant  mass  flow  rate 
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which  ean  be  pumped  through  the  core  using  the  existing 
pumps . 

5.4.7  Analysis  of  Results 

D20  pressure  measurements  around  the  primary  coolant 
loop  of  the  existing  MITR  show  that  there  exists  a  23 
psig  pressure  drop  between  the  primary  heat  exchanger 
exits  and  the  core  inlet  piping  (3&).      The  pressure  loss 
in  the  piping  between  these  two  points  should  be  no  more 
than  1-2  psi.  The  23  psig  pressure  drop  can  be  explained 
by  the  fact  that  between  the  heat  exchanger  exits  and  the 
core  inlet  are  valves  which  are  kept  partly  closed  in 
order  to  regulate  the  primary  coolant  flow.   Were  these 
valves  completely  opened,  it  is  reasonable  to  assume  that 
an  additional  pressure  head  on  the  order  of  21  psig  would 
be  available.   Consequently,  core  pressure  drop  could  be 
increased  by  approximately  21  psig  while  maintaining  a 
flow  rate  of  2000  gpm  with  the  existing  primary  coolant 
pumps.   However,  a  pressure  drop  of  21  psig  through  a 
finned  channel  corresponds  to  a  flow  rate  of  only  4  gallons 
per  minute.   Were  10.95  gpm  to  be  pumped  through  a  smooth 
channel,  equivalent  heat  transfer  would  be  realized. 
Implicit  in  the  above  information  is  the  fact  that  there 
exists  a  certain  cutoff  core  size,  assuming  no  additional 
pumping  power  is  to  be  installed. 
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Smooth  coolant  channels  will  give  better  heat  in 
cores  of  smaller  size,  whereas  it  is  advantageous  to  use 
the  finned  coolant  channels  in  cores  of  larger  size. 
Table  5.4.7.1  below  gives  a  brief  analysis  of  two  coolant 
channel  sizes  of  possible  interest. 

Table  5.4.7.1 


2.5"x0.090M 
channel 


2.1"x0.090" 
channel 


Max.  finned  channel 
flow  rate* 


4  gpm 


3.36  gpm 


Equivalent  smooth 
channel  flow  rate 


10.95  gpm 


9.20  gpm 


Break  even  point 


183  channels,  or 
11-17  channel 
elements 


218  channels,  or 
15-15  channel 
elements 


Actual  heat  transfer 

t««v.«~o~  hAffinned 
Increase 


hA( smooth 


1.2  for  a  12 
element  core 


1.84  for  a  26 
element  core 


*Core  AP    =21  psi 
max      e 


Thus,  for  purposes  of  the  MITR  redesign,  the  benefit  to  be 
derived  from  fuel  plate  finning  varies  with  the  size  of 
the  core  under  consideration. 

One  core  configuration  considered  for  the  phase  II 
High  Flux  MITR  contains  26  fuel  elements,  with  each  element 
containing  15-2.1"  x  0.090"  coolant  channels.   Comparing 
smooth  channels  to  finned  channels  for  this  core  configura- 
tion shows  that  the  pressure  drop  limited  finned  channel 
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will  provide  a  factor  of  1.84  better  heat  transfer  than  a 
smooth  channel.   However,  since  the  heat  transfer  area  of 
the  finned  channel  is  a  factor  of  I.67  greater  than  the 
smooth  channel,  the  heat  transfer  coefficient  is  in  fact 
increased  by  only  a  factor  of  1.1.   More  important, 
though,  is  the  fact  that  the  factor  increase  of  1.84  is 
greater  than  the  factor  of  1.8  required  for  the  light  water 
core.   While  the  optimum  plate  roughness  condition  for 
heat  transfer  augmentation  has  not  yet  been  achieved,  it 
has  been  demonstrated  that  the  currently  installed  process 
system  is  capable  of  removing  the  heat  from  the  hottest 
fuel  element  of  the  proposed  High  Flux  MITR  core. 

For  a  flow  of  2000  gpm  through  a  26  element  core, 
5.15  gpm  would  flow  through  each  coolant  channel.   From 
Figure  5.4.6.6  this  flow  rate  would  cause  a  pressure  drop 
across  the  core  on  the  order  of  33  psi,  or  12  psi  above 
the  existing  system  capability.   However,  if  one  were  to 
either  install  new  main  coolant  pumps  or  install  booster 
pumps  downstream  of  the  heat  exchangers  capable  of  provid- 
ing an  additional  12  psi  of  pressure  head,  the  finned 
channel  pressure  drop  limitation  could  be  eliminated,  thus 
allowing  one  to  fully  realize  the  2.2  factor  increase  in 
heat  transfer  for  a  finned  channel  over  a  corresponding 
smooth  channel.   In  addition,  if  new  coolant  pumps  were  to 
be  purchased,  pumps  of  not  only  higher  head  capacity  but 
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also  greater  flow  rate  capacity  could  be  installed,  thus 
adding  even  a  little  more  to  the  degree  of  heat  transfer 
augmentation  possible. 

5.4.8  Further  Investigations 

At  the  suggestion  and  with  the  assistance  of  this 
author,  the  heat  transfer  and  pressure  drop  characteris- 
tics of  a  longitudinally  finned  coolant  channel  have  been 
investigated  by  Mr.  Donald  Uhl.   A  longitudinally  finned 
(fins  running  parallel  to  the  coolant  flow  path)  fuel 
plate  enjoys  three  advantages  over  a  transverse  finned 
(fins  perpendicular  to  the  coolant  flow  path)  fuel  plate 
such  as  investigated  in  the  previous  sections:   (1)  the 
longitudinal  fins  provide  little  additional  turbulence  and 
thereby  create  less  pressure  drop  across  the  core,  (2) 
fabrication  of  the  fins  is  less  time  consuming  because 
fewer  milling  cuts  are  required,  and  (3)  fuel  plate  buck- 
ling cannot  completely  block  a  coolant  channel  because 
flow  can  always  occur  in  the  grooves  between  fins. 

It  was  originally  intended  that  the  coolant  channel 
plate  be  milled  so  as  to  produce  fins  20  mils  wide,  20 
mils  in  height,  and  spaced  20  mils  apart,  thus  increasing 
surface  area  over  that  for  a  smooth  plate  by  a  factor  of 
2.   However,  due  to  difficulties  experienced  in  milling 
the  fins  into  the  stainless  steel  plate  (difficulties 
which  should  not  be  encountered  when  milling  fins  into 
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the  actual  aluminum  fuel  plates),  the  spacing  between 
fins  was  increased  from  20  mils  to  40  mils,  thus  providing 
an  area  increase  of  I.67. 

Pressure  drop  data  versus  flow  rate  for  the  longitu- 
dinally finned  coolant  channel  is  shown  on  Figure  5.^.6.6. 
Although  larger  than  that  for  a  smooth  channel,  pressure 
drop  is  shown  to  be  considerably  less  than  that  for  a 
transverse  finned  channel  and  will  not  impose  a  limiting 
condition  on  primary  coolant  flow  rate. 

Experimental  heat  transfer  data  is  presented  in 
Figure  5.^.8.1.   Each  point  shown  is  an  average  of  the 
data  from  each  of  the  10  thermocouple  readings  recorded 
for  each  run.   The  experimental  data  is  fit  by  the  rela- 
tion Nuf  =  0.037  Ref0,8  Prf0,3.   Thus  the  longitudinally 
finned  channel  augments  heat  transfer  by  a  factor  of  1.6 
over  that  from  a  smooth  plate.   Since  the  area  ratio, 
finned  to  smooth,  was  I.67,  it  is  seen  that  the  increase 
in  heat  transfer  was  0.958  times  the  increase  in  surface 
area.   Thus,  it  is  reasonable  to  expect  that  if  the 
longitudinal  fins  were  to  be  spaced  20  mils  apart,  thereby 
doubling  the  heat  transfer  area,  one  could  expect  heat 
transfer  to  be  increased  by  a  factor  of  1.92. 
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Chapter  6 
CONCLUSIONS 

It  is  not  intended  that  the  work  which  has  been 
presented  represent  a  final  design  or  a  final  design 
analysis,  but  rather  that  the  preceeding  work  make  some 
contribution  toward  that  final  design.   Chapter  2  has 
shown  that  while  much  work  remains  to  be  done  in  the 
area  of  mechanical  design,  there  seem  to  be  no  prohibi- 
tively difficult  or  unsurmountable  design  problems. 
While  the  shield  design  and  supporting  calculations  of 
Chapter  3  must  certainly  be  refined  based  on  final 
shielding  requirements,  it  has  been  shown  that  shielding 
will  not  pose  a  difficult  design  problem.   Chapter  4  has 
indicated  the  need  to  find  a  means  of  augmenting  the  con- 
vective  heat  transfer  from  the  High  Flux  MITR  fuel  plates. 
In  response,  the  experimental  work  of  Chapter  5  has  demon- 
strated that  by  modification  of  the  fuel  plate  surface  the 
required  heat  transfer  augmentation  can  be  achieved. 

In  surveying  the  work  pertinent  to  the  contents  of 
this  thesis  which  must  be  accomplished  prior  to  setting 
a  final  High  Flux  MITR  design,  it  is  apparent  that  mechani- 
cal design  work  is  being  thoroughly  covered  by  Mr.  E. 
Barnett.   Final  shield  design  can  be  obtained  by  simple 
extrapolation  of  the  work  of  Chapter  3.   It  is  in  the 
area  of  core  heat  removal  that  the  most  intriguing  area 
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for  future  work  lies,  that  of  optimizing  the  means  by 
which  fuel  plate  heat  transfer  is  to  be  augmented.   In 
order  to  make  a  tentative  decision  based  on  information 
available  at  this  time,  the  following  summary  analysis 
based  on  heat  transfer  coefficient,  h,  is  presented  for 

the  phase  II  core  design. 

13 
For  a  maximum  homogeneous  power  density  of  1.3  x  10 

o 

fissions/cm  -sec  along  the  base  of  the  core  (equivalent 
to  that  assumed  in  previous  analyses),  the  maximum  heat 
flux  through  fuel  plates  spaced  0.090"  apart  will  be 
2.12  x  10^  BTU/hr-ft2.   At  the  coolant  entrance,  the 
saturation  temperature  of  HpO  is  228. 5°F.  Assuming  the 
reactor  were  just  about  to  scram  due  to  high  coolant  out- 
let temperature,  bulk  inlet  coolant  temperature  would  be 
ll4.5°F.   Thus,  since  fuel  plate  wall  temperature  is  the 
sum  of  bulk  coolant  temperature  plus  the  film  temperature 
rise  from  the  bulk  fluid  to  the  fuel  plate  wall  (T  =  T,  + 
Q/hA) ,  in  order  to  prevent  boiling  the  quantity  Q/hA  must 
be  less  than  ll4°F.   Employing  hot  channel  factors  and 
solving  for  heat  transfer  coefficient  yields  a  maximum 
required  value  of  h  =  3750  BTU/hr-ft2-°F.   Were  hot 
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channel  factors  not  employed,  the  required  value  of  h 
would  have  been  1870  BTU/hr-ft2-°F.   From  the  experimental 
results  of  Chapter  5  at  a  flow  rate  of  2000  gpm,  "h"  is 
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calculated  to  be  2430  BTU/hr-ft  -  F  for  a  smooth  channel, 
3880  BTU/hr-ft2-°F  for  a  longitudinally  finned  channel 
(fins  40  mils  apart),  and  5300  BTU/hr-ft2-°F  for  a  trans- 
verse finned  channel.   Repeating  in  tabular  form: 

^reauired  (wi-tn  not  channel  factors)         3750 
h    .   ,  (without  hot  channel  factors)      1870 

hsmooth  channel  2^30 

hlongitudinal  fins  <sPaced  4°  mils  aPart)  388° 

longitudinal  fins  (sPaced  20  mils  aPart)  466° 

transverse  fins  (requires  booster  pump)  5300 

It  is  apparent  that  if  hot  channel  factors  are  not  employ- 
ed even  the  smooth  channel  is  capable  of  providing  adequate 
heat  removal.   The  ratio 

i_ 
smooth  channel    _  ,  ~ 


required  without  H.C.F. 


However,  a  safety  factor  of  1.3  is  not  considered  adequate 
due  to  the  uncertainties  involved  in  this  "first -of -a-kind" 
design. 

On  the  basis  of  the  above  analysis,  it  is  recommended 
that  longitudinal  fins  be  milled  into  the  fuel  plates  of 
the  High  Flux  MITR.   The  fins  should  be  20  mils  in  height, 
40  mils  in  width,  and  should  be  spaced  20  mils  apart. 
Should  milling  difficulty  require  that  the  fin  spacing  be 
increased,  the  spacing  may  be  increased  up  to  but  not 
beyond  40  mils. 
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,  Appendix  A 

EXPERIMENTAL  DATA 
In  the  experimental  data  which  follows,  the  run 
numbers  beginning  with  MA"  are  the  smooth  channel  runs, 
run  numbers  beginning  with  "B"  are  the  initial  transverse 
finned  channel  runs,  and  run  numbers  beginning  with  "c" 

are  the  transverse  finned  channel  runs  after  the  channel 

0  3 
was  inverted.  Average  Re  and  average  Nu/Pr  *   values 

are  an  average  of  calculations  based  on  centerline 

thermocouple  readings  only  for  the  finned  channel  runs. 


Data  Sheet 
Run  no.  A-l 
Flow  Rate  =  1.66  gpm 
Electrical  Power  =  43,900  BTUAr 
Q/A  =  52,800  BTU/hr-ft2 
Inlet  Pressure  =6.1  psi 
Exit  Pressure  =4.6  psi 
Inlet  Temperature  =  98°F 
Exit  Temperature  =  151.5°F 

Thermo-  <p     T     T 
couple    wo    wi    b 
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AP  =  1.5  psi 


h 


Re    Nu/Pr0'3 


2 

201 

197 

145 

172 

1010 

8430 

27.9 

2R 

193 

189 

145 

167 

1200 

8150 

32.6 

3 

192 

188 

138 

163 

1050 

7920 

28.5 

4 

190 

186 

131 

159 

960 

7760 

25.8 

4R 

186 

182 

131 

157 

1030 

7590 

27.6 

5 

173 

169 

125 

147 

1200 

7642 

31.5 

5R 

166 

162 

125 

144 

1420 

6920 

37.3 

6 

170 

166 

125 

146 

1285 

6950 

33.6 

7 

174 

170 

125 

148 

1170 

7760 

30.6 

7R 

169 

I65 

'125 

145 

1320 

6930 

34.5 

8 

165 

161 

118 

139 

1240 

6650 

32.1 

8r 

168 

164 

118 

141 

1145 

6660 

29.8 

9 

177 

173 

111 

142 

860 

6680 

22.3 

10 

157 

152 

105 

129 

1120 

6l45 

28.5 

10R 

166 

162 

105 

134 

920 

6390 

23.7 

Average 

Re  = 

7230 

Average 

Nu/Pr 

,0.3 

29.7 
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Data  Sheet  2 

Run  no.  A-2 

Flow  Rate  =2.5  gpm 

Electrical  Power  =  60,460  BTU/hr 

Q/A  =  72,550  BTU/hr-ft2 

Inlet  Pressure  =  9.0  psi 

Exit  Pressure  =  7.0  psi 

Inlet  Temperature  =  91°F 

Exit  Temperature  =  139. 5°F 

Thermo-  T 

couple    wo    wi    b 


AP  =  2  psi 


Re 


Nu/Pr0'3 


'r 

39.4 

53.2 
42.8 
32.8 
39.4 
37.9 
41.9 
38.2 
39.4 

34.7 
41.8 
38.1 
35.4 
37.8 
30.6 


2 

2R 

3 

4 

4R 

5 

5R 

6 

7 
7R 

8 
8r 

9 

10 

10R 


189 
176 

179 
186 
176 

171 
166 

171 
169 

175 
160 
164 
161 
153 
163 


183 
170 
173 
180 
170 
165 
160 
165 
163 
169 
154 
158 
155 
147 
157 


134 
134 
128 
122 
122 
115 
115 
115 
115 

115 

1 

109 

109 

103 

97 

97 


158 
152 
150 
150 
146 
140 
138 
140 
140 
142 
132 
134 

129 
122 
127 


1475 
2010 
1620 
1240 
1510 
1470 
1630 
1480 
1520 
1345 
1640 
1490 
1400 
1465 
1210 


11400 
11000 

10900 
10900 
10600 
10050 
9820 
10050 

10000 

10200 

9400 

94oo 

9300 

8650 
9050 


Average  Re  =  10,050 
Average  Nu/Pr0,3  =38.9 
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Data  Sheet  3 
Run  no.  A-3 
Flow  Rate  =  3.4  gpra 
Electrical  Power  =  53,860  BTU/hr 
q/A  =  64,630  BTU/hr-ft2 
Inlet  Pressure  =  15.2  psi 
Exit  Pressure  =  13.5  psi 
Inlet  Temperature  =  101°F 
Exit  Temperature  =  133°F 
Thermo- 


AP  =  1,7  psi 


J.  liC  J.  iuv-»  — 

couple 

T 
wo 

T  . 

Wl 

Tb 

T 

h 

Re 

Nu/Pr0,3 

2 

165 

159 

129 

144 

2170 

14700 

56.3 

2R 

161 

156 

129 

143 

2460 

13650 

64.0 

3 

157 

152 

125 

139 

2430 

13200 

62.5 

4 

158 

153 

121 

137 

2060 

13200 

53.1 

4R 

156 

151 

121 

136 

2200 

13000 

58.4 

5 

149 

144 

117 

130 

2450 

13200 

62.2 

5R 

151 

145 

117 

131 

2300 

12540 

58.8 

6 

150 

144 

117 

131 

2400 

13200 

61.0 

7 

145 

l4o 

117 

128 

2940 

12500 

74.7 

7R 

161 

156 

117 

137 

1670 

13200 

43.1 

8 

142 

137 

113 

125 

2740 

11800 

69.O 

8r 

154 

148 

113 

131 

1840 

12500 

46.6 

9 

141 

136 

109 

123 

2410 

11670 

60.5 

10 

135 

130 

105 

117 

2620 

11100 

65.1 

10R 

148 

143 

105 

124 

1705 

11760 

43.0 

Average 

Re  =  12 

,740 

Average 

Nu/Pr°* 

3  =  58 

.6 

Data  Sheet 

Run  no.  A -4 

Flow  Rate  =  3.5  gpm 

Electrical  Power  =  73,090  BTUAr 

Q/A  =  87,770  BTU/hr-ft2 

Inlet  Pressure  =  l4  psi 

Exit  Pressure  =  11.6  psi 

Inlet  Temperature  =  91°F 

Exit  Temperature  =  133. 5°F 

Thermo-  m 

couple    wo    wi 
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AP  =  2.4  psi 


Re    Nu/Pr0,3 


2 

170 

163 

128 

146 

2510 

15700 

65.5 

2R 

165 

158 

128 

143 

2920 

i44oo 

76.0 

3 

162 

155 

123 

139 

2740 

i4ooo 

71.4 

4 

163 

156 

118 

137 

2320 

13600 

60.6 

4r 

159 

151 

118 

139 

2700 

13970 

70.3 

5 

153 

146 

112 

129 

2580 

12850 

65.5 

5R 

154 

147 

112 

130 

2510 

12900 

63.6 

6 

154 

147 

112 

130 

2510 

12900 

63.6 

7 

156 

149 

112 

131 

2370 

12900 

60.6 

7R 

162 

155 

1 

112 

134 

2040 

13250 

52.3 

8 

148 

141 

107 

124 

2580 

12300 

65.O 

8R 

157 

150 

107 

129 

2040 

12850 

51.8 

9 

146 

139 

102 

121 

2370 

11800 

59.2 

10 

143 

136 

96 

116 

2190 

11340 

54.4 

10R 

152 

145 

96 

121 

1790 

11800 

44.7 

Average  Re  =  13,100 
Average  Nu/Pr0'3  =61.6 


Data  Sheet  5 
Run  no,  A -5 
Flow  Rate  =4.0  gpm 
Electrical  Power  =  62,530  BTU/hr 
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Q/A  =  75,040  BTU/hr -ft 2 
Inlet  Pressure  =  16.7  psi 
Exit  Pressure  =  13.7  psi 
Inlet  Temperature  =  100°F 
Exit  Temperature  =  131. 5°F 


AP  =  3.0  psi 


Thermo- 
couple 

T 
wo 

T  . 

Wl 

Tb 

Tf 

h 

Re 

Nu/Pr0'3 

2 

166 

160 

128 

144 

2345 

16640 

61.3 

2R 

162 

156 

128 

142 

2680 

16360 

69.7 

3 

160 

154 

124 

139 

2500 

16080 

64.7 

4 

159 

153 

120 

137 

2270 

15700 

58.6 

4r 

155 

149 

120 

135 

2590 

15460 

66.3 

5 

150 

144 

116 

130 

2680 

14870 

68.1 

5R 

151 

145 

116 

131 

2590 

14870 

65.7 

6 

150 

144 

116 

130 

2730 

14370 

69.3 

7 

150 

144 

116 

130 

2730 

14870 

69.3 

7R 

159 

153  , 

116 

135 

2030 

15460 

52.0 

8 

146 

l4o 

112 

126 

2680 

14530 

67.4 

8r 

153 

147 

112 

130 

2140 

14870 

54.5 

9 

146 

l4o 

108 

124 

2345 

14020 

58.7 

10 

142 

136 

104 

120 

2340 

13620 

58.5 

10R 

150 

144 

104 

124 

1900 

14020 

47.9 

Average 

Re  =  15 

,080 

Average 

Nu/Pr0, 

3  =  62 

.1 

Data  Sheet  6 
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Run  no*.  A -6 

Flow  Rate  =4.2  gpm 

Electrical  Power  =  55,500  BTU/hr 

Q/A  =  66,590  BTU/hr -ft2 

Inlet  Pressure  =  15.8  psi 

Exit  Pressure  =  14  psi 

Inlet  Temperature  =  119  F 

Exit  Temperature  =  l46°F 

Thermo-   t     T     T 
couple     wo    wi    h 


AP  =  1.8  psi 


h 


Re    Nu/Pr 


Average  Re  =  18,990 

0  3 
Average  Nu/Pr    =  65.5 


0.3 


2 

177 

172 

143 

157 

2230 

20500 

59.8 

2R 

168 

163 

143 

152 

3300 

19700 

87.I 

3 

172 

167 

139 

153 

2390 

19750 

63.1 

4 

172 

167 

136 

151 

2130 

19600 

56.4 

4r 

163 

158 

136 

147 

3010 

19000 

78.8 

5 

165 

160 

133 

146 

2480 

19200 

64.8 

5R 

157 

153 

133 

142 

3310 

18300 

86.2 

6 

165 

160 

133 

146 

2480 

19200 

64.8 

7 

165 

160 

133 

146 

2480 

19200 

64.8 

7R 

169 

164  ' 

133 

148 

2160 

19000 

56.3 

8 

165 

160 

129 

144 

2210 

19400 

57.7 

8R 

163 

158 

129 

143 

2360 

18300 

61.4 

9 

162 

157 

126 

141 

2180 

18200 

56.3 

10 

156 

151 

123 

136 

2350 

17600 

60.5 

10R 

154 

149 

123 

136 

2490 

17850 

63.9 

Data  Sheet 
Run  no*  A -7 
Flow  Rate  =  4.95  gpm 
Electrical  Power  =  86,460  BTU/hr 
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Q/A  =  103,750  BTU/hr-ftT 
Inlet  Pressure  =  23.2  psi 
Exit  Pressure  =  19.5  psi 
Inlet  Temperature  =  103°F 
Exit  Temperature  =  138. 5°F 


AP  =  3.7  psi 


Thermo- 
couple 

T 
wo 

T  . 

Wl 

Tb 

Tf 

h 

Re 

Nu/Pr0*3 

2 

180 

171 

134 

153 

2780 

22050 

74.0 

2R 

176 

168 

134 

151 

3110 

21600 

82.2 

3 

173 

164 

130 

147 

2990 

21050 

78.7 

4 

173 

164 

125 

145 

2650 

20700 

69.2 

4R 

168 

160 

125 

142 

3030 

19900 

78.7 

5 

161 

152 

121 

136 

3340 

19300 

86.0 

5R 

173 

164 

121 

142 

2410 

19900 

62.6 

6 

163 

155 

121 

138 

3060 

19450 

79.4 

7 

162 

153 

121 

138 

3230 

19400 

83.2 

7R 

173 

164  , 

121 

142 

2410 

19900 

62.6 

8 

160 

152 

116 

134 

3000 

19000 

75.9 

8R 

169 

160 

116 

138 

2390 

19500 

61.8 

9 

160 

152 

112 

132 

2630 

18400 

67.I 

10 

151 

143 

108 

125 

2960 

17800 

74.9 

10R 

159 

150 

108 

129 

2420 

18350 

62.1 

Average 

Re  =  19 

,750 

Average 

Nu/Pr0, 

3  =  73 

.2 
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Data  Sheet  8 
Run  no.  B-l 
Flow  Rate  =  1.38  gpm 
Electrical  Power  =  34,800  BTU/hr 
Q/A  =  41,800  BTU/hr -ft2 
Inlet  Pressure  =   9.5  psi 
Exit  Pressure  =  4.5  psi 
Inlet  Temperature  =  101°F 


AP  =   5  psi 


Exit  Temperature  = 

152°F 

Thermo- 
couple 

T 
wo 

Twi 

Tb 

T 

h 

Re 

Nu/Pr0,3 

2 

181 

178 

146 

162 

1300 

6650 

35.1 

2R 

175 

172 

146 

159 

1595 

6500 

42.8 

3 

169 

166 

139 

153 

1570 

6190 

41.7 

4 

158 

155 

133 

144 

19^0 

5760 

50.8 

4R 

157 

154 

133 

143 

2020 

5710 

52.6 

5 

139 

135 

127 

131 

4920 

5170 

125.4 

5R 

137 

134 

127 

130 

5970 

5140 

152.0 

6 

148 

145 

127 

136 

2300 

5^00 

59.1 

7 

142 

139 

127 

133 

3480 

5290 

89.2 

7R 

138 

134 

127 

130 

5430 

5l4o 

138.2 

8 

143 

139 

120 

130 

2200 

5l4o 

56.0 

9 

136 

133 

114 

123 

2235 

4820 

56.1 

10 

126 

123 

107 

115 

2770 

4490 

68.5 

10R 

127 

124 

107 

116 

2560 

4530 

63.6 

Average 

Re  =  5: 

,520 

Average 

Nu/Pr0, 

,3  =  52.6 

' 
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Data  Sheet  9 


Run  no.  B-2 

Flow  Rate  =  2.0  gpm 

Electrical  Power  =  48,400  BTU/hr 

Q/A  =  58,100  BTU/hr-ft2 

Inlet  Pressure  =  13.8  psi 

Exit  Pressure  =  6.5  psi 

Inlet  Temperature  =  94.5°F 

Exit  Temperature  =  l43°F 


AP  =  7.3  psi 


Thermo- 
couple 

T 
wo 

T  . 

Wl 

Tb 

Tf 

h 

Re 

Nu/Pr0,3 

2 

176 

171 

137 

154 

1720 

9140 

45.6 

2R 

173 

168 

137 

152 

1900 

8980 

50.3 

3 

165 

160 

131 

145 

2040 

8520 

53.4 

4 

152 

147 

125 

136 

2620 

7900 

67.3 

4R 

153 

148 

125 

137 

2470 

8000 

63. a 

5 

133 

128 

119 

123 

6455 

7050 

162 

5R 

133 

128 

119 

123 

6455 

7050 

162 

6 

146 

141 

119 

134 

2665 

7800 

68.4 

7 

137 

132 

119 

125 

4305 

7200 

108 

7R 

134 

129 

119 

124 

5480 

7140 

138 

8 

138 

133 

113 

123 

2890 

7060 

72.5 

9 

130 

125 

107 

116 

3160 

6630 

78.3 

10 

120 

115 

101 

108 

4010 

6100 

97.5 

10R 

124 

119 

101 

110 

3230 

6250 

79.2 

Average 

Re  =  7: 

,640 

Average 

Nu/Pr0, 

'3  =  67 

.6 

Data  Sheet  10 
Run  no*.  B-3 
Flow  Rate  =  2.6  gpm 
Electrical  Power  =  80,000  BTU/hr 
Q/A  =  96,100  BTU/hr-ft2 
Inlet  Pressure  =   19.5  psi 
Exit  Pressure  =     9  psi 
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AP  =10.5  psi 


Inlet  Temperature  = 

95"F 

Exit  Temperature  = 

157°F 

Thermo- 
couple 

T 
wo 

T  . 

Wl 

Tb 

T 

h 

Re 

Nu/Pr0*3 

2 

197 

189 

149 

169 

2440 

12900 

66.6 

2R 

194 

186 

149 

163 

2590 

12350 

70.7 

3 

185 

177 

142 

159 

2690 

12000 

72.3 

4 

171 

163 

134 

149 

3260 

11100 

85.8 

4R 

171 

163 

134 

149 

3260 

11100 

85.8 

5 

146 

138 

126 

132 

7810 

$650 

200 

5R 

148 

140 

126 

133 

7120 

9780 

182 

6 

158 

150 

126 

138 

4055 

10150 

105 

7 

152 

144 

126 

135 

5340 

9930 

137 

7R 

148 

140 

126 

133 

6720 

9780 

172 

8 

153 

145 

118 

132 

3560 

9650 

90.9 

9 

143 

135 

111 

123 

3990 

8900 

100 

10 

130 

122 

103 

112 

5005 

8050 

123 

10R 

135 

127 

103 

115 

4020 

8280 

99.5 

Average  Re  =  10,440 
Average  Nu/Pr0'3  =90.0 
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Data  Sheet  11 
Flow  no.  B-4 
Flow  Rate  =  3.0  gpm 
Electrical  Power  =  79,500  BTU/hr 
Q/A  =  95,400  BTU/hr -ft2 
Inlet  Pressure  =  25  psi 
Exit  Pressure  =  11.8  psi 


AP  =  13.2  psi 


Inlet  Temperature  = 

99"F 

Exit  Temperature  = 

152°F 

Thermo- 
couple 

wo 

T  . 

Wl 

\ 

T 

h 

Re 

Nu/Pr0,3 

2 

191 

183 

145 

I69 

2540 

15300 

69.2 

2R 

185 

177 

145 

161 

3020 

14350 

81.3 

3 

179 

171 

139 

155 

2940 

13800 

78.4 

4 

167 

159 

132 

146 

3530 

12850 

93.7 

4R 

165 

157 

132 

145 

3820 

12700 

99.8 

5 

146 

138 

126 

132 

7630 

11450 

195 

5R 

144 

136 

126 

131 

9090 

11300 

232 

6 

157 

149 

126 

137 

4150 

11950 

107 

7 

150 

142 

126 

134 

5780 

11650 

148 

7R 

147 

139 

126 

132 

7070 

11450 

181 

8 

151 

143 

119 

131 

3980 

11300 

101 

9 

143 

135 

112 

123 

4260 

10550 

107 

10 

133 

125 

106 

115 

5000 

9800 

124 

10R 

134 

126 

106 

116 

4680 

9900 

116 

Average 

Re  = 

12,250 

Average 

Nu/Pr 

0.3  = 

97.7 

Data  Sheet  12 

Run  no.  B-5 

Plow  Rate  =  3.5  gpm 

Electrical  Power  =  80,000  BTU/hr 

Q/A  =  96,100  BTU/hr- ft2 

Inlet  Pressure  =  32  psi 

Exit  Pressure  =  14.8  psi 

Inlet  Temperature  =  94.5°F 

Exit  Temperature  =  l4o.5°F 

Thermo-  T     T     m     , 
couple    wo    wi 
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AP  =  17.2  psi 


Average  Re  =  13,220 
Average  Nu/Pr0,3  =  100.4 


Re    Nu/Pr0,3 


2 

177 

I69 

135 

152 

2870 

15650 

75.9 

2R 

173 

I65 

135 

150 

3200 

15400 

84.6 

3 

167 

159 

129 

144 

3170 

14700 

83.0 

4 

156 

149 

123 

136 

3770 

13750 

97.1 

4R 

157 

149 

123 

136 

3700 

13750 

95.2 

5 

139 

131 

118 

124 

7120 

12450 

179 

5R 

i4o 

132 

118 

125 

6720 

12500 

169 

6 

150 

142 

118 

130 

3920 

13100 

99.9 

7 

142 

134 

118 

126 

5820 

12600 

147 

7R 

l4o 

132 

118 

125 

6630 

12500 

167 

8 

143 

135 

112 

124 

4110 

12450 

103 

9 

135 

127 

106 

117 

4580 

11650 

114 

10 

126 

118 

100 

109 

5490 

10800 

134 

10R 

130 

122 

100 

111 

4410 

10950 

118 
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Data  Sheet  13 
Run  no,  B-6 
Flow  Rate  =3.9  gpm 
Electrical  Power  =  79,400  BTU/hr 


Q/A  -  95,300  BTU/hr-ft^ 
Inlet  Pressure  =  38.5  psi 
Exit  Pressure  =  17.5  psi 
Inlet  Temperature  =  101°F 
Exit  Temperature  =  l4l°F 


AP  =   21  psi 


Thermo- 
couple 

T 
wo 

T  . 

Wl 

Th 

Tf 

h 

Re 

Nu/Pr0,3 

2 

179 

171 

136 

153 

2760 

17350 

73.3 

2R 

171 

163 

136 

149 

3600 

16900 

94.7 

3 

168 

160 

131 

145 

3310 

16350 

86.6 

4 

156 

148 

126 

137 

4330 

15350 

112 

4R 

156 

148 

126 

137 

4410 

15350 

114 

5 

i4o 

132 

121 

127 

8360 

14100 

212 

5R 

i4o 

132 

121 

126 

8910 

13900 

223 

6 

150 

142 

121 

132 

4540 

14700 

116 

7 

132 

135 

121 

128 

6810 

14200 

173 

7R 

141 

133 
1 

121 

127 

3080 

i4ioo 

204 

8 

145 

137 

116 

127 

4540 

i4ioo 

115 

9 

139 

131 

111 

121 

4890 

13300 

122 

10 

129 

121 

106 

114 

6230 

12500 

154 

10R 

133 

125 

106 

116 

5020 

12700 

124 

Average  Re  =  l4,86o 
Average  Nu/Pr0,3  =  111.1 


Data  Sheet  14 
Run  no.  B-7 
Flow  Rate  =  4.3  gpm 
Electrical  Power  =  80,000  BTU/hr 
Q/A  =  96,200  BTU/hr -ft2 
Inlet  Pressure  =  44  psi 
Exit  Pressure  =  20  psi 
Inlet  Temperature  =  97.5°F 
Exit  Temperature  =  135°F 
Thermo - 
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AP  =  24  psi 


j.uc.Lim-»- 

couple 

'  T 
wo 

T  . 

Wl 

Tb 

Tf 

h 

Re 

Nu/Pr0,3 

2 

163 

155 

130 

142 

3850 

17600 

100 

2R 

160 

152 

130 

141 

4370 

17500 

114 

3 

159 

151 

126 

138 

3850 

17100 

99.4 

4 

150 

142 

121 

131 

4510 

16100 

115 

4R 

149 

141 

121 

131 

4810 

16100 

123 

5 

137 

129 

116 

122 

7580 

14850 

190 

5R 

137 

129 

116 

123 

7400 

15000 

185 

6 

145 

137 

116 

126 

4740 

15400 

120 

7 

139 

131 

116 

123 

6420 

14950 

161 

7R 

137 

129 

114 

123 

7570 

14950 

190 

8 

140 

132 

112 

122 

4690 

14850 

118 

9 

134 

126 

107 

116 

5200 

14100 

129 

10 

126 

118 

102 

110 

6000 

13250 

147 

Average  Re  =  15,780 
Average  Nu/Pr0,3  =  118,2 
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Data 

Sheet  1 

5 

Run  no. 

1 

C-l 

Flow  Rate  =1.5 

1  gpm 

Electri< 

:al  Powe 

r  =  77,500 

BTU/hr 

Q/A  =  93,100  BTU/hr- 

ft2 

Inlet  Pressure 
Exit  Pressure 

=  9.7 
=  4.5 

psi 
psi 

AP  = 

5.2  ps 

i 

Inlet  Temperature  = 

97°F 

Exit  Temperature  = 

199°F 

Thermo- 
couple 

T 

wo 

T  . 

Wl 

Tb 

T 

h 

Re 

Nu/Pr0,3 

2 

250 

242 

186 

214 

1660 

10100 

44.8 

2R 

252 

244 

186 

215 

1605 

10100 

43.4 

3 

235 

227 

174 

200 

1760 

9330 

47.5 

4 

212 

204 

161 

183 

2160 

84oo 

58.4 

5 

178 

170 

148 

159 

4230 

7050 

114 

5R 

170 

162 

148 

155 

6650 

6800 

179 

6 

190 

182 

148 

165 

2740 

74oo 

74.0 

7 

170 

162 

148 

155 

6650 

6800 

179 

7R 

172 

164 

148 

156 

5800 

6900 

156 

8 

182 

174 

135 

155 

2390 

6800 

64.6 

8R 

178 

170 

135 

153 

2660 

6700 

71.8 

9 

I69 

161 

123 

142 

2450 

6150 

66.2 

10 

148 

i4o 

110 

125 

3150 

5330 

85.O 

10R 

152 

144 

110 

127 

2780 

5430 

75.0 

Average 

Re  =  7570 

Average 

Nu/Pr°* 

3  =  63.1 
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Data 

Sheet 

16 

Run  no.  C-2 

• 

Flow 

Rate  =  2. 

0  gpm 

Electrical  Pow 

er  =  ' 

78,900  BTU/hr 

Q/A  = 

9^,800  BTU/hr 

-ft2 

Inlet 
Exit 

Pressure 
Pressure 

=  13 
=  6.! 

.8  psi 
5  psi 

AP 

=  7.3 

psi 

Inlet 

Temperature  = 

98°F 

Exit 

Temperature  = 

176°F 

Thermo-  T 
couple    wo 

T  . 

Wl 

Tb 

Tf 

h 

Re 

Nu/Pr0,3 

2 

220 

212 

166 

189 

2060 

11600 

54.5 

2R 

220 

212 

166 

I89 

2060 

ll600 

54.5 

3 

208 

200 

157 

178 

2180 

10800 

57.8 

4 

192 

184 

147 

165 

2600 

9860 

68.9 

5 

164 

156 

137 

147 

4990 

8540 

130 

5R 

158 

150 

137 

144 

7300 

8350 

193 

6 

176 

168 

137 

153 

3050 

8950 

80.8 

7 

158 

150 

137 

144 

7300 

8350 

193 

7R 

160 

152 

137 

145 

6300 

8400 

167 

8 

169 

161 

127 
1 

144 

2790 

8350 

73.9 

8R 

165 

157 

127 

142 

3160 

8230 

83.7 

9 

159 

151 

118 

134 

2830 

7700 

75.0 

10 

142 

134 

108 

121 

3710 

6850 

98.3 

10R 

146 

138 

108 

123 

3160 

6970 

83.7 

Average  Re  =  9090 
Average  Nu/Pr0,3  = 


73.1 
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Data 

Sheet 

17 

Run  no.  C-3 

• 

Flow 

Rate  =  2 

.5  gpm 

| 

Elect 

rical  Power  = 

78,900 

BTU/h 

r 

1 

Q/A  = 

9^,800  BTU/hr 

-ft2 

Inlet 
Exit 

Pressure  =  19 
Pressure  =  9 

psi 
psi 

AP  = 

10  psi 

Inlet 

Temperature  = 

98°F 

Exit  ' 

Temperature  = 

l6l°F 

Thermo-  T 
couple   wo 

T  . 

Wl 

Tb 

Tf 

h 

Re 

Nu/Pr0'3 

2 

197 

I89 

153 

171 

2630 

12800 

68.4 

2R 

198 

190 

153 

171 

2560 

12800 

66.6 

3 

188 

180 

145 

163 

2710 

12100 

70.5 

4 

174 

166 

137 

152 

3270 

11100 

85.O 

5 

155 

147 

130 

138 

5410 

9960 

141 

5R 

151 

143 

130 

136 

7020 

9750 

183 

6 

166 

158 

130 

144 

3390 

10400 

88.1 

7 

151 

143 

130 

136 

7020 

9750 

183 

7R 

152 

143 

130 

137 

6530 

9850 

170 

8 

160 

152 

122 

137 

3160 

9850 

82.1 

8R 

157 

149 

122 

136 

3510 

9750 

91.2 

9 

151 

143 

114 

129 

3270 

9200 

85.O 

10 

137 

129 

106 

118 

4120 

8330 

107 

10R 

141 

133 

106 

120 

3510 

8500 

92.1 

Average  Re  =  10,480 
Average  Nu/Pr0,3  =83.6 


Data  Sheet  18 

Fun  no„  C-4 

Flow  Rate  =  3.0  gpm 

Electrical  Power  =  78,900  BTU/hr 

Q/A  =  94,800  BTU/hr-ft2 

Inlet  Pressure  =  25  psi 

Exit  Pressure  =  11.8  psi 

Thermo-  T     T     T 
couple   wo    wi    b 
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AP  =  13.2  psi 


Average  Re  =  11,800 
Average  Nu/Pr0*3  =92.1 


Re    Nu/Pr0,3 


2 

183 

175 

143 

159 

2960 

i4ioo 

75.5 

2R 

183 

175 

143 

159 

2960 

14100 

75.5 

3 

173 

165 

137 

151 

3380 

13200 

86.2 

4 

165 

157 

130 

144 

3510 

12500 

89.5 

5 

147 

139 

124 

131 

6110 

11200 

156 

5R 

144 

136 

124 

130 

7580 

11100 

193 

6 

156 

148 

124 

136 

3870 

11700 

98.6 

7 

144 

136 

124 

130 

7580 

11100 

193 

7R 

144 

136 

124 

130 

7580 

11100 

193 

8 

152 

144 

117 

131 

3510 

11200 

89.5 

8R 

149 

141 

117 

129 

3950 

11050 

101 

9 

145 

137 

110 

124 

3510 

10600 

89.5 

10 

133 

125 

104 

115 

4510 

9700 

115 

10R 

136 

128 

io4 

116 

3950 

9800 

101 
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Data 

Sheet 

19 

Run  no.  C-5 

• 

Flow 

Rate  =  3 

.5  gpm 

Elect 

rical  Power  = 

78,900 

BTU/h 

r 

Q/A  = 

94,800  BTU/hr 

-ft2 

Inlet 
Exit 

Pressure  =  32 
Pressure  =  14. 

psi 
8  psi 

AP 

=  17.2 

psi 

Inlet 

Temperature  = 

98°F 

Exit 

Temperature  = 

142°F 

Thermo-  T 
couple   wo 

T  . 

Wl 

Tb 

T 

h 

Re 

Nu/Pr0'3 

2 

174 

166 

137 

151 

3260 

15400 

83.1 

2R 

174 

166 

137 

151 

3260 

15400 

83.1 

3 

166 

158 

131 

145 

3505 

14700 

89.4 

4 

157 

149 

126 

137 

4030 

13850 

103 

5 

145 

137 

120 

129 

5560 

12900 

142 

5R 

142 

134 

120 

127 

6760 

12700 

172 

.6 

152 

144 

120 

132 

3950 

13250 

101 

7 

142 

134 

120 

127 

6760 

12900 

172 

7R 

141 

133 

120 

127 

7280 

12900 

185 

8 

148 

140 

114 

127 

3640 

12900 

92.8 

8R 

145 

137 

114 

125 

4120 

12450 

105 

9 

142 

134 

109 

122 

3790 

12100 

96.6 

10 

131 

123 

104 

113 

4850 

11100 

124 

10R 

134 

126 

104 

115 

4210 

11350 

107 

Average  Re  =  13,250 
Average  Nu/Pr0#3  =  98.4 


Data  Sheet  20 
Run  no.  C-6 
Flow  Rate  =  3.5  gpm 
Electrical  Power  =  137,300  BTU/hr 
Q/A  =  165,000  BTU/hr-ft2 
Inlet  Pressure  =  32  psi 
Exit  Pressure  =  14.8  psi 
Inlet  Temperature  =  93°F 
Exit  Temperature  =  170°F 
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AP 


17.2  psi 


Average  Re  =   15,250 
Average  Nu/Pr0,3  =   116.7 


Thermo- 
couple 

T 
wo 

T  . 

Wl 

Tb 

Tf 

h 

Re 

Nu/Pr0,3 

2 

219 

205 

160 

183 

3670 

19600 

96.5 

2R 

218 

204 

160 

182 

3750 

19500 

98.5 

3 

205 

191 

151 

171 

4130 

17950 

109 

4 

195 

181 

141 

161 

4130 

16700 

109 

5 

168 

154 

132 

143 

7340 

14450 

193 

5R 

163 

149 

132 

140 

9430 

14150 

248 

6 

179 

165 

132 

148 

4920 

13750 

130 

7 

164 

150 

132 

141 

8930 

14250 

235 

7R 

163 

149 

132 

140 

9440 

14150 

248 

8 

174 

160 

122 

141 

4350 

14250 

115 

8R 

170 

156 

122 

139 

4860 

14000 

128 

9 

164 

150 

112 

131 

4350 

13100 

115 

10 

148 

134 

103 

119 

5330 

11750 

l4o 

10R 

151 

137 

103 

120 

4860 

11900 

128 
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Data  Sheet  21 
Run  na.  C-7 
Flow  Rate  =  3.5  gpm 
Electrical  Power  =  35,^50  BTU/hr 
Q/A  =  42,550  BTU/hr-ft2 
Inlet  Pressure  =  32  psi 
Exit  Pressure  =  14.8  psi 
Inlet  Temperature  =  95°F 
Exit  Temperature  =  115. 5  F 
Thermo- 


AP  =  17.2  psi 


T     T 
couple   wo    wi 


h 


Re     Nu/Pr0,3 


2 

134 

131 

113 

122 

2430 

12100 

59.5 

2R 

134 

131 

113 

122 

2430 

12100 

59.5 

3 

130 

127 

110 

118 

2640 

11600 

64.8 

4 

125 

122 

108 

114 

3110 

11250 

76.1 

5 

121 

118 

105 

111 

3490 

10900 

85.4 

5R 

119 

116 

105 

110 

4170 

10800 

102 

6 

125 

122 

105 

113 

2630 

11100 

64.3 

7 

119 

116 

105 

110 

4170 

10800 

102 

7R 

118 

115 

105 

110 

4620 

10800 

113 

8 

123 

120 

103 

111 

2550 

10900 

62.4 

8R 

119 

116 

103 

109 

3350 

10700 

82.0 

9 

118 

115 

100 

107 

2955 

10500 

72.4 

10 

113 

110 

98 

104 

3580 

10200 

87.5 

10R 

115 

112 

98 

105 

3060 

10300 

75.0 

Average  Re  =  11,075 
Average  Nu/Pr0,3  =70.4 
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Data  Sheet  22 


Run  no.  C-8. 

Flow  Rate  =  4.0  gpm 

Electrical  Power  =  78,900  BTU/hr 

Q/A  =  94,800  BTU/hr-ft2 

Inlet  Pressure  =  39.5  psi 

Exit  Pressure  =  18  psi 

Inlet  Temperature  =  98°F 

Exit  Temperature  =  137°F 


AP  =  21.5  psi 


Thermo- 
couple 

T 
wo 

Twi 

Tb 

Tf 

h 

Re 

Nu/Pr0,3 

2 

165.5 

157.5 

132 

145 

3710 

16800 

93.5 

2R 

165.5 

157.5 

132 

145 

3710 

16800 

93.5 

3 

158 

150 

127 

139 

4120 

16000 

104 

4 

151.5 

143.5 

122.5 

133 

4500 

15300 

113 

5 

140 

132 

117.5 

125 

6530 

14200 

164 

5R 

138 

130 

117.5 

124 

7560 

14100 

190 

6 

146 

138 

117.5 

128 

4610 

14550 

116 

7 

138 

130 

117.5 

124 

7560 

14100 

190 

7R 

137 

129 

117.5 

123 

8220 

13950 

207 

8 

143.5 

135.5 

112.5 

124 

4110 

14100 

103.5 

8r 

141 

133 

112.5 

123 

4610 

13950 

116 

9 

138 

130 

107 

119 

4no 

13450 

103.5 

10 

129 

121 

103 

112 

5250 

12550 

132 

10R 

132 

124 

103 

114 

4510 

12800 

114 

Average  Re  =  14,630 
Average  Nu/Pr0#3  =  IO8.9 


Data  Sheet  23 
Run  no*.  C-9 
Flow  Rate  =  4.5  gpm 
Electrical  Power  78,900  BTU/hr 
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Q/A  =  94,800  BTU/hr-ftd 
Inlet  Pressure  =  48.5  psi 
Exit  Pressure  =  22  psi 
Inlet  Temperature  =  97°F 
Exit  Temperature  =  132°F 


AP  =  26.5  psi 


Thermo- 
couple 

T 
wo 

T  . 

•  Wl 

Tb 

T 

h 

Re 

Nr/Pu0'3 

2 

161 

153 

128 

141 

3790 

18300 

94.7 

2R 

160 

152 

128 

140 

3950 

18150 

98.7 

3 

154 

146 

123 

135 

4120 

17500 

103 

4 

147 

139 

119 

129 

4730 

16550 

108 

5 

138 

130 

114.5 

122 

6110 

15600 

153 

5R 

135 

127 

114.5 

121 

7570 

15400 

189 

6 

144 

136 

114.5 

125 

4410 

16000 

110 

7 

136 

128 

114.5 

121 

7010 

15400 

175 

7R 

135 

127 

114.5 

121 

756o 

15400 

189 

8 

141 

133 

110 

121 

4120 

15400 

103 

8R 

137 

129 

110 

120 

4990 

15300 

125 

9 

136 

128 

106 

117 

4310 

14800 

108 

10 

127 

119 

101 

110 

5260 

13900 

131 

10R 

129.5 

121.5 

101 

111 

4610 

14000 

115 

Average 

Re  = 

15,990 

Average 

Nu/Pr 

.0.3  = 

109.6 
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